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March  1.  1978 
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8:40  Workshop  Details  and  Overview 

8:55  Workshop  Objectives 

REQUIREMENTS  IN  AIRCRAFT  SURFACE  REPRESENTATION 

9:15  "Geometric  Modeling  in  Conceptual  Design" 

9:35  "Airplane  Surface  Descriptive  Geometry  Requirements" 

10:10  "Aircraft  Configuration  Modeling" 

10:30  "Geometric  Effect  on  Internal  Flow  Computations" 
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11:10  "Flow  Field  Grids" 

11:30  "Numerical  Aerodynamic  Simulation  Facility" 
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4:10  "A  Solution  to  the  Surface  Intersection  Problem" 

4:30  "Generation  of  Geometric  Input  for  30  Potential  Flow  Programs" 

4:50  Adjourn 
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8:30  "Body  Oriented  Mesh  Generation  for  3~0  Flow  Fields" 

8:50  "Boundary  Fitted  Coordinate  System  Using  Tension  Splines" 

9:10  "Tchebycheff  Approximations  for  Surface  Definition" 


I. .  Roberts/Ames 

J.  Ashbauqh/Ames 
T.  J.  Gregory /Ames 


P.  E.  Divan/Rockwell 

R.  Wallace/Boeing 

D.  R.  Clark/Analytical  Methods 

H.  Kao/Lewis 

C.  Yates/Langley 

J.  L.  Steger/Ames 

R.  L.  Carmichael /Ames 

W.  R.  Mann/Ames 

R.  E.  Miller/Coeing 


S.  M.  Staley/Univ.  of  Conn. 
D.  Anderson/Purdue 

T.  R.  Rau/Langley 

D.  P.  Roland/Informatics 
H.  Timmer/McDonnel 1- Doug las 
J.  L.  Hess  & D.  Halsey/ 
McDonnel 1- Doug las 


I.  C.  Bhately/General  Dynamics 
L.  L.  Pres ley /Ames 
R.  E.  Smith/Langley 
H.  Hoy /Ames 


Program  (Contd, ) 


2 


GEOMETRY  SYSTEMS  IN  USE 

10:50  "IPEGS-Interactive  Parametric  Equation  System?' 

11:10  "I CAD- Interactive  Computer  Aided  Design" 

11:30  "QUICK  GEOMETRY  Representation  of  Surfaces" 
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2:10  "Interfacing  Aerodynamic  Programs  to  AVID" 

2:50  Panel 

3:50  Summary 
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This  discussion  will  describe  the  Workshop  objectives, 
souse  of  the  reasons  for  holding  the  Workshop,  and  the 
emphasis  that  we'd  like  to  maintain  throughout  the 
program. 
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The  Workshop  origins  began  about  a year  ago  and  became  clearer  last 
September  in  a meeting  of  NASA,  Air  Force  and  Navy  representatives  in 
Washington.  At  this  meeting,  the  planning  for  the  new  aerodynamic 
paneling  technique  embodied  in  a system  called  FAN  AIR  (Paneling 
aerodynamics)  was  becoming  finalized.  Basically,  it  was  recognized  by 
those  involved  that  the  new  aerodynamics  techniques  could  use  more 
detailed  surface  definition  beyond  that  which  was  typically  used.  The 
next  slice  indicates  tlat  a realistic  and  complex  aircraft  configuration 
can  be  defined  in  terms  of  small  quadrilateral  panels  that  are  the  basis 
for  the  new  panelinc  techniques.  This  technique  is  based  on  the  fact 
that  pressure  on  each  panel  can  be  computed  as  a function  of  the  free 
stream  flow  conditions  and  the  inclination  of  the  panel  as  well  as  the 
i nc 1 i nation  of  the  surrounding  panels.  The  result  is  - large  linear 
algebra  problem  that  is  sclvatle  by  the  more  powerful  computers.  Referring 
back  to  the  last  slide,  it  was  apparent  that  significant  resources  were 
going  into  geometrV  definition  and  that,  in  fact,  more  resources  were  planned 
and  being  requested.  Also,  there  were  rany  alternatives  in  the  proposed 
approaches  and  the  extensions  and  enhancements  to  these  alternatives  provided 
a wide  variety  of  options.  At  this  point  the  picture  was  unclear  with  regard 
to  selecting  particular  options  and  it  was  felt  that  it  was  time  to  stop  and 
survey  the  whole  activity.  Hence,  the  interest  in  holding  this  Workshop. 


T.  J.  Gregory 
Vugraph  #2 


iMiBt IHBBaHK 


•he  specific  objectives  of  NASA  Hoaw 

would  help  coordinate  the  activity  that  that  “°rkSh°P 

between  aii  of  the  participants  i„™,ve, 

for  aerodynamic  computation  . surface  representation 

«-«nputation,  and  to  posslbiv  aa;n  = 

preferred  approaches  or  a tentative  on-  « “"«hsus  as  to 

between  the  approaches  if  ....  ° re9ardi"9  comnonality 

discussions  of 

Asaln.  we  wanted  industry  participatin,  Z'ZZJT"'''' 

Secondly,  „e  f.It  fte:  Workshop 

and  issues  surroundin9  the  techni  1 “"derstand  the  elements 

prepare  . reas^T  “ " “ ^ •«- 

of  surface  representation  techno^  ^ ^ ^ 


Jo  Gregory 
Vugraph 


ORIGINAL  PAGE  IS 
OF  POOR  QUALITY  ; 


The  Workshop  objectives  at  NASA  Ames  include  those  just  discussed 
and  additional  ones.  First,  our  primary  function  is  to  develop  new 
technology  in  aerodynamics  and  if  this  involves  geometry  or  surface 
representation  then  that  Is  an  area  of  interest  and  activity  for  us. 
We'd  also  like  to  have  the  aerodynamic  and  geometry  technology  used 
by  the  aerospace  community  and  that  means  having  their  opinions  and 
ideas  included  at  the  early  stages  of  this  development  process. 
Finally,  I think  the  key  motivation  for  the  Ames  aerodynmists  is  to 
have  a "painless"  geometry  package  which  v-<e  can  use  for  technology 
development  in  aerodynamics.  Certainly  if  the  Workshop  can  provide 
us  with  new  ideas,  direction,  and  approaches  that  would  lead  to 
acconpl ishing  this  latter  objective,  then  the  Workshop  would  be 
most  beneficial  for  us. 


T.  J.  Gregory 
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The  cent  slide  suggests  that  all  the  attendees  at  the  Workshop 
uould  benefit  from  a survey  of  the  field  and  would  provide  ar 
opportunity  to  show  the  capability  in  each  organization.  Both  of 
these  are  important  for  those  interested  in  participating  in 
NASA  sponsored  development  of  the  technology.  In  general 
probably  the  most  direct  benefit  for  all  attendees  at  the  Workshop 
will  be  to  gain  Information  or  even  software  that  may  be  of 
iimediate  value  to  their  own  efforts  in  aircraft  surface  definition. 
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Prior  to  starting  the  first  panel  session,  I'd  like  to  describe  some 
terminology  that  will  be  used  throughout  the  next  two  days  and 
emphasize  what  we  mean  by  aircraft  surface  representation.  There  are 
two  types  of  surfaces  of  interest,  aircraft  surfaces  and  surfaces 
within  the  flow  field.  Aircraft  surfaces  can  be  defined  by  a heirarchy 
of  elements.  The  first  element  is  a component  such  as  a wing,  body, 
nacelle,,  etc.  that  in  turn  can  be  described  by  surface  patches.  These 
are  described  by  either  systems  of  equations,  points,  or  functions. 

The  patches  can  be  further  subdivided  into  panels  as  indicated  on  the 
earlier  vugraph  of  a ca.^lete  aircraft  configuration.  These  panels  and 
the  patches  can  be  described  in  terms  of  the  edges  or  curves  along  the 
boundaries,  but  the  panels  are  sufficiently  described,  for  aerodynamic 

paneling  computations,  in  terms  of  points  at  the  intersections  of  their 
edges. 

Flow  field  surfaces  are  needed  to  define  such  items  as  shock  waves,  vortex 
sheets,  separation  bubble  areas,  etc.  These  are  of  major  importance  to 
the  aerodynamists  and  will  become  more  important  as  we  get  further  into 
this  technology.  The  same  heirarchy  of  elements  mentioned  above  can 
apply  to  flow  field  surfaces. 
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t There  is  another  field  of  major  importance  to  computational  aerodynamics 

and  that  is  the  definition  of  a flow  field  volume  (i.e.  solid).  This  is 
usually  done  by  means  of  grids  and  meshes  in  the  flow  field  that  are 
divided  by  either  uniform  spacing  in  the  simpliest  case,  or  by  stream- 
lines or  other  distributions.  These  are  used  to  make  finite  difference 
, computations  using  the  fundamental  partial  differential  equations  in 

"X  aerodynamics.  These  finite  differences  computations  are  an  emerging 

field  in  aerodynamics  and  generating  significant  interest  within  NASA. 
Again,  the  emphasis  in  this  Workshop  is  on  surface  definition  and  not 
on  flow  field  volume  (solid)  definition.  Perhaps  the  specialists  in 
this  latter  field  will  generate  a workshop  in  the  future. 
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To  reiterate,  the  Workshop  emphasis  is  on  surface  representation 
and  ,ts  integration  with  aerodynamics.  computers,  graphics  and 
Wind  tunnel  model  fabrication  as  well  as  flow  field  grid  generation, 
but  none  of  these  Items  per  se.  It's  our  intention  to  try  and 
focus  the  oiscussions  and  papers  at  the  Workshop  on  surface 

representation  and  to  defer  detailed  discussions  of  these  other 
items  to  other  workshops  or  conferences. 
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The 'geometric  description  of  an  airplane  usually  starts  with  a character- 
ization of  the  configuration.  There  are  only  a few  elements  of  a con- 
figuration, such  as  the  payload,  number  of  crew,  number  of  engines,  engine 
location,  which  are  fixed.  The  remainder  of  the  configuration  character- 
istics, including  fuselage  geometry,  wing  geometry,  control  surfaces,  . 
empennage  geometry,  landing  gear  arrangement,  wing  high-lift  devices,  etc., 
are  all  variables  that  get  perturbed  during  design  evolution.'  Therefore,  the 
first  requirement  on  any  airplane  configuration  geometry  description  system  is 
flexibility.  Easy  modifications  to  geometry  of  a component , or  the  relations 
of  components  with  respect  to  each  other  Is  essential i 
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This  result  of  a preliminary  design  exercise  for  a transonic  airplane  con- 
figuration shows  many  of  the  elements  described  by  the  prior  charts.  A 
transonic  airplane,  typified  by  this  NASA  contract  study  configuration,  has 
the  added  complex  requirement  of  satisfying  an  area  rule  distribution. 
Therefore,  any  relatively  minor  modification  such  as  wing  area,  nacelle 
placement,  etc.  has  major  significance  to  the  total  design  compatability. 
Simultaneous  satisfaction  of  all  design  requirements  and  performance  criteria 
demands  very  sophisticated  and  efficient  geometry  processing.  This  figure 
illustrates  many  of  the  variables  that  have  a first  order  constraining  effect 
on  airplane  geometry  description. 
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As  indicated  earlier  by  the  airplane  design  block  diagram,  one  of  the  geometry 
requirements  on  the  complete  configuration  are  the  pitch  and  roll  ground 
clearances  as  shown  by  this  line  drawing  of  the  747.  These  geometric  con- 
ditions are  necessary  for  flight  control  evaluations,  damage  determination 
studies,  and  pilot  vision  polars.  Again,  if  any  of  these  critical  conditions 
are  not  met  by  the  airplane  configuration,  then  the  geometry  must  be  changed. 
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During  configuration  design  or  analysis,  it  is  necessary  to  perform  extensive 
analytical  and  experimental  studies.  Ideally,  the  same  geometry  Is  used  for 
generating  mathematical  representation  of  an  airplane  configuration  and  its 
corresponding  flow  fields  as  is  used  to  generate  wind  tunnel  models  for 
testing.  This  figure  shows  one  kind  of  mathematical  representation  of  an 
engine  nacelle  and  its  associated  exhaust  flow.  This  representation  is  typ- 
ical of  analyses  used  to  determine  nacelle  interference  with  wing  air  flow 
properties. 
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PANEL  AERODYNAMICS  SINGULARITY  REPRESENTATION 

OF  DUAL  BARREL  FLOW  NACELLE 
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Once  the  nacelle  geometry  details  have  been  adequately  modeled,  with  suitably 
dense  singularity  arrays,  then  there  is  a gradual  buildup  to  the  complete  con- 
figuration. This  figure  shows  the  addition  of  the  nacelle  strut  and  adjacent 
wing,  to  the  nacelle.  This  three-view  layout  Illustrates  the  complex  local 
tailoring  associating  with  a close  coupled  nacelle.  The  nacelle  is  canted 
with  respect  to  the  strut  and  the  wing  and  the  strut  is  tailored  to  minimize 
local  interference.  This  complexity  is  readily  seen  in  these  orthogonal  views. 
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[ " A perspective  view  of  the  wing  segment,  nacelle  strut,  nacelle  and  its  asso- 

F ciated  flow  control  surfaces  is  the  best  way  to  appreciate  its  geometric 

intricacy.  Notice  that  the  exhaust  flow  tube  is  controlled  well  past  the  wing 
trailing  edge  and  there  is  an  intake  flow  control  surface  located  within  the 
nacelle  under  the  leading  edge  of  the  nacelle  strut.  Ill  of  these  singularity 
' surfaces  are  necessary  to  accurately  evaluate  the  engine  installation  inter- 
im ference  flow  phenomena. 


s 


I 


I 


I 


f 


14*  UkiLtce. 


original  page 
OF  POOR  QUALITY 


This  two  view  drawing  shows  the  complete  model  used  for  the  study.  The  body 
was-  simplified,  because  of  its  remote  distance  from  the  region  of  interest. 
There  are  approximately  2700  singularities  associated  with  the  analytical 
model,  which  produces  a like  number  of  simultaneous  equations  for  solution. 
Therefore,  the  computing  cost  of  such  a complex  mathematical  problem  requires 
very  accurate  geometrical  modeling  with  adequate  visibility  of  that  geometry 
prior  to  committing  to  the  computing  process.  One  set  of  computed , solutions 
can  cost  as  much  as  many  wind  tunnel  model  test  runs. 
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Before  describing  geometry  system  requirements  as  they  are  viewed  today,  it  is 
well  to  understand  current  lofting  practice  of  The  Boeing  Commerical  Airplane 
Company.  This  wing  perspective  drawing  shows  the  lofting  views  constructed 
for  generating  a stored  definition  in  two  views;  plan  view  and  rear  view. 

This  process,  of  course,  has  its  heritage  in  the  ship  building  industry.  Two 
types  of  curves  have  been  prevalent  in  this  current  practice:  (1)  conic 
chains,  and  (2)  cubic  chains.  The  conic  chain  has  been  the  traditional  fav- 
orite, because  it  is  simple  to  generate  by  most  mechanical  drawing  processes 
and  to  check  by  most  manual  computing  methods.  It  is  also  easy  to  control 
inflection  points.  This  curve  type  is  used  for  the  master  definition  of  pro- 
duction configurations. 

The  cubic  chain  with  its  point,  slope,  and  curvature  continuity  at  given  nodes 
has  usually  been  best  suited  for  design  and  development  purposes.  It  is 
analagous  to  the  process  of  using  ducks  and  a spline  on  the  drawing  board.  It 
is  practical  to  embellish  this  mathematical  representation  with  both  point 
enrichment  and  smoothing  processes.  However,  cubics  have  the  disadvantage  of 
causing  ripples  or  inflection  points  for  data  sets  that  are  not  smooth 
initially.  These  traits  make  more  complex  algorithms  necessary  for  generating 
good  cubic  chain  curves. 
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CURRENT  LOFTING  PRACTICE 


CONIC  CHAIN  CURVE 


• BEST  FIT  TO  POINTS  WITH 
SLOPE  CONTINUITY 

• MASTER  DEFINITION 
CUBIC  CHAIN  CURVE 


• POINT  AND  SLOPE 
CONTINUITY 

• DESIGN  AND  DEVELOPMENT 


2 1 UJalW 


Current  practices  using  conic  and  cubic  curves  present  significant  limitations 
for  the  geometry  description  of  an  airplane  conf iguration  and  its  components. 
Multiple  curve  types  cause  difficulties  with  geometry  automation  and  for  the 
design  user.  The  single-valued  surfaces  associated  with  current  lofting  prac- 
tices are  usually  imposed  by  the  extraction  process  which  cannot  adequately 
distinguish  between  multivalued  components.  Another  difficulty  encountered  is 
mating  surfaces  of  various  components,  because  the  surfaces  are  defined  as 
projected  control  curves.  This  representation  makes  the  definition  (stored  as 
equation  coefficients)  very  expensive  to  uniquely  transform  between  skewed 
coordinate  systems. 


All  'of  these  cited  factors  complicate  the  design  and  analysis  processes, 
because  the>  are  not  very  flexible  and  are  costly  to  use.  In  many  cases,  the 
systems  are  designed  mostly  for  geometry  extraction  and  provide  little  flex- 
ibility for  geometry  generation  processes,  which  often  require  data  enrichment 
and  extensive  three-dimensional  smoothing.  „ 


LIMITATIONS  OF  CURRENT  PRACTICES 

MULTIPLE  CURVE  TYPES 

• COMPLICATES  SMOOTHING 

• COMPLEX  DEFINITION  STORAGE 

• DIFFICULT  DATA  EXTRACTION 

SINGLE-VALUED  SURFACES 

• D'FFICULT  JOINT  SMOOTHING 
•COMPLICATES  MULTI-VALUED  COMPONENTS 

SURFACE  MATING 

•DIFFERENT  AXIS  SYSTEMS 

• PROJECTED  CONTROL  CURVES 

LIMITED  TRANSFORMATIONS 

•NO  DIRECT  TRANSFER  BETWEEN  SKEWED  AXES 

• REFIT  EXTRACTED  POINTS 

DESIGN  AND  ANALYSIS 

•DESIGN  REQUIRES  ENRICHING  AND  SMOOTHING 
•ANALYSIS  DATA  EXPENSIVE  TO  EXTRACT 


Basic  geometry  system  requirements  contain  two  principal  factors,  the  first  of 
which  is  mathematical.  The  mathematical  factors  directly  lead  to  the  con- 
clusion that  surface  representation  is  best  accomplished  with  parametric, 
biquintic  patches.  This  makes  mathematical  practice  reasonably  consistent 
with  past  cubic  concepts,  except  that  the  higher  order  polynomial  provides  the 
essential  element  of  local  character.  The  parametric  form  provides  the  nec- 
essary capability  for  handling  multivalued  surfaces  and  performing  smoothing 
and  extraction  processes  with  more  consistency,  since  quintic  equations  are 
used  over  an  entire  surface.  Where  simple  curve  types,  such  as  straight  lines 
and  circles  become  necessary,  the  quintic  equations  are  perfectly  adequate  for 
defining  these  surface  regions  to  well  within  data  extraction  computation 
accuracy. 


BASIC  GEOMETRY  SYSTEM  REQUIREMENTS 


MATHEMATICAL  FACTORS  ARE- 


ll)  ONE  CURVE  TYPE  PATCHES 

(2)  CONTINUOUS  CURVATURE  SURFACES  • 

(3)  LOCAL  CHARACTER  PATCHES 

4 

(4)  REAL,  MULTIVALUED  SURFACES 

(5)  ENRICHING,  SMOOTHING  & EXTRACTING 


BIQUINTIC 

. . V 

1 PARAMETRIC  FORM 


2 S lUfiiace 
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The  second  major  factor  associated  with  basic  geometry  system  requirements  is 
the  user  working  environment.  The  two  primary  innovations  that  have  recently 
improved  the  geometry  working  environment  are  the  minicomputer  and  its  asso- 
ciated micro  processors  for  use  in  interactive  graphics  devices.  With  th4s 
hardware  capability  it  is  practical  to  provide  excellent  accuracy  with  low 
computing  costs  in  a geometry  system  that  can  be  used  all  the  way  from  pre- 
liminary design  through  to  detailed  design  activities.  It  is  essential  for 
all  engineering  technologies  to  have  access  to  and  influence  the  design  geo- 
metry evolution  to  adequately  reflect  their  responsibilities.  Conversational 
interactive  graphics  gives  the  average  user  a reasonably  acceptable  working 
environment  which  will  not  overwhelm  him  with  the  necessity  for  training  that 
makes  him  a computer  expert.  The  last  environmental  requirement  is  commun- 
icating geometry  with  accuracy  and  speed  between  all  involved  developers  as 
well  as  the  ultimate  users,  who  are  charge  with  building  the  airplane. 


BASIC  GEOMETRY  SYSTEM  REQUIREMENTS 


ENVIRONMENTAL  FACTORS  ARE  - 

(1)  GOOD  ACCURACY  & LOW  COMPUTING  COST 

12)  ONE  SYSTEM  FOR  PRELIMINARY  & DETAILED  DESIGN 

13)  EQUALLY  USEFUL  FOR  ALL  ENGINEERING  TECHNOLOGIES 
M)  USE  CONVERSATIONAL  INTERACTIVE  GRAPHICS 

(5)  AUTOMATE  GEOMETRY  DATA  BASE  INTERFACES 


This  figure  shows  the  implications  of  smoothing  on  analytical  results.  A 
simple  exercise  was  performed  to  manually  record  the  coordinates  of  an  airfoil 
to  compute  its  pressure  distribution.  As  seen  from  the  pressure  coefficient 
graph,  the  manual  unsmoothed  data  caused  severe  adverse  pressure  gradients 
near  the  nose  of  the  airfoil,  which  would  probably  lead  to  local  adverse 
effects  on  the  boundary  layer,  if  not  separation.  Similarly,  at  the  trailing 
edge  there  was  an  added  adverse  gradient  due  to  the  data  irregularity.  Simple 
two-dimensional  analytical  smoothing  produced  the  smoother  more  satisfactory 
dashed  line  results.  In  all  honesty,  this  was  not  a rigged  case,  but  simply 
an  illustration  of  an  everyday  event  when  accomplished  without  proper  atten- 
tion to  geometric  properties.  The  results  cause  poor  aerodynamic  performance. 
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PRESSURE  DISTRIBUTIONS  COMPUTED  FROM  MANUAL 

& SMOOTHED  ANALYTICAL  DATA 


MANUAL,  UNSMOOTHED  DATA  INPUT 
4 ANALYTICAL,  SMOOTH  DATA  INPUT 


PRESSURE 

COEFFICIENT 


y/c(ANALYTICAL) 

0.02110 

0.03070 

0.03690 

0.04190 

0.04905 

0.05340 

0.05647 


0.4  , 0.6 

x/c 

CHORD  FRACTION 
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The  'two-parameter  biquintic  surface  representation  takes  the  equation  form 
shown  on  this  chart.  Using  the  previous  wing  illustration,  the  patch  shown  m 
real  geometry  is  mathematically  handled  in  parametric  form  and  the  definition 
is  stored  as  derivatives  at  the  corner  points  rather  than  as  coefficients  of 
the  respective  bounding  quintic  lines.  There  are  many  advantages  to  t^iis  type 
of  representation  when  put  into  practice  as  computerized  methodology. 


TWO-PARAMETER  SURFACES 
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y ■ y (s.  t)  WHERE  xls.  t)  • E Ex-.sV 
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PARAMETRIC  B IQUINTIC  EQUATIONS 


PATCH  REAL  GEOMETRY 
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AT  , 
PAT 
COW 
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PATCH  PARAMETRIC  REPRESENTATION 
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To  illustrate  the  data  for  representing  a typical  airplane  configuration  these 
data  array  schematics  are  shown.  At  the  juncture  of  every  pair  of  lines  there 
is  a data  set  corresponding  to  that  corner  of  the  patch.  Note  that  most  of 
the  data  arrays  are  rectangular.  Except  for  the  cutout  regions  where  spatial 
fairing  properties  are  necessary,  the  configuration  paneling  is  straight- 
forward. Where  intersections  cause  local  need  for  fairing  or  irregular  boun- 
daries, then  spatial  handling  techniques  are  required. 


Developing  the  paneling  representations  of  a configuration  without  consid- 
erable automation,  is  a very  tedious  and  time-consuming  task.  It  is  not 
uncommon  for  an  engineer  to  expend  a man  month  in  developing  such  extensive 
paneling  schemes.  Usually  it  is  necessary  to  build  these  representations  in 
an  component-by-component  fashion.  By  running  simple  evaluations  of  isolated 
components,  it  then  becomes  possible  to  develop  confidence  that  the  final 
results  sought  will  be  computed  with  good  accuracy.  It  is  very  common  to  find 
that  people  have  been  unable  to  perform  a satisfactory  analytical  evaluation 
of  a configuration,  simply  because  there  was  insufficient  time  to  develop  the 
geometry  and  its  associated  paneling  scheme. 
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SUMMARY 


1 - BASIC  REQUIREMENTS  TOUGH  TO  SATISFY 

• ONE  CURVE  TYPE  WITH  LOCAL  CHARACTER 

• MULTIVALUED  , REAL  SURFACES 

• ENRICHING  & SMOOTHING  WITHOUT  DISTORTION 

• VERY  ACCURATE , YET  EASY  TO  CONTROL  & CHEAP  TO  USE 

2 - TECHNOLOGY  IS  AVAILABLE  FOR  SYSTEM  DEVELOPMENT 

• MATHEMATICS  OF  BIQUINTICS  & TOPOLOGY 

• MINICOMPUTER  REFRESH  GRAPHICS 

• LARGE  SCIENTIFIC  COMPUTERS 

• EXPERIENCE  WITH  LESS  CAPABLE  GEOMETRY  SYSTEMS 
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INPUT  OF  THE  CONFIGURATION 

PANEL  BY  PANEL  CORNER  POINTS 


ACCEPTABILITY 

UNACCEPTABLE 


SECTION  BY  SECTION 


BLOCK  BY  BLOCK 


EACH  POINT  DEFINED. 

MACHINE  CONNECTS  ADJACENT  FAIR  1 

SECTIONS  TO  MAKE  PANELS 

DEFINE  CURVE.  MACHINE 

DIVIDES  AND  CONNECTS.  G00D  ' 

DEFINE  END  CURVES. 

MACHINE  DIVIDES  AND  600D  1 

CONNECTS  FOR  PANELS. 


***  DATA  PREVIEW  IS  ESSENTIAL  FOR  EFFECTIVE  MACHINE  USE  *** 


1.  USE  OF  A MACHINt  DIGITISER  RAISES  THESE  TO  GOOD  AND  EXCELLENT. 
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GEOMETERY 

+ AERODYNAMICS  ( AIRLOADS, VELOCITIES  ETC,  ) 

+ STREAMLINE  BEHAVIOR 
+ BOUlIUARYLAYER  behavior 
+ ETC. 

= LARGE  VOLUMES  OF  PRINTOUT* 

CAPABILITY  TO  REVIEW  OUTPUT  DATA  INTERACTIVELY  IS  ESSENTIAL. 
MACHINE  PLOT  AND  PRINT  ONLY  REQUIRED  DATA.  ( SAVE  THE  REST  ) 


MOST  OF  WHICH  NEVER  GETS  READ. 


OUTPUT 


/)M1 


THE  USE  OF  INTERACTIVE  COMPUTING  WITH  GRAPHICS  CAN  DRAMATICALLY 

SPEED  UP  THE  DATA  FLOW 


DATA  OUT 


DATA  OUT 


I = INTERACTIVE;  B = BATCH 


AIRCRAFT  CONFIGURATION  MODELING 


FROM  TERMINAL  INTERACTIVE 


OUTPUT  RESTART 


THIS  PROGRAM  PLOTS  SHAPES  FROM  INPUT  DATA  DECK 

OF  THE 

ttt  UBAERO  PROGRAM  *** 


INPUT  CAN  BE  FULL  UBAERO  DATA  SET  OR  PARTIAL, STARTING  AT  CARD  7 
KEY  IN  BAUD  RATE  AND  RETURN 

? 300 


YOU  HAUE 
1 
2 
3 
A 
5 


THE  FOLLOWING  PLOT  OPTIONS 
PLOT  ALL  SECTIONS 

PLOT  ALL  SECTIONS  IN  A DESIGNATED  BLOCK 
PLOT  ANY  DESIGNATED  SINGLE  SECTION 
COMPARE  UP  TO  10  DESIGNATED  SECTIONS 
TERMINATE! 


KEY  IN  DESIRED  OPTION  AND  RETURN! 

? 2 

WHAT  BLOCK  DO  YOU  UISH  TO  PLOT? 
KEY  IN  BLOCK  NUMBER  AND  RETURN* 


r 


yi 


THIS  PROGRAM  PLOTS  THE  OUTPUT  DATA  FROM 


THE 

UBAERO  FAMILY  OF  PROGRAMS 


YOU  HAUE  THE  FOLLOUIHG  OPTIONS. 

1 GEOMETERY  GROUP 

2 AERODYNAMICS  GROUP 

3 STPEAMLINE  GROUP 

4 BOUNDAPYLAYER  GROUP 

KEY  IN  SELECTED  OPTION  AND  RETURN* 

<eg.3RETURN) 


YOU  HAVE  SELECTED  PLOTS  FROM  THE  STREAMLINE  GROUP 

********** 

SELECT  PARAMETERS  TO  BE  PLOTTED  FROM  THE  MENU  BELOU 


1. 

X 

S.  ux 

0 

u • 

Y 

6.  UY 

3 . 

Z 

7.  UZ 

4. 

s 

3.  U RESULTANT 

9.  CP 

10.  K1 

11.  K2 

12.  K3 


NOTE 


VOU  CAN  PLOT  ANY  TIJO  PARAMETERS  AGAINST 


EACH  other 


KEY  IN  PARAMETERS  AND  RETURN!  (eg.  4,9, 1#3 

THERE  ARE  **  9 * t STREAMLINES.  DO  YOU  UISH  TO  PLOT  THEM  ALL7 
KEY  IN  YES  OR  NO  AND  RETURN! YES 


THIS  PROGRAM  GENERATES  TEKTRONIX  PLOTS  OF  GENERAL  3D  BODIES 

INPUT  NAY  BE  READ  FROM  UBAERO  TAPE11  OR  UIA  THE  KEY  BOARD 

TX***XXXX* 


SELECT  INPUT  NODE 

USING  TAPE11?  KEY  IN  1 AND  RETURN 
USING  KEYBRD?  KEY  IN  2 AND  RETURN, 

***NOTE  t*  IF  1 SELECTED  TAPE/FILE  MUST  BE  IDENTIFIED  TX* 


? 1 

DO  YOU  UANT  TO  DPAU  THE  UHOLE  BODY? 
KEY  IN  YES  OR  NO  AND  RETURN! 

? no 


fi£TI0N  0F  *0DY  T0  BE  PLOTTED  MUST  HAUE  A CONNECTED 
STRING  OF  PANEL  INDICES. 

KEY  IN  INDICES  OF  FIRST  AND  LAST  PANELS  AND  RETURN1 
(EG.  172,431  RETURN') 

? 22,44 


* if 


KEY  IN  NUU, IPRINT, IHIDE, I BUG 

NUU  • NUNBER  OF  UIEUS  **  NO  HORE  THAN  10  * 
IPRINT  « 1 PRINTS  INPUT  DATA 

- 0 NO  DATA  PRINT  OUT 
IHIDE  = 0 ELLIMINATES  HIDDEN  LINES 
« 1 LEAUES  HIDDEN  LINES 
IBUG  * 0 SUPPRESSES  DEBUG  PRINT  OUT 
» 1 FULL  DEBUG  PRINT  OUT 


? 2,0, 1,0 


IS  BODY  SYP1ETRICAL  ABOUT  X-Z  PLANE? 


KEY  IN  YES  OR  NO  AND  RETURN! 


9 yes 


DO  YOU  UANT  TO  PLOT  REFLECTED  BODY? 


KEY  IN  YES  OR  NO  AND  RETURN! 


? no 


key  in  XUUE.YUUE.ZUUE  - ORDINATES  OF  UIEU  POINT 
***  POINT  MUST  LIE  OUTSIDE  BODY  *** 


O BO  105  BOTTOM  CENTERLINE 


(a)  Top  Centerline 

Figure  19.  Comparison  of  Calculated  and  Experimental  Pressure  Distributions  on  the  E0105 
Fuselage. 


DRAG  PROGRAM 


Model  2.  Vokticity  Panels 


DOWNSTREAM  PARTS  OF  "SPLIT"  PANELS  HAVE 
PANELS  HAVE  LINEAR  VORTICITY  ZER0  V°RTICITY  VALUE 

DISTRIBUTION  AS  WELL  AS  , 

UNIFORM  SOURCE 


Figure  1: 


The  method  chosen  for  computation  generally  dictates  the  requirement 
for  geometric  smoothness.  In  the  case  of  viscous  n«  coloration 

" “ aXiSy“i,etriC  diffus<--*-  *«*  rtth  n center., ody,  one  of  the 
approaches  is  to  parabolise  the  Xavier -stoles  options  in  . 
streamline  orthogonal  coordinate  system  and  then  to  perform 
computation  by  marching  In  this  coordinate  system  « observe 
that,  as  sheen  in  item  (e) , to  insure  stable  computation  the 

Streandine  coordinates  have  to  be  accurately  determined  up  to 
third  derivatives. 


VISCOUS  FLOW  IN  A DIFFUSER  DUCT 

(Ol)  Transjorm  Navier-  Stokes  FjuaNons  +o 
Streamline  Orthogonal  Coordinate  System 


(b)  “Parabolize  the  Resulting  Fauations  -for  Marching 
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Figures  2,  3 and  U 


Figure  2 shows  the  mesh  distribution  c-f  a smoothed  annular 
diffuser  duct  and  the  calculated  velocity  profiles  at  each 
station.  The  computation  was  terminated  normally  at  the  end 
of  the  duct.  However,  if  would  use  an  identical  computer 
program  but  with  an  unsmoo*  ed  geometry,  computation  would 
become  unstable  and  brc«_h  „uwn  shoi'tly  after  the  entrance 
section.  The  configuration  shown  in  Figure  'I  is  a partially 
smoothed  one,  but  the  computation  still  aborted  before  it  reached 
the  exit  section.  (An  indication  of  irregularity  may  be  seen  at 
the  waist  of  the  centerbody.) 


Figure  5: 

All  three  segments  of  (1) , (2)  and  (3)  are  the  traces  of  a cubic 

eejuation,  and  they  all  satisfy  the  tangcncy  continuity  condition 

A and  R as  shouTi  (that  is,  their  unit  tangent  vectors  at  A 

and  B match  with  the  unit  tangent  vectors  of  the  neigliboring 

segments)  . The  striking  difference  in  appearance  is  solely  d\je 

to  the  difference  in  magnitude  of  the  tangent  vectors  in  terms 

of  u.  Thus,  when  a general  parameter  other  than  the  arc  length 

is  used  for  curve  fitting,  the  first  derivatives  on  both  sides 

of  a data  point  need  not  equal.  This  property  also  appears  in 

second  and  higher  derivatives.  (The  smoothness  requirement  for 

second  derivatives  is  to  match  the  curvature  on  both  sides  of 

d2 

a data  point,  which  is  1 , with  s being  the  arc  length.) 

ds 

Since  the  degenerated  case  of  a Coons’  surface  patch  equation 
(say,  w=o)  is  a cubic  or  a fifth  degree  polynomial,  the  property 
of  discontinuity  in  parametric  derivatives  of  u and  w also  pre- 
vails in  Coons’  boundary  matrix. 
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Fifth  degree  surface  patch 
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Fi<5,  G basic  equations  for  Coons’ 
5th  degree  surface  patch 


GEOMETRY  REQUIREMENTS  FOR  UNSTEADY  AERODYNAMICS 


E.  Carson  Yates,  Jr.  and  Luigi  Mori  no 
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Accurate  calculation  of  aeroelastic  characteristics  required  for  the  analysis  and  design  of  high- 
performance  aircraft  requires  accurate  and  efficient  evaluation  of  steady  and  unsteady  aerodynamic  loads 
on  aircraft  having  arbitrary  shapes  and  motions,  including  structural  deformations.  This  presentation 
will  address  the  aircraft  geometry  requirements  for  unsteady  aerodynamic  computations  and  will  emphasize 
differences  between  requirements  for  steady  and  unsteady  flow. 


* # 


Requirements  for  aeroelastic  analysis  and  design  are  in  several  respects  ore  complicated  and  more 
severe  than  for  the  more  conventional  steady-state  aerodynamics . For  example:  (1)  The  aeroelastician 

deals  with  flexible  structures  so  that  even  in  steady-state  conditions,  the  aerodynamic  load  is  a function 
of  structural  deformation,  and  vice  versa.  (2)  The  unsteady  aerodynamic  formulations  required  in  dynamic 
aeroelasticity  involve  complex  quantities  (e.g.,  normalwash  velocities,  aerodynamic  influence  functions, 
and  pressure)  that  manifest  time-  or  frequency-dependent  attenuations  and  phase  shifts  relative  to  steady 
state.  (3)  In  dynamic  a<  roe  1 nsti city — flutter,  for  example — the  aeroelastician  must  evaluate  pressure 
distributions  for  vibration  mode  shapes  that  are  much  more  wiggly  than  a typical  steady-state  mean-camber 
surface.  The  corresponding  pressure  distributions  will  also  be  more  wiggly  than  those  for  steady  state 
so  that  computational  convergence  requirements  are  usually  more  severe  than  for  steady  state.  (4)  Flutter 
analyses,  as  well  as  iterative  structural  resizing,  require  evaluation  of  pressure  distributions  for  a 
multiplicity  of  mode  shapes,  frequencies,  aircraft  loading  conditio,.-,  etc.  Consequently,  computational 
efficiency  is  vital,  and  it  is  essential  to  minimize  the  amount  of  recomputation  required  when  mode  shapes 
and/or  frequencies  are  changed. 

Kith  these  thoughts  in  mind,  we  shall  discuss  geometry  requirements  within  the  framework  of  the 
SOUSSA  aerodynamic  formulation  because  it  is  the  most  general  potential-flow  program  that  we  now  have 
under  development  (with  regard  to  aircraft  geometry,  motion  and  deformations,  and  speed  ranges)  and 
because  present  and  future  SOUSSA  geometry  requirements  are  as  stringent  as  those  for  anv  aerodynamic 
program  that  we  now  anticipate.  Geometry  required  is  considered  to  be  comnesed  of  three  parts:  (1) 

shape  of  vehicle,  (2)  orientation  of  vehicle,  (3)  deformat  ion (s) . Orientation  involves  little  more 
than  a rotation  of  coordinate  axes  and  consequently  will  not  be  emphasized  here.  Defcrmaticns  can  be 
finite  but  are  more  usually  taken  to  lie  infinitesimal  .and  .approximated  by  a linear  combination  of  the 
natural  undamped  vibration  modes  of  the  aircraft.  As  m.anv  ns  two  dozen  modes  or  more  may  be  required 
to  converge  the  aeroelastic  solution.  A corollary  geometrical  requirement  is  determination  of  wake 
shape  which  is  not  known  a priori  although  it  may  be  assumed  to  be  flat  for  many  applications.* 


UNSTEADY  COMPLICATIONS: 

o load  is  function  of  deformation  and  vice  versa 

o COMPLEX  QUANTITIES 
o WIGGLY  DEFORMATION  MODES 

0 multiplicity  of  modes,  frequencies,  etc. 

GEOMETRY  REQOIRED: 
o SHAPE 
0 ORIENTATION 
o DEFORMATI ON( S) 


To  set  the  stage,  a brief  review  of  SOUSSA  formulation  is  in  order.  Application  of  Green's  theorem 
leads  to  an  integral  equation  for  the  perturbation  velocity  potential  at  any  point  P in  the  flow 

or  on  the  flow  boundary  (i.e.,  on  S).  Note  that  the  second  integral  contains  only  linear  terms  which 
are  integrated  over  the  boundary  surface-  S,  whereas  the  first  integral  contains  nonlinear  terms  F, 
involving  products  of  derivatives  qf  which  must  be  Integrated  over  the  fluid  volume. 


The  boundury  condition  clearly  shows  the  effect  oi  time  variation  of  S. 


harmonic,  for  example,  the 


ds 

It 


term  becomes  iuiS  so  that  the  normalwash  at 


it  the  variation  is 
S - 0 becomes  complex. 


The  imaginary  part,  however,  involves  only  surface  ordinates  (including  displacements  and  deformations), 
whereas  the  real  (steady-state)  part  involves  derivatives  of  surface  ordinates.  Thus,  introduction 
of  unsteadiness  does  not  i upose  more  stringent  requirements  on  surface  definition  as  far  as  quantities 
required  are  concerned.  It  .;ay,  however,  require  greater  accuracy  and  greater  amounts  of  geometrical 
information  (e.g.,  for  more  points  on  the  surface)  in  order  to  define  adequately  the  wiggly  modes  of  defoi 
mation  referred  to  previously. 


The  influence  of  nonlinear  terms  F is  being  studied  in  the  development  of  SOUSSA  aerodynamics  for 
the  transonic  range.  However,  these  terms  are  not  included  in  the  present  computer  program. 
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STEADY,  OSCILLATORY  AND  UNSTEADY  SUBSONIC  AND  SUPERSONIC  AERODYNAMICS 
• • (SOUSSA) 


Objective:  An  accurate,  general,  unified- method  for  calculating  steady  and  unsteady  loads  on 
complete  aircraft  with  arbitrary  shape  and  motion  in  subsonic  or  supersonic  flow, 
with  emphasis  on  application  in  computer-aided  structural  design 

Approach:  Green's  theorem  is  used  to  formulate  exact  integral  equation  for  potential. 

• * 

f(p‘V  = jfffo  Fd^/T,  + fff[v,s{&  vj-  «>?,&)  -4  jfc  11  - f^)]Jas\'  dSdt, 

where  9 * = perturbation  velocity  potential 

G = Green's  function 

' F = nonlinear  terms  • % 

S(x.y.z.t)  = 0 defines  body  surface 

|os|  = + s/  + s/ 

Exact  boundary  condition  on  body  is 

§j=  ?df7$‘vs=U+(/~(?+w)-vs=o 
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Surface  paneling  and  Laplace  transform  solution  yield  a matrix  equation  relating  the  unknown  potential 
f on  the  vehicle  surface  to  the  normalwash  ip.  Elements  of  the  coefficient  matrices  are  Independent 
of  normalwash  (and  hence  deformation)  and  are  simple  functions  of  the  Laplace  variable  s.  For  a given 
paneling  arrangement  they  depend  only  on  Mach  number. 

Use  of  arbitrary  nonplanar  quadrilateral  panels  permits  matching  nodes  of  the  aerodynamic  panels  to 
the  nodes  of  a structural  finite-element  model,  if  desired,  in  order  to  use  the  nodal  coordinates  and 
calculated  displacements  directly  without  requirement  for  Interpolation.  In  general,  however,  solution 
for  the  velocity-potential  matrix  requires  the  following  geometrical  input:  (1)  Coordinates  of  panel 

nodes  usually  obtained  by  interpolation  (lofting)  from  aircraft  shape  information.  (2)  Time-dependent 
normalwash  at  control  points  which  usually  do  not  coincide  with  panel  nodes.  Normalwash  Involves 
coordinates  and  slopes  obtained  from  aircraft  shape  plus  orientation  plus  deformation.  Note  that 
increasing  the  number  of  deformation  modes  used  involves  only  adding  columns  to  the  4)  and  f matrices, 
and  updating  the  entire  set  of  deformation  modes,  as  in  a structural  design  application,  involves  only 
replacing  the  matrix.  The  Y and  Z matrices  are  unaffected  in  either  case. 

Surface  pressures  are  obtained  from  Bernoulli's  equation.  Generalized  aerodynamic  forces,  including 
aerodynamic  coefficients  and  stability  derivatives,  are  determined  from  weighted  integrals  of  the 
pressure  which  require  values  of  surface  displacement  (due  to  rigid-body  rotation  and/or  modal  deformation) 
at  a set  of  integration  points  which  may  not  coincide  with  the  panel  nodes  nor  normalwash  control  points. 

The  geometrical  information  required  by  SOUSSA  can,  of  course,  be  generated  with  any  suitable 
geometry  preprocessor  as  long  as  the  results  are  cast  in  required  SOUSSA  input  format.  It  is  evident, 
however,  that  automatic  paneling  capability  is  essential  to  the  efficient  processing  of  complicated 
shapes  and  deformations  that  may  require  many  hundreds  of  panels.  Such  capability  should  include  not 
only  automatic  calculation  of  the  coordinates  of  nodes,  normalwash  control  points,  and  integration  points, 
but  also  automatic  identification  numbering  for  these  points  as  well  as  for  the  panels  and  systematic 
identification  of  which  nodes  go  with  which  panels. 


Solution  by  spatial  discretization  with  arbitrary  nonplanar  quadrilateral 
sur  ? ce  panels  and  time  solution  by  Laplace  transform  results  in 

where  (pj,  '=  Laplace  transform  of  perturbation  velocity  potential 
% = Laplace  transform  of  normalwdSli 

" SjA  ~(<VA 

s = Laplace  transform  variable 

Bjh»  cjh>  Djh>  Fjn’  Gjn  “ integrals  over  surface  panels,  independent 

of  normalwash  and  s 

0jh>  Tljh  " la9  ‘functions 
Snh  = +-1 


Surface  pressures  are  obtained  from  Bernoulli's  equation. 
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GENERAL  POTENTIAL-FLOW  AERODYNAMICS 

(SOUSSA) 

I 

GENERAL  FINITE-ELEMENT  METHOD: 

o ARBITRARY  COMPLETE  A/C  CONFIGURATION 

© STEADY  AND  GENERAL  UNSTEADY  MOTION 

o SUBSONIC  AND  SUPERSONIC 

° COMPUTATIONAL  EFFICIENCY 

CURRENT  DEVELOPMENTS: 

©NONLINEAR  EFFECTS  (TRANSONIC  FLOW. 

WAKE  DEFORMATION) 

t 

o IMPROVED  FINITE  ELEMENTS  (HIGHER  ORDER 
SPECIAL  PURPOSE) 

Q ROTATIONAL  FLOW  (TURBULENCE.  VISCOSITY) 
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In  the  integral  equation  for  the  velocity  potential  (previously  shown)  the  surface  integration  extends 
over  the  surface  of  the  aircraft  plus  its  wake,  and  the  no-penetration  boundary  condition  DS/Dt  = 0 
applies  over  both.  Moreover,  the  pressure  must  be  continuous  across  the  wake  although  the  potential  is 
discontinuous.  The  forward  edge  of  the  wake  of  a lifting  surface  is,  of  course,  always  located  at  the 
lifting-surface  trailing  edge,  but  the  position  of  the  rest  of  the  wake  is  not  known  a priori,  is  variable, 
time-dependent,  and  must  be  determined  in  the  calculation.  This  variability  requires  relocation  and 
reorientation  of  the  wake  and  its  panels  during  the  calculation,  perhaps  many  times  if  the  calculation 
is  iterative. 

This  figure  also  shows  a shockwave  which  is  isolated  from  the  flow  field  by  a portion  of  the  surface  S. 
Over  this  portion  of  S the  no-penetration  boundary  condition  must  be  replaced  by  Rankine-Hugoniot  conditions 
which  quantify  shock-induced  discontinuities  in  derivatives  of  the  potential  although  the  potential  itself 
is  continuous  across  the  shock.  These  discontinuities  make  it  desirable  to  have  panel  edges  lie  along 
the  foot  of  the  shock.  But  shock  location,  shape,  extent,  strength,  and  velocity  relative  to  the  vehicle 
surface  are  time  dependent.  Moreover,  motion  of  finite  amplitude — even  small  amplitude — that  is  needed 
to  investigate  limit-cycle  aeroelastic  response  can  lead  to  large-amplitude  shock  motion  and  even 
discontinuous  shock  location.  Consequently,  requiring  panel  edges  to  coincide  with  the  foot  of  the  shock 
can  require  extensive  repaneling  in  the  vicinity 'of  the  shock  during  calculations  for  unsteady  motion.  In 
contrast,  nonlinear  calculations  for  shock-free  transonic  flow  require  no  repaneling  and  impose  no  special 
requirements  for  surface  geometry. 


Shocks  and  wakes  from  lifting  surfaces  Impinge  upon  fuselages  or  other  portions  of  the 
vehicle  along  lines  that  are  «time  dependent.  Because  of  the  discontinuities  in  potential  or 
its  derivatives  at  these  impingement  lines.  It  is  desirable  that  panel  edges  coincide  with 
them.  Hence,  time  dependent  repaneling  in  these  vicinities  is  also  indicated. 

For  simplicity  in  its  development,  the  present  SfHJSSA  program  contains  zeroth-order 
(constant-potential)  aerodynamic  elements.  However,  it  has  been  intended  from  the  beginning 
that  the  program  would  employ  higher-order  elements  in  order  to  reduce  the  number  of  elements 
required  to  converge  the  solution.  Such  elements  have  been  developed  and  will  soon  be 
incorporated  into  the  program.  In  addition,  special-purpose  elements  are  being  developed 
for  paneling  in  regions  where  correct  variation  of  potential  is  theoretically  known.  These 
elements  have  built-in  shape  functions  to  produce  the  correct  variation  of  potential,  for 
example,  adjacent  to  normnlwash  discontinuities  such  as  control-surface  hinge  lines  and  side 
edges,  or  correct  variation  of  potential  derivatives  as  at  subsonic  trailing  edges.  Tn  addition, 
flow-through  elements  arc  required  to  model  engine  thrust  in  nacelles  and  lo  panel  shockwaves. 
Such  elements  impose  no  new  requirements  for  surface  geometry  information.  Hole  that  no  special 
panels  are  required  adjacent  to  shock  or  wake  Impingement  on  the  body  surface.  Use  of  higher- 
order  and  special-purpose  elements  should  reduce  computer  time  and  sLorage  requirements  but 
probably  will  do  little  to  reduce  the  amount  of  geometrical  input  information  required. 

Although  fewer  elements  are  used,  more  information  is  required  per  element.  Detailed  accuracy 
of  information  out  requires  detailed  accuracy  of  information  in  Irrespective  of  the  level  of 
sophistication. 

Finally,  incorporating  t^ie  effects  of  viscosity  and  rotational  flow  will  Impose  a 
requirement  for  relatively  high  accuracy  of  computed  pressure  gradients  and  hence  will  require 
higher-order  elements  (at  least  third  order)  than  would  he  required  for  most  potential-flow 
problems.  Alternatively,  it  is  possible  that  required  accuracy  and  order  of  continuity  may  be 
attained  from  solutions  using  lower-order  elements  followed  by  spline  (or  other)  Interpolation  of 
the  calculated  potential.  • 

Specific  treatment  of  aircraft  geometry  in  the  SOUSSA  program  itself  will  next  be  described  by 
Dr.  Morino. 
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The  rest  of  this  presentation  is  devoted  to  the  specific  geometry 
requirements  for  the  program  SOUSSA  P (Steady,  Oscillatory  and  Unsteady, 
Subsonic  and  Supersonic  Aerodynamics;  Production  Version).  As  presented 
abov^  the  Green’s  function  method  yields  an  integral  equation  over  the 
surface  of  the  aircraft  and  its  wake  (with  differential -del ay  dependence 
on  time).  Dividing  the  surfaces  in  quadrilateral  elements  and  assuming 
the  potential,  the  normalwash  and  the  potential  discontinuity  to  be  constant 
within  each  tlement  one  obtains  Eq.  (1). 

The  coefficients  B.,,  C..,  etc.,  are  evaluated  analytically,  with 
the  original  surface  approximated  by  a hyperboloidal  paraboloid 
(hyperboloidal  element).  Numerical  quadrature  is  used  for  distant 
element . 


In  order  to  complete  the  formulation  three  additional  relationships 
are  required:  * 

1.  Boundary  conditions,  relating  normalwash  [p  to  the  generalized 
coordinates  q^  (Eq.  2) 

2.  Bernoulli's  theorem  relating  pressure  coefficient  C to 

potential  $ (Eq.  3)  ^ 

3.  Definition  of  generalized  forces,  e , as  functionals  of  the 

pressure  coefficient  C (Eq.  4) 

F 

Finally  combining  Eqs.  1 to  1*  one  obtains  the  matrix  E relating  the 
generalized  forces,  e^,  to  the  generalized  coordinates,  q^. 
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This  elide  lows  the  flow  chart  for  the  program  SOUSSA,  and  It  la 
presented  In  order  to  Indicate  how  the  geometric  Information  la  used  In 
the  program.  The  check  points  will  be  discussed  later.  Here  only  the 
function  of  each  module  Is  briefly  described. 

Interfaces 


BODYG.  CONTG  amd  WAKEG:  Elaborate  the  geometry  Input  of  checkpoints  1 
and  2 (user  oriented)  Into  the  checkpoints  5,  6 and  7 as  needed  In  the 
rest  of  the  program. 

Potential-normal  wash  relationehlp  (mode  Independent ) 


COEFB: 

evaluates 

the  body  coefficients 

Bjh*  Cjh*  Djh  and 

V 

COEFW: 

evaluates 

the  wake  coefficients 

Fjn*  Gjn*  Snh*  6jn 

and 

YZMOD:  combines  the  above  frequency-independent  coefficients  to  yield 

the  frequency-dependent  matrices  [Y^]  and  UjjJ* 

Boundary  conditions  (mode  dependent) 

~ ^ 

E1MOD:  evaluates  the  matrix  E^  relating  jj;  to  ^ 

Bernoulli's  Theorem  (mode  independent) 

E3MOD:  evaluates  the  matrix  relating  Cp  to  ^ 

Generalized  Forces  (mode  dependent) 

EAMOD:  evaluates  the  matrix  E^  relating  e_  to 

Combining 

EMOD:  evaluates  the  generalized-aerodyr.amic-foree  matrix 

- = 

ADMOD:  implements  an  aerodynamic  design  method  which  yields  the  shape 

from  a prescribed  pressure  distribution. 
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This  slide  presents  the  contents  of  Checkpoint  ill  (input  to  module 
BODYG),  which  consists  of  Information  describing  the  geometry  of  the 
aircraft  body.  It  is  user-oriented  in  that  the  quantities  required  are 
compatible  with  the  output  of  state-of-the-art  geometry  preprocessors. 
Also,  if  the  aircraft  is  symmetric  with  respect  to  the  x-z  plane,  then 
only  the  right  half  need  be  supplied.  The  same  is  true  for  the  x-y  plane. 

Regarding  the  individual  components  of  Checkpoint  ill : 

o The  Cartesian  coordinates  of  the  nodes  are  assumed  to  be  already 
rotated;  that  is,  the  aircraft  is  oriented  as  desired  by  the 
user. 

o Referring  to  the  example  depicted  on  the  slide,  element  number  1, 
corner  1 yields  node  numbet  2. 

o The  body-symmetry  code  numbers  reflect  whether  symmetry  is 
considered  with  respect  to  the  x-z  and/or  v-z  planes. 

o The  element  code  numbers  provide  information  "uch  as  whether 

or  not  a wake  emanates  from  an  edge  of  an  element,  or  if  an  edge 
coincides  with  a hinge  line,  etc. 
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This  slide  presents  the  contents  of  Checkpoints  0 2 (input  to 
Module  WAKEG) , 03  (input  to  Module  E4M0D)  and  04  (input  to  Module  E1M0D) . 

Checkpoint  02  consists  of  information  describing  the  geometry  of  the 
wake.  By  describing  the  wake  as  a collection  of  strips,  many  different 
forms  of  input  can  easily  be  made  compatible.  If  a wake  strip  is  symmetric 
with  respect  to  the  x-z  plane,  then  only  the  right  half  need  be  supplied 
(same  for  the  x-y  plane) . 

Also, 

o Desired  orientation  of  the  wake  with  respect  to  the 
aircraft  is  assumed  to  already  be  satisfied. 

o The  matrix  that  relates  each  wake  strip  no.  with  the 

corresponding  four  trailing-edge  element  numbers  is  used  in 
evaluating  the  trailing  edge  values  of  the  potential  and 
for  determining  the  values  of  the  pressure  discontinuity 
at  the  centroid  of  the  trailing-edge  elements. 

Checkpoint  03  corresponds  to  the  generalized-forces  deformation 
modes,  and  Checkpoint  04  corresponds  to  the  boundary-condition  deformation 
modes . 
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CHECKPOINTS  ■*  3 AND  4 


This  slide  presents  the  contents  of  Checkpoints  05  (input  to  module 
COEFB)  and  06  (input  to  modules  COEFB  and  COEFW) . These  checkpoints  are 
at  a lower  level  than  checkpoints  01-4  if  one  views  the  SOUSSA  P flow 
diagram  as  a top-down  representation.  The  implications  of  this  are  that 
these  checkpoints  are  not  as  "user-oriented"  as  higher-level  checkpoints, 
since  program  execution  has  progressed  to  this  point.  This  is  evidenced 
by  the  fact  that  for  Checkpoint  05,  the  same  quantities  as  Checkpoint  01 
are  required  except  that  symmetry  conditions  (and  their  advantages  in 
preparing  geometrical  input)  are  not  considered.  Furthermore,  geometrical 
quantities  such  as  the  base  vectors  and  normals  of  rface  elements  are 
not  as  readily  available  from  geometry  preprocessc  as  the  information 
contained  in  Checkpoint  01.  These  considerations  must  be  accounted  for 
by  those  users  desiring  to  begin  execution  of  SOUSSA  P at  this  level. 

For  version  1.1  of  SOUSSA  P,  the  locatio'n  of  the  control  points  must 
be  specified  as  the  geometrical  centroids  of  the  body  elements.  For 
future  versions  (first-order  finite  element  formulation),  the  location 
of  the  nodes  will  be  the  necessary  input. 
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This  slide  presents  the  contents  of  Checkpoint  07  (Input  to  module 
COEFW) . This  Is  also  not  s top-level  checkpoint,  hence,  its  contents  may 
not  be  "user-oriented"  as,  say,  Checkpoint  02.  For  Instance,  at  this 
level : 

o The  coordinates  of  the  wake  elements  (as  opposed  to  the  wake 
strips)  are  required,  and  no  symmetry  conditions  may  be  taken 
advantage  of  (l.e.,  all  the  elements  must  be  Input). 

o The  matrix  used  in  correction  for  the  trailing-edge  potential 
values,  for  example,  must  be  given  for  the  elements  comprising 
the  complete  wake. 

o Most  geometry  preprocessois  would  not  provide  the  matrix  of  the 
coefficients  of  influence  of  the  trailing-edge  elements  that 
determiner  the  value  of  the  potential  discontinuity  for  each  wake 
element.  Note  for  SOUSSA  P 1.1  these  coefficients  are  simply  1 
and  -1,  but  for  later  versions,  splines  will  be  used  to  determine 
these  coefficients. 
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FLOW  FIELD  GRID  GENERATION 


1.  DISCUSSION  OF  CURVILINEAR  GRID  GENERATION 

2.  EFFECTS  OF  INACCURATE  SURFACE  REPRESENTATION 


The  accurate  three-dimensional  simulation  of  flow 
fields  by  finite  difference  or  related  methods  will 
require  a very  accurate  representation  of  the  surface 
geometry.  Continuity  of  surface  slope  and  curvature 
is  needed  unless  the  configuration  is,  indeed,  discon- 
tinuous in  these  features.  To  support  this  view  the 
problem  of  curvilinear  grid  generation  for  finite  dif- 
ference (or  finite  volume  and  some  finite  element) 
procedures  is  briefly  sketched.  The  sensitivity  of 
the  numerical  solution  to  inaccurate  surface  represen- 
tation is  then  illustrated  with  examples  (numerous 
examples  of  inaccuracy  have  been  generated  over  the 
years,  but  they  tend  to  be  quickly  discarded). 


Steger  1. 


Finite  difference  methods  are  used  when  nonlinear 
effects  such  as  compressibility  and  strong  viscous 
interaction  are  important.  In  a finite  difference 
problem  a solution  is  obtained  over  the  entire  flow 
field  domain. 


Steger  2. 


FINITE  DIFFERENCE  SIMULATION 


TRANSITION 


In  the  finite  difference  method  the  flow  field  is 
discretized  (or  meshed)  and  derivatives  are  replaced 
by  difference  approximations.  This  process  results  in 
a large  nonlinear  system  of  algebraic  equations  which 
may  require  simultaneous  solution.  Generally  the  equa- 
tions are  sparse  and  well  ordered  so  that  efficient  solu 
tion  methods  can  be  devised  that  are  often  amenable  to 
vectorized  computer  processing.  This  is  especially  true 
if  the  aerodynamic  configuration  is  forced  to  coincide 
with  a grid  surface. 

To  maintain  accuracy,  grid  points  should  be  clustered 
to  the  action  regions  of  the  flow  field.  In  this  exampl 
points  are  clustered  to  the  leading  and  trailing  edge 
of  the  airfoil.  Viscous  layers  are  resolved  by  cluster- 
ing to  the  airfoil  surface. 


Steger  3. 


In  one  form  of  grid  generation,  which  is  illustrated 
here  in  two  dimensions,  points  are  distributed  on  the 
body  and  outer  boundaries.  Curvilinear  coordinates  f, 
and  n are  then  generated  by  the  solution  of  an  elliptic 
equation  that  satisfies  a maximum  principle.  By  properly 
choosing  minimum  and  maximum  values  of  f,  and  n on  the 
boundaries  (see  sketch),  contour  levels  of  monotonicallv 
increasing  values  of  f,  and  n can  be  found  that  trace 
out  a curvilinear  coordinate  system. 

In  practice,  the  elliptic  equations  are  solved  in  a 
transformed  plane  along  specified  f,  and  n coordinate 
lines.  The  solution  for  x and  y along  f, ,n  coordi- 
nates in  the  transformed  plane  then  automatically  finds 
constant  lines  of  f,  and  n in  the  physical  plane. 

Accurate  surface  representation  enters  into  this 
process  only  once.  When  the  grid  points  are  specified 
along  the  body  surface,  they  must  lie  on  or  very  near 
the  correct  boundary  curve.  Otherwise  an  error  will 
result,  not  in  generating  a grid,  but  in  later  solving 
for  the  flow  field  about  the  correct  configuration. 


Steger  4. 


A grid  generated  by  the  previously  described 
procedure  prior  to  viscous  layer  reclustering. 


Steger  5. 


Grid  detail  near  body  surface  after  viscous 
layer  reclustering. 


Steger  6. 
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Small  inaccuracies  in  the  surface  representation 
of  a configuration  can  lead  to  much  larger  errors  in 
predicted  aerodynamic  quantities.  In  this  example 
(furnished  by  David  Nixon),  airfoil  ordinates  are 
slightly  altered  by  placing  sine-wave  bumps  on  the 
upper  surface  of  the  profile.  The  perturbed  ordinates 
are  always  within  1/2  percent  of  their  correct  local 
value,  yet  the  percentage  error  in  the  Cp  distribu- 
tion is  in  places  much  greater. 

In  general  the  computer  processing  work  of  gen- 
erating an  accurate  surface  representation  is  much 
less  than  the  work  in  obtaining  an  accurate  finite 
difference  simulation.  Consequently,  the  geometry 
should  be  much  more  accurately  represented  than  the 
estimated  accuracy  of  the  finite  difference  method. 


Steger  7. 
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In  this  example  of  viscous  transonic  flow  the 
radius  of  curvature  of  the  airfoil  actually  changes 
near  the  leading  edge.  The  finite  difference  scheme 
responds  to  the  change  with  the  peak  in  pressure 
distribution  shown  at  the  leading  edge.  A discon- 
tinuity in  curvature  due  to  inaccurate  surface  repre- 
sentation will  result  in  similar  peaks  in  the  pres- 
sure distribution. 


Steger  8. 


One  possible  mapping  scheme  for  three 
dimensions. 


Steger  9. 


Another  possible  mapping  scheme  for  three 
dimensions.  The  axis  singularity  is  definitely 
not  a problem  for  certain  formulations  of  the 
transformed  flow  equations. 


Steger  10. 
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WELL  ORDERED  GRID  MAPPINGS 


WARPED  SPHERICAL  MAPPING 
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FIGURE  2 - NASF  CONFIGURATION 


NASA  ,s  at  the  preliminary  definition  phase  of  a project  to  create  a 
numerical  calculation  of  turbulent  flow  over  complex  shapes.  A 
centra,  element  of  the  facility  is  the  high-speed  parallel  p-ocessor 
capable  of  computing  speeds  in  the  rang,  of  10  floating  point 
operat rons/sec.  (I  gigaflop).  An  essential  feature  of  the  facility 

'S  3 9raPMCS  Station'  Th‘  °»i«tive  this  presentation  is  to 
out  ;„e  the  genera,  reguirements  for  this  station  and  to  relate  the 
qualitative  nature  of  the  displays  desired  as  wel I as  the 
quantitative  leveis  of  data  required  to  create  such  displays 
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FIGURE  3 - SURFACE  GEOMETRY 


An  obvious  use  of  the  station  is  the  display  of  surface  geometry.  Several 

b a ,'-„Tdare-aVai'ab'e  5“Ch  diSD'aVS-  ^ conf  '9urat  ion  may  be  th0MI 
tha!  1,7  ^Stem  " 3 W.I  display  systems  exist 

l SUC  P'Ctures  to  be  rotated,  zoomed,  and  clipped  very  rapidly 
and  g,ve  the  impression  of  motion  on  the  screen.  This  pictare  pains  in 
rea 1 r sm  and  depth  by  dispiav  if  on,y  the  visibie  sarface  boundaries  A. 

-s  time  ,t  is  not  possible  to  nan i pa  late  sach  a display  in  real-tine 
becaase  of  the  heay y compatat iona . ,oad  repaired  to  sort  the  sarfaces 
into  hidden  and  visible.  A third  type  of  display,  the  shaded  sarface  is 
Possible  as ing  video  technigaes  rather  than  line  drawing.  Whi.e  these 
p.ctares  gam  a great  deal  of  realism,  they  lack  the  fine  isolation  of 

7 g SVSten'  ThiS  a“th°r's  •»'■"•»  - that  the  wire  frame 

type  display  is  the  mose  useful  for  this  station. 


FIGURE  4 - 2-D  AND  3-D  FLOW  FIELD  GRIDS 
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Previous  speakers  have  alluded  to  the  problem  of  computing  grid  points 
at  which  finite  difference  methods  are  used  to  compute  solutions  to 
non-linear  partial  differential  equations.  One  use  of  the  graphics 
station  is  the  display  of  such  grids.  While  this  is  clearly  feasible 
for  the  2-D  grid,  the  3-D  grid  contains  so  many  lines  that  it  is 
confusing.  Clearly,  some  innovative  techniques  will  be  required  to 
allow  the  user  to  understand  the  network. 


FIGURE  5 - PHYSICAL  RESULT 


The  results  of  the  aerodynamic  calculations  must  be  presented  to  the 
user  of  the  simulation  facility.  As  with  the  3-D  flow  field  qrids. 
considerable  innovation  will  be  required  to  present  meaningful  displays. 
A combination  of  dynamic  displays  with  variable  intensity  and  color 
will  probably  be  required. 
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FIGURE  6 - PRESSURE 


displays 


A number  or  well-known 
distributions  on  wings, 
blended  configurations. 


Th  are  aVanabU  f°r  of  pressure 

=se  are  not  easily  adapted  to  fuse, apes  and 

Aga.n,  considerate  innovation  „il,  be  repaired. 


FIGURE  7 - SURFACE  GEOMETRY  - WING 


The  next  7 figures  are  an  outline  of  the  data  requirement  for  surface 
geometry  and  flow  field  grids.  A wing  of  reasonable  complexity  could 
be  represented  by  either  of  1)  a dense  set  of  data;  2)  a set  of  spline 
curves;  or  3)  a set  of  parametric  patches.  Regardless  of  the 
representation  used,  the  data  requirements  are  approximately  103  points 


FIGURE  8 - SURFACE  GEOMETRY  - NACELLE 


a nacelle  requires  also  about  103  points. 
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FIGURE  9 - SURFACE  GEOMETRY  - FUSELAGE 


There  are  different  types  of  fuselages  - either  transport  or 
fighter.  It  appears  that  transport  fuselages  will  require  about 

1000  points  while  fighter  fuselages  will  require  about  5000 
points. 
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FIGURE  10  - SURFACE  GEOMETRY  - INTERSECTION 


In  the  region  of  component  intersections,  some  detailed  definition  is 
required.  A typical  intersection  could  have  approximately  3750  points. 


FIGURE  11  - SURFACE  GEOMETRY 


The  results  of  the  previous  U 

l *-•  - «» - a rat(Kr 

' °ff.2  P°dded  -"I  require  about  5000  poiou  I""""!* 

configuration  with  complex  fuselage,  several  i„t  ^ 

of  wings  and  tails  would  require  about  ,0=  points"^ 'CT  ^ ““  fU" 

airplane  configuration  with  control  surfaces  h h , ’ * ‘rU'y 

’nd  as  a wind-tunnel  i^del)  ciiuld  ea  ' ' SyS,em•  S*»™ 

Points  for  its  definition.  However  it  ■ ' . ^ °r  ev'n  '°7 

complexity  can  be  simulated  numerically  Conf^ ''Oorations  of  this 

me*  Ca  y the  foreseeable  future. 
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FIGURE  12  - FLOW  FIELD  GRID  GEOMETRY 


The  size  of  the  flow  field  grid  Is  strongly  related  to  th.  , 

different  types  of  grids  are  utilized  H „■  9rid‘  T"° 

, . * i zed,  depending  upon  whether 

be. ng  solved  are  inviscld  or  viscous.  The  ,*•  . equations 

grid  is  approximately  106  points.  SU  S'Ze  °f  ^ f,OW  f!e,d 


figure  13  - data  volumes 


From  the  previous  2 figures,  the  surface  geometry  and  grid  geometry  require- 
ments are  seen  to  be  about  104  and  106  , respectively.  AH  other  results 
are  related  to  these  figures.  About  10  physical  quantities  are  computed  at 
each  computation  point,  thereby  giving  107  points  in  the  field  and  105  on 
the  surface.  3y  integration  into  quantities  such  as  spare  loadings,  moments, 
etc.,  this  data  is  further  compressed. 
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FIGURE  1*»  - TIME  TO  MOVE  DATA 


Now  that  the  sizes  of  the  various  pieces  of  data  are  known,  the  times 
required  to  move  data  from  one  computer  to  another  can  be  estimated.  The 
speed  limits  associated  with  various  facilities  are  shown  on  the  horizontal 
axis.  Several  points  are  of  interest.  If  telephone  grade  lines  are  used, 
it  is  impractical  to  transmit  more  than  104  words.  Even  with  high  speed  lines, 
the  input/output  limitations  of  most  minicomputers  are  limtied  to  the  10**— 10s 
bit/second  range,  making  the  process  of  transmitting  a 3-D  flow  field  grid  a 
matter  of  an  hour  or  so. 
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MACHINE  SHOP  AND  WIND  TUNNEL 
MODEL  REQUIREMENTS 


WALTER  MANN  - NASA-AMES  RESEARCH  CENTER 
March  1,  1978 


Mann  1 


In  this  talk  It  Is  desired  to  focus  attention  on  and  place  In  perspective  the  machine  shop's  role 
In  wind  tunnel  research.  Clearly  this  role  Includes  the  fabrication  of  airfoil  shapes  that  are 

to  be  Instrumented  and  tested.  In  many  cases  It  also  Includes  fabrication  of  wind  tunnel  components 
such  as  rotor  blades. 


I 

i 


4 n 


Mann  3 


To  fabricate  the  airfoil  shapes  the  machine  shop  will  usually  employ  numerically  controlled  (N/C) 
machine  tools.  The  data  that  controls  the  operation  of  these  machine  tools  will  usually  be  gener 
ated  on  computers  by  N/C  programs  of  which  APT  Is  the  best  known  but  not  the  only  example.  It 
Is  at  this  level  that  I wish  to  focus  attention  - on  the  requirements  and  performances  oi  the  N/C 
programs  that  produce  the  cutter  path  data  for  control  of  the  fabrication  process. 
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instructions  for Tn  p r TTa 7 T “*■  **  Sh0P  **  -V  then  write 
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In  practice.  how»uo»  **®nn  7 

e*  nwever,  there  are  tun  ro,u»i 

The  flr  “',t,eS  Wm  *“*  **  -*  -t  contend. 

to  the  N/C  pro'graVanVmtu  We^  ftr^hT*5  fr°*’  the  engineer  may  not  be  In  a f 

<>e  translated,  rotated  or  scaled  i»*  '**  ™SSa9eii  or  "formulated.  if  the  °™  ““Stable 

generated,  say  for  a fillet  bet  ’ Pr<iSentS  Sreat  prob,em-  However  If  T °"'y 

titatlon  from  an  accor  c,  eT"  f"  .?  **  '*  * * 

scale  drawing  or  actua,  j.,  ‘ * extr«ted  by  d1„- 

The,  y we77  present  difficulty 

h second  is  that  the  cutter  Dath,i  i ** 

their  development,  are  not  yet  100  De9^  ^ °f  ***  N/C  programs » even  at  this  lat 

and  It  i*  tha  * , yet  100  Percent  reliable  th1s  late  stage  in 

form  appropr1auSl°fthhe  Mrt  ’‘r°Sr*mer-  -Here  he 'has  a chofcT  fT  ”°re  re',al,,e  others 

;:Ha? a,9°^-  - ^ zzzz the  data  *-  ti* 

three  years  It  win  un<ferway  to  resolve  this  1«t  “ Mtter  ,s  **•*  a 

"°t  * a consideration  In  part  fabrication.  an<l  h0pefu,1>  two  or 


TW  gEAUT/ES  MACHINE  SHOPMUST CONTEND 


Now  let  us  return  to  the  first  stated  reality  - the  need  to  reformulate  Input  geometry  to  accept- 
able form. 

Present  N/C  programs  taken  In  their  totality  allow  geometric  Input  of  three  basic  types: 

1.  Pure  analytical  - planes,  conics,  quadrics,  general  parametrlcs,  etc. 

2.  Pure  discrete  - points  and  vectors  for  tabulated  cylinders  and  meshes  of  points  and 
vectors  for  sculptured  surfaces. 


3.  Composite  of  both  analytical  and  discrete. 


N/c  PROGRAMS  ALLOW  THREE  R4SIC  TVP£S 

OF  GEOMETRY  PEHN  IT/ON 
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Mann  11 


If  the  engineer  describes  a shape  by  pure  analytic  means  (please  no  differential  equations)  then 
he  may  expect  an  adequate  representation  of  the  shape  In  the  fabricated  part.  However,  If  the 
engineer  describes  his  shape  with  discrete  data  (less  than  a semi-dense  set)  then  he  must  be 
aware  that  this  data  must  be  fitted  to  equations  by  the  part  programmer  before  it  Is  input  to  the 
N/C  program  or  that  Its  fitting  Is  a function  of  the  N/C  program. 


DISCRETE  DATA  MUST  BE  FITTED 
BY  EQUATIONS' 
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If  the  data  adequately  characterizes  the  shape  In  terms  of  the  fitting  procedure,  no  problem  arises 
is  at  rr  ^ q“eSt,0n  '*  a,W>yS  Whet"er  *he  d*ta  U ,Mde<|u,te  or  the  ^ttlng  procedure 


Thus  whenever  discrete  data  Is  presented  to  the  machine  shop  the  engineer  Is  not  out  of  the  woods 
until  he  can  accept  the  final  fabricated  shape  or  some  Intermediate  verification  data. 


WITH  DISCRETE  DATA  DEFINITION 
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What  then  Is  the  best  shape  definition  for  the  rnchlne  shop?  Analytic  data. 

best?  Discrete  or  composite  data  for  which  N/C  program  fitting  procedures  a 
able.  What  Is  least  desirable?  Just  points.  Procec-jre, 


best  geometry  definitions 


WORST 
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IPAD  - GENERAL  SYSTEM  ARCHITECTURE 


IPAD  as  a system  has  been  designed  to  this  point  within  a total  corporate  complex  and  Its  use  of  computers. 
Part  A of  this  vlewfoll  shows  the  topmost  consideration  in  this  general  system  architecture  and  envisions 
that  a corporate  complex  is  comprised  of  several  CAO/'.AM  complexes  which  communicate  with  each  other 
through  a corporate  network. 

The  CAD/CAM  complex  in  turn  is  composed  of  several  computer  complexes,  each  of  which  have  their  local 
data  base  and  they  in  turn  communicate  with  each  other  through  a local  network  which  has  attached  to  the 
local  network,  a global  data  base  and  an  IPAD  global  data  base. 

In  general,  these  two  networks  have  markedly  different  comnuni cat ions  characteristics.  The  corporate 
network  is  characterized  by  the  utilization  of  normal  communications  or  microwave  media  and  is  character- 
ized by  spted  in  the  order  of  kilo  baud.  On  the  other  hand,  the  local  network  ucllitizes  specialized 
communications  technology  and  a series  of  microprocessors.  These  microprocessors  provide  the  means  for 
attaching  heterogeneous  computers  to  the  local  network.  The  local  network  operates  in  distances  in  the 
order  of  thousands  of  feet  and  at  speeds  in  the  order  of  megabaud. 

Within  any  one  of  these  computer  complexes,  we  find  a host  computer,  Its  local  data  base  and  an  IPAD  local 
data  base.  An  Interface  between  the  host  and  the  IPAD  system,  the  IPAD  system  and  a body  of  non-IPAD 
program  application  programs  as  well. 


INTEGRATION  OF  GEOMETRY-GRAPHICS  INTO  IPAD  SYSTEM 


This  slide  repeats  some  of  the  elements  of  a computer  complex  and  some  of  the  events  of  the  IPAD  system 
architecture.  In  addition,  as  shown  in  the  highlighted  portion  of  the  slide,  the  particular  components 
that  are  associated  with  the  Geometry-graphics  system  are  indicated.  Shown  here  are:  1)  IPAD  Geometry 
Standard  Utility,  High-Level  Graphics  (HLGR)  which  is  optional.  HLGR  provides  high-level  FORTRAN-callable 
subroutines  that  is  not  purely  graphics  (e.g.,  data  modeling,  plotting  formats).  They  can  be  considered 
graphics-related  macros  which  are  used  to  generate  frequently  used  displays;  and  2)  the  IPAD  graphics 
primitives. 

The  IPAD  graphics  primitives  and  the  IPAD  geometry  standard  utility  will  be  discussed  in  detail  in  the 
following  slides.  It  should  be  noted  from  this  slide  that  the  graphics  user  interfaces  with  the  host 
operating  system  (using  graphics  interface  software)  who  in  turn  interfaces  with  the  IPAD  executive. 
Through  the  IPAD  executive,  the  graphics  user  has  access  to  the  analysis  programs  and  the  various  compo- 
nents of  the  IPAD  system. as  shown.  The  graphics  user  also  has  access  to  the  data  base  through  the  IPAD 
information  processor,  and  access  to  the  local  network  and  all  the  resources  that  cascade  outward  from 
the  local  network  as  shown  on  the  earlier  architecture  slides. 


REM 

| 3-1-73 

PAGE  5 

1 — 
i — 1 

U- 

O 

INTEGRATION  OF  GEOMETRY  - GRAPHICS  INTO  IPAD  SYSTEM 


DISTRIBUTED  IPAD  LOCAL  NETWORK  EXAMPLE 


As  mentioned,  the  graphics  user  has  access  through  the  host  operating  system  to  the  local  network.  The 
IPAD  architecture  anticipates  as  typical  Implementations  of  this  distributed  architecture,  a variety  of 
processers  (computers)  which  might  be  specialized  to  particular  tasks  or  processing  requirement  ( i . e . , 
graphics  nodes).  Such  an  example  of  processing  requirements  is  shown  on  this  slide  where,  associated 
with  the  prime  might  be  finite  element  modeling,  associated  with  a larger  processor  like  a CDC  CYBER 
finite  element  analysis  and  so  forth,  as  illustrated.  It  is  the  intent,  presently,  of  the  IPAD  develop- 
ment, to  allocate  a CAD  system  such  as  AD2000  to  a DEC  PDP/11-70  processor  and  to  associate  this  processor 
with  other  processors  in  the  network  as  shown  here. 
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IPAD  GEOMETRY  STANDARDS 


This  slide  discusses  various  standards  items  and  their  impact  on  the  IPAD  system  design.  Standards  would 
be  concerned  with  primitives  covering  point,  curve  and  so  forth  as  well  as  primitives  for  text  and  dlrnen- 
tioning,  etc.  Such  primitives  thus  cover  the  aspects  of  geometry  as  well  as  drafting.  Secondly,  rela- 
tionships between  the  primitives  are  of  concern  covering  items  such  as  union  replication.  Thirdly,  geome- 
try model  management  is  also  a concern.  Model  management  cover  aspects  such  as  save,  restore,  and  copy 
for  dealing  with  models  which  have  been  constructed  through  the  use  of  primitives  and  primitive  relation- 
ships. Such  models  represent  the  needs  of  analysis,  parts  for  manufacturing,  kinematic  studies,  and  so 
forth. 

These  standards  impact  various  aspects  of  the  IPAD  system  design.  As  an  example,  it  is  anticipated  in  the 
language  area  of  the  IPAD  system  that  a geometry  data  definition  language  and  a geometry  data  management 
language  will  be  required  to  be  implemented,  and  in  the  case  of  the  geometry  definition,  all  of  the  primi- 
tives are  o*  direct  impact  upon  this  language.  Similarly,  the  primitive  relationships  have  a direct  im- 
pact on  the  manipulation  language. 

The  data  management  aspects  of  the  IPAD  system  are  impacted  very  heavily  by  the  requirements  of  model 
management.  Thus  IPIP,  the  data  base  management  system  utility  in  IPAD,  will  be  directly  influenced  by 
the  data  structures  that  are  implied  in  the  model  management  aspects  of  the  standard  as  well  as  the  parti- 
cular data  structures  associated  with  the  primitives  and  the  primitive  relationships.  This  portion  of  the 
IPAD  data  manager  will  have  to  be  tailored  for  high  performance  to  meet  the  response  time  requirements  of 
a highly  interactive  user  engaged  in  geometry  modeling. 

The  data  communication  implications  of  the  standards  will  have  an  effect  upon  both  the  local  data  base  and 
global  data  base  aspects  of  the  architecture  as  shown  on  the  earlier  configuration  charts.  The  highly 
interactive  dialog  work  is  expected  to  be  handled,  on  a data  communication  basis,  through  the  local  net- 
work and  to  be  handled  between  not  only  the  local  data  base  to  the  local  data  base,  but  also  be  handled 
Internally  within  any  given  processor  to  the  local  data  base  directly.  Large  data  volumes  and  traffic 
over  longer  distances  which  might  take  place  between  the  various  CAD  complexes  as  shown  on  the  earlier 
architectural  slides,  would  be  handled  at  lower  speeds  on  a corporate  network. 
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STANDARD  ITEM 


IPAD  SYSTEM  DESIGN  IMPACT 


PRIMITIVES.  — 

POINT 
CURVE 
SURFACE 
SOLID 
TEXT 

DIMENSIONS 


GEOMETRY 


DRAFTING 


♦ 


PRIMITIVE  RELATIONSHIPS 
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UNION 

REPLICATION 

DELETE 


LANGUAGES, 

GDDL  - GEOMETRY  DATA  DEFINITION  LANG. 

GDML  - GEOMETRY  DATA  MANIPULATION  LANG. 
' QUERY  - PRIMITIVE  OR  MODEL  LEVEL 
MODEL  MGT. 

DATA  MANAGEMENT 
IPIP  (DBMS)  UTILITY 

- DATA  STRUCTURES  SUPPORTED 

- TAILORED  FOR  HIGH  PERFORMANCE 


SAVE 

RESTORE 

COPY 


DAJAXQMUNICAILQN 

- DATA  VOLUMES  AND  TRAFFIC  ON 
LOCAL  NETWORK 

- RESPONSE  FOR  HIGHLY  INTERACTIVE 
DIALOGUE 
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GEOMETRY  MODELING  PROCESS 


The  various  elements  of  the  IPAD  system  and  the  geometry  standard  utilities  are  shown  in  relation  to  the 
geometry  modeling  process. 

A user,  through  his  terminal  and  the  host  system,  is  attached  to  the  IPAD  executive  IPEX.  Through  IPEX 
the  user  has  1)  access  to  IPIP  for  direct  access  to  GDDL,  GDML,  Query  and  Model  Management  activities  in 
the  Local  IPAD  D.B.  or  the  Global  IPAD  D.  B.,  and  2)  access  to  the  IPAD  CAD  system  and  the  IPAD  geometry 
standard  utility  for  constructing,  manipulating,  editing  and  managing  3D  geometry  models.  During  the 
modeling  process,  the  model  will  be  in  the  Local  IPAD  D.B.  However,  the  model  management  utilities  of 
IPIP  provide  for  functions  between  the  Local  and  Global  Data  Bases.  Finally,  the  various  geometry  models 
of  interest  to  a community  of  users  are  stored  in  the  Global  IPAD  D.B. 
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USER  TERMINAL 


GEOMETRY  MODELING  PROCESS 


I P E X 


CAD  SYSTEM 


-SUBROUTINE  CALLS 


I P I P 

USER  INTERFACE 


DIRECT  USER  INPUTS 


1 

G D D L 

G D M L 

QUERY 

P T 1 = POINT  (x,Y,Z) 

SET  1 = UNION  (LN  1,  LN  2,  LN  3) 

PT  3 = ASSOC  (point) 

L N Z = LINE  (PTl>  P T z) 

SET  5 = UNION  (SET  L SET  2) 

ARC  3 = CANON  (A  R C) 

• 

• 

DELETE  = (POINT,  PT  1) 

A 

• 

• 

• 

W 

• 

• 

3-d  geometry  model 


LOCAL  DATA  BASE 


MODEL  MANAGEMENT 
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■GLOBAL  DATA  BASE 
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IPAD  GRAPHIC  STANDARD 


This  chart  illustrates  the  various  elements  that  are  utilized  when  displaying  pictures  of  models  which  have 
been  created,  and  illustrates  in  addition,  the  location  of  the  IPAD  graphic  standard  elements.  The  user 
would  utilize,  from  the  IPAD  data  base,  various  analyses,  graphics  utilities,  and  high-level  graphics  rou- 
tines to  construct,  through  the  use  of  the  graphics  primitives,  a set  of  pictures,  finally  building  up  a 
virtual  display  file.  This  display  file  could  be  saved  for  other  postprocessing  needs.  The  display  file 
is  passed  to  the  device  interface  software.  The  display  file  could  be  passed  directly  to  the  smart  termi- 
nal. In  the  case  of  dumb  terminals,  the  device  interface  software  would  map  the  display  file  into  the 
routines  associated  with  displaying  the  picture.  Similar  processes  would  take  place  to  off-line  plotting 
devices. 

The  core  primitives  are  the  simplest  graphics  tools.  Primitives  include:  Input,  output,  segments,  attri- 

butes, view  transformations,  and  control.  They  embody  the  concept  of  portability  and  graphics  standardi- 
zation. The  graphics  primitives  are  designed  to  allow  easy  creation  and  component  modification  of  graphic 
displays.  Primitives  are  not  designed  to  modify  the  user's  data  structures  directly  when  CRT  changes  are 
made.  The  responsibility  for  modifying  data  structure  resides  with  the  application  program  responsible 
for  the  data  structure  (e.g.,  CAD).  However,  "hooks"  are  provided  by  the  primitives  to  allow  a program  to 
propagate  changes  made  on  the  CRT  in  the  data  structure  produced  by  the  application  program. 
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GEOMETRY  AND  GRAPHICS  STANDARDS  IN  CAD  SYSTEMS 


Generically,  a CAD  system  has  such  major  elements  as  geometry-related  analysis  routines,  transformation 
capability  between  the  geometry  and  the  graphics  routines,  provision  for  user  interaction  data  management, 
and  contains  procedures  for  bridling  both  geometry  and  graphics.  Typical  of  such  systems  are  CADD,  AD2000, 
and  CADAM.  Each  CAD  system  has  its  own  unique  method  of  handling  geometry.  Therefore,  corrmuni cation  be- 
tween CAD/CAM  systems  and  related  D.B.'s  become  a major  problem.  The  IPAD  geometry  standard  will  provide 
a common  basis  to  which  each  can  be  interfaced  and  thus  reduces  the  number  of  translators  required.  Such 
systems  are  characteristically  used  to  accept  data  from  a user  and  to  produce  either  models,  geometric 
models  of  particular  entities  of  interest  and  also  to  display  pictures  of  those  models.  Modeling  process 
uses  geometry  standards.  In  this  picture  we  show  that  the  IPAD  geometry  standard  is  intended  to  encom- 
pass the  work  to  date  of  the  ANSI  Y14.26.1  committee,  plus  additional  items  such  as  drafting  standards, 
which  were  mentioned  earlier.  Similarly,  there  is  an  IPAD  graphic  standard  which  is  based  primarily  on 
the  work  of  an  ACM  SIGGRAPH  and  augmented  by  basic  drawing  elements.  It  is  the  intent  of  the  IPAD  program 
to  develop  an  IPAD  geometry  standard  based  upon  the  ANSI  work  and  an  IPAD  graphics  standard,  based  upon 
ACM  SIGGRAPH  work. 
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FUTURE  USE  OF  IPAD  GEOMETRY  AND  GRAPHICS  STANDARDS 


IPAD  graphics  standard  has  the  potential  of  being  adopted  by  such  standards  bodies  as  the  ACM  SIGGRAPH 
and  the  ANSI  committees.  In  the  future  we  would  expect  that  systems  would  be  designed  and  built  with 
this  standard  as  a basis  and  that  for  the  Interim,  the  existence  of  the  IPAD  standards  would  expedite 
♦he  integration  of  such  systems  Into  an  IPAD  environment. 

From  the  point  of  existing  systems,  two  different  strategies  are  available  to  us.  With  respect  to  geo- 
metry, we  have  the  option  of  translating,  as  shown  on  the  previous  slide,  between  the  IPAD  geometry  stan- 
dard and  the  geometry  standard  of  the  particular  commercial  CAD  system  in  use,  thus  creating  an  interface. 
Or,  directly  replacing  those  particular  geometry  portions  of  the  conmercial  CAD  system  and  associated  data 
management.  In  the  case  of  graphics,  we  anticipate  that  it  will  be  necessary  to  replace  the  graphics  rou- 
tines directly  with  the  ' 'AD  graphics  primitives.  The  geometry  and  graphics  associated  routines  In 
AD2000  are  being  investigated  at  our  development  site,  to  determine  the  relative  effort  to  interface  or 
replace. 


FUTURE  USE  OF  IPAD  GEOMETRY  AND  GRAPHICS  STANDARDS 
FUTURE  SYSTEMS 

o DESIGN  AND  BUILD  WITH  STANDARDS  AS  A BASIS, 
o EXPEDITE  INTEGRATION  INTO  AN  IPAD  ENVIRONMENT, 


existing  systems 

GEOMETRY  ASPECTS 
U)  INTERFACE  TO 

OR 

(2)  REPLACE  WITH 

o IPAD  GEOMETRY  DATA 

o MANAGEMENT  UTILITY  AND 

ASSOCIATED  DATA  BASE 


GRAPHICS  ASPECTS 

o REPLACE  WITH  IPAD 
GRAPHICS  PRIMITIVES 


QUESTIONS  THAT  MIGHT  BE  DISCUSSED  IN  PANELS 


PANEL  I:  REQUIREMENTS 

1.  How  accurate  must  be  the  surface  location,  slope  and  curvatures  to  give 
acceptable  aerodynamic  calculation  results? 

2.  Spline  methods  (linear  or  higher  order  interpolation)  introduce  some  degree 
of  artificial  waviness.  At  what  point  does  this  introduce  significant 
errors  in  aerodynamic  calculation  (or  measurement) 7 

3.  Real  aerodynamic  shapes  usually  have  discontinuities  in  the  surface  such 
as  sharp  edges,  corners,  gaps  and  inlet/exhaust  holes.  How  precisely  do 
these  need  to  be  specified? 

*».  Aerodynamic  calculations  often  involve  several  iterations  or  cases  that 
call  for  minor  adjustments  in  surface  shape,  flow  conditions  or  paneling 
density/distribution.  Whet  is  a “reasonable  turnaround  time"  between 
aerodynamic  analyses?  This  question  could  be  stated  as:  What  "turnaround 

time"  would  cause  ae.-odynjmic  analysts  or  researchers  to  abandon  or  hardly 
use  an  available  computerized  geometry  method?  i'ow  much  input  and  re- 
input would  be  inlii  biting? 

5-  The  original  aerodynamic  shape  specification  can  come  from  a number  of 
sources;  drawings  at  several  levels  of  detail,  lists  of  coordinates, 
systems  of  equations  or  coefficients,  actual  wind  tunnel  models  or  aircraft. 
Should  a surface  representation  system  accomodate  all  these  sources?  What 
checking  methods  are  needed? 

6.  Is  it  time  to  discuss  standardization? 
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Analytical  Surface  Description 


Parametric  representation  of  unit  spherical  surface 


X ■ cos  9 sin  1 
Y ■ sin  9 sin  ® 
z * cos  f 


Two  parameters  8,®  required  to  generate  surface  in  3-D 


Vector  parametric  equation 


S ■ (sin®  cosfl)e,  ♦ (sin®  sln8)et  + (cns$)Cj 


Parametric  tangents 

?i  , *4 
* a® 

Tv/ 1 s t vectors#  (Cross  derivatives) 

iL  £L 

5 Pa®  * 9®d& 

Surface  normals  (Cross  rroduct  of  two  tanqent  vectors) 


N - 9 s x »S 


Analytical  Surface  Description 


ADVANTAGES 

Exact  mathematical  description 

Mo  storing  large  quantities  of  data  points 

Secondary  surface  properties  easily  derived 

t.e.  surface  areas,  volumes,  curvatures,  etc. 
No  Interpolation  schemes  necessary 


DISADVANTAGES 

Limited  class  of  surfaces  addressable 
Mathematics  of  composite  objects  relatively  complex 
Problems  not  yet  solved  Include  boundary  conditions, 
smootheness,  oblique  angle  intersections,  fillets 


DESIGN  SYSTEMS  / I'TILITY 

PADL  system  (Part  and  Assembly  Description  Language) 
Discrete  part  description 


PARE  l» 


Bilinear  Surface  Description 


MATHEMATICS 


Q(u,w)  «[(l-u)  ul  fp(0/0) 

|P(1'n) 


9(0,1)]  [i-wl 

9(1,1)]  [w  J 


ADVANTAGES 

Simple  to  construct/ Implement 

Points  on  surface  linear  combination  of  patch  endpoints 
I’ser  not  Involved  with  parametric  tangents,  etc. 


DISADVANTAGES 

Not  flexible  multi-purpose  technloue 
A single  natch  Is  not  a curved  surface 
Mo  rroro  than  first  order  continuity  anywhere 
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Lofted  or  Ruled  Surface  Description 
MATHEMATICS 

Q(u,w)  - P(u,0)  (1-w)  ♦ P(u,l)  w 


ADVANTAGES 
General  Scheme 

Can  accommodate  different  curve  types(pol ynomf al ,spl Ine,etc) 


DISADVANTAGES 

Contains  curves  on  only  two  sides  of  patch 
Linear  Interpolation  scheme  In  one  direction 
C fairness  across  patch  boundaries 


Linear  Coons'  Surface 


MATHEMATICS 


Q(u/w) 


ADVANTAGES 

Allows  curve  on  all  four  sides  of  natch 
Supports  different  curve  descriptions 
Linear  blending  easy  to  Implement 


DISADVANTAGES 

Some  curve  types  could  be  Involved  mathematical! 
Linear  blending  functions  determine  Internal 
surface  shape  (regardless  of  curve  type) 


Overhauser-foons*  Surface 


PARE  7 


mathematics 


Q(u,w) 


- p ( 0^'  0 ) ) ^(  ° 5 0 : w ’ u >*p  ( U a ) n,(  w ) ♦ 


P(1,w)B4(u) 


whe  re : 


B.(t)  • 1 - 3t*  ♦ 2t* 
*,<*>  - 3t*  - 2t 


,SMU,n*  functions 


"™f"J  ^e“^tcahr%‘^-^-ser  curves 


Overhauser-Coons ' Surrace 


ADVANTAGES 

Improved  algorithm  for  fast  computation 
Facilitates  the  use  of  'Shaping  Tools' 

Better  local  control  than  cubic  curve  techniques 


DISADVANTAGES 

Large  node  displacement  may  cause  spurious  wiggles 
Adding  a point  to  curve  affects  3 curve  segments 
Moving  a point  affects  4 neighboring  curve  segments 
Normally  12  adjacent  patches  are  affected  by  the 

displacement  of  a single  point  In  common  to  4 patches 
Complex  data  structure/data  handling  facilities 
Storage  requirements  - data  points  and  tangents 
c'  fairness  across  patch  boundaries 


DESIGN  SYSTEM  / UTILITY 

Three  Dimensional  Design  System  --  J.  Brewer  ( 1977-Purdue) 
Total  software  implementation  (no  special  hardware  required) 


Bicubic  Surface  Description 


MATHEMATICS 


Q(u,w) 


[ 


B,(u)  B/u)  tyu) 


P 


B»Cw> 
V w) 
Vw) 
Itfw) 


where  Bl(t)  arc  the  cubic  'Blendlnjr  Functions' 


[b,( t ) B/t)  B,(t)  tyt)]  - [tJ  t t 
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and  P Is  the  boundary  condition  matrix 
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w-tangent 
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Plcublc  Surface  Description 


ADVANTAGES 

Well  proven  and  widely  used  method 


DISADVANTAGES 

\ 

Not  readily  adaptable  to  surface  'Shaping  Tools' 

Three  ouantltles  (position  , tangent  vectors, twl st  vectors) 
all  must  be  worked  with  to  manipulate  surface  shape 
Hard  to  have  Intuitive  feel  for  surface  shape  control 
Calculation  difficulty  Is  high  (matrix  Inversion) 

Poor  local  curve  control 

Must  split  curve  to  represent  knuckles 

C*  fairness 

Large  storage  requirements 

Spurious  wiggles  from  parametric  cuhlcs 


DESIGN  SYSTEM  / UTILITY 
NASA  / Ames  system 

HIM  DEF  --  Naval  Ships  Engineering  Center 


Planned  Design 
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Berler  Curve  anH  Surface  Description 
MATHFMATICS 


Turves 

«(s)  -Si  Pt 

where  fj(s)  - (l)  s‘ (l-s)"'* 
5 - polyron  vertices 


Cartesian  Product  (Tensor  Product)  Surface 


Vs'*’  ■ ££><*>  f>(t>  ?, 
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Rezler  Curve  and  Surface  Description 
ADVANTAGES 

Low  storage  requirements  --  polygon  points 
High  order  curve  and  surface  fairness 
Intuitive  user  Interface  to  surface  control 
tangent  and  twist  vector  control  parameters 
manipulated  by  placement  of  polygon  points 
Interactive  curve  approximation  shown  effective 
Ah  Initio  surface  design  capabilities 
Variation  diminishing  properties 


DISADVANTAGES 
Local  curve  control  poor 

Concatenating  Rezler  curves  and  surfaces  not  well  developed 

Using  high  degree  polynomials,  some  smoothness  lost 
computation  costs  Increase  proport lonatel y 


DFSIGN  SYSTEM  / 1'TILITY 


UNISJIRF  System  --  Regie  Renault 


B-Spllne  Curve  end  Surface  Description 


MATHEMATICS 


Curves 


C(t)  - SljlLIt) 


where  Pj  * polygon  vertices 
M.  ■ basts  functions 


the  N matrix 


N 


Ni-«  j . . . 


# 


0 


Cartesian  Product  Surface 


Ctfu,w) 


■SSis 


% (u)  Mm(w) 
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B-Spllne  Curve  and  Surface  Description 
ADVANTAGES 

Well  conditioned  for  curve  order  less  than  20 
Local  basis  s.t.  at  every  point  only  k basis  A 0 
Low  storage  requirements  --  polygon  points 
R-Spllne  formulation  contains  Bezier  as  a special  case 
Good  user  control  handles  for  manipulating  surface 
No.  and  placement  of  polygon  vertices 
Multiple  (repeated)  polygon  vertices 
Order  of  curve  2^  k$  no.  of  vertices 
Ab  Initio  curve  and  surface  design 

Follows  polygonal  form  more  closely  than  Bezier  curve 
Good  local  curve  control 
Variation  diminishing  property 


DISADVANTAGES 

Mot  well  developed  as  a curve  fitting  technique 
Calculation  difficulty  moderate  to  high 


DESIGN  SYSTEM  / UTILITY 

University  of  Utah  System  — J.H.  Clark  (1974) 
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SUMMARY 


In  conjunction  with  research  being  conducted  by  the 
Mechanical  Engineering  Department  at  that  University  of 
Connecticut  In  the  area  of  three  dimensional  model 
generation,  a review  of  the  literature  has  been  conducted. 
The  topics  of  primary  Interest  were  (1)  Mathematical 
techniques  for  three  dimensional  surface  representat Ion  and 
design  and  (2)  State  of  the  art  In  three  dimensional 
computer  model  generation.  This  presentation  to  the  NASA  / 
Ames  workshop  on  Aircraft  Surface  Representation  for 
Aeordynamtc  Computation  Is  a summary  of  the  Investigation 
Into  topic  (1)  above. 
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Introduction 


Curve  and  surface  algorithms  have  received  ever-increasing 
attention  during  recent  years.  Bezier  and  Riesenfeld  introduced 
algorithms  that  allow  more  intuitive  control  of  shape  than  the 
earlier  developed  Coons  formulation.  A lesser  known  parametric 
curve  description,  known  as  parabolic  blending,  was  developed 
by  A.  W.  Overhauser  in  1968  at  Ford  Motor  Company.  This  formu- 
lation offers  unique  advantages  to  interactive  curve  and  surface 
manipulation  because  it  is  based  solely  on  coordinates  on  the 
curves,  and  not  parametric  derivatives. 

Formulation  of  Overhauser  Curves  from  Parabolic  Blending. 

c(+)»  0-+)f(r) 

Parabolic  blending  consists  of  a parametric  blend  of  two 
parabolas  p(r)  and  q(s). 
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p(r)  and  q(s)  are  plane  parabolas  constructed  from  the  first 
three  points,  and  the  second  three,  respectively. 


The  curve  is  a blend  of  these  parabolas  between  the  inner 
two  points. 
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Each  parabola  can  be  expressed  as  a quadratic  in  its 
parameter . 
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These  can  be  combined  via  a linear  relationship  between 
r and  t,  and  s and  t,  resulting  in  a cubic  space  curve. 


t 

/ 

/ Ml 


✓ 

» 

Ts 


q (•) 


FI**1*  1 A n Ovtrhauicr  Cur v« 


l*)  * I"*-’  -t*  -t  ^ 


% -v» 

, -*■/.  2 

p* 

O ■'*  ° 

?, 

Of  O O _ 

A 

Pq 

-a. 


By  arbitrarily  assigning  a value  of 
compact,  easy  to  construct  formulation  is 
on  the  space  coordinates  on  the  curve. 


r ■ 1/2  at  t • 0,  a 
found,  depending  only 


A series  of  blended  curves  is  constructed  easily  by 
sequencing  the  selection  of  the  four  points.  First  derivative 
continuity  is  maintained  as  a consequence  of  the  c(t)  formulation. 


If  a discontinuity  in  the  first  derivative  is  desired,  off- 
curve  points  can  be  introduced.  In  this  manner,  these  points 
allow  flexible  shape  control  similar  to  Bezier’s  design  points. 


t 

< 

qU) 

H fur*  4 An  Ov*rhowt*r  Cur**  *IU 
Iud  Point* 


This  formulation  also  allows  two  points  to  be  coincident 


without  numerical  difficulties. 


A surface  can  he  constructed  from  Overhauser  curves  by 
several  methods.  Overhauser  suggested  that  sets  of  four  points 
on  "parallel"  curves  be  used  to  form  a parabolic  "blend"  across 
an  interior  patch.  The  "traces"  could  be  calculated  in  the 
opposite  parametric  direction,  or  the  two  could  be  averaged  if 
desired.  In  any  case,  the  surface  shape  depends  upon  sixteen 
adjacent  coordinates. 

I nt  cract i on 

Overhauser  curves  and  surfaces  were  added  to  an  interact  iv 
30  design  system  being  developed  in  the  Computer- Oes ign  and 
Graphics  Laboratory  in  the  School  of  Mechanical  engineering  at 
Purdue  University.  Ihe  system  also  allowed  Coons'  bicubic 
surfaces  to  be  interactively  manipulated  using  an  Imlac  PUS-1 
refresh  display  connected  to  a POP  11/40  minicomputer. 


An  Overhauser  Curve  is  constructed  by  placing  points  in 
3 space  with  a dynamic  3D  cursor.  The  curve  is  m nipulated  by 
"picking- up"  a point  and  dragging  it.  Affected  segments  of 
the  curve  are  intensified  by  the  system.  New  points  can  be 
easily  added  to  the  curve  in  the  same  manner. 


rtAj-Dces«vs/  - S' 


A group  of  points  forming  a curve  or  surface  can  be 
manipulated  using  a shaping  tool,  such  as  a plane.  The  3D 
plane  is  defined  by  the  designer,  and  points  "stick"  to  it 
as  it  moves  in  space.  This  illustrates  one  flexibility 
afforded  by  the  Overhauser  formulation. 

The  ten  minute  film  shows  the  3P  design  system  in  action, 
including  the  manipulation  of  Coons'  and  Overhauser  curves 
and  surfaces. 
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An  Overhauser  curve  is  created  using  the  3D  cursor 
a planar  shape  tool  is  used  to  "flatten"  the  curve. 


(•)  (b) 


A cylinder  consist’ng  of  eight  Overhauser  surfaces  is 
generated  and  truncated  with  the  planar  shaping  tool. 


and 
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REPRESENTS  GRAPHICAL  DISPLAYS 


Graphical  displays  can  take  many  different  shapes  and  forms.  The  ultimate  goal 
Is  to  obtain  as  high  a quality  representation  of  the  actual  data  as  possible  as 
cheaply  and  gulckly  as  we  can.  In  doing  this,  we  should  consider  the  ultimate 
potential  of  the  plotting  device  as  well  as  the  data  Itself. 
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RASTER  space-line  plotting  concept 


All  raster-oriented  plotting  devices  use  some  sort  of  , 
method.  The  raster  cnaro  K curve  approximation 
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POINT-THINNING  CONCEPT 

A curvilinear  line  may  be  broken  down  Into  line  segments  at  the  programming  level 
rather  than  the  plotting  device  level.  This  Is  especially  beneficial  In  working  with  very 
dense  data  representations,  since  many  calls  to  the  actual  plotting  routines  may  then  be 
bipassed.  In  this  example,  a reasonable  curve  representation  could  be  specified  with  just 
nine  calls  (circled  points)  to  the  plotting  routines. 
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FORTRAN  IMPLEMENTATION  OF  POINT-THINNING  LOGIC 


X 


FORTRAN  code  Implementing  a point-thinning  technique  Is  Illustrated  by  the 
DO  loop  as  shown.  Input  x-y  data  to  be  plotted  Is  provided  In  the  XLOC  rnd  YLOC 
arrays.  The  AMULT  factor  Is  the  number  of  rasters  per  unit  of  plotting  length, 
o.g.,  200  rasters  per  Inch  for  most  CALCOMP  devices.  The  I FIX  function  converts 
a floating  point  number  to  the  Integer  mode.  The  CALPLT  routine  Is  the  actual 
plotting  routine.  In  Its  call,  the  third  parameter  specifies  pen  up  (a3)  or  pen 
down  (*2). 
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FORTRAN  IMPLEMENTATION  OF  POINT-THINNING  LOGIC 


FORTRAN  code  Implementing  a point-thinning  technique  Is  Illustrated  by  the 
DO  loop  as  shown.  Input  x-y  data  to  be  plotted  Is  provided  In  the  XLOC  rnd  YLOC 
arrays.  The  AMULT  factor  Is  the  number  of  rasters  per  unit  of  plotting  length, 
c.g.,  200  rasters  per  Inch  for  most  CALCOMP  devices.  The  I FIX  function  converts 
a floating  point  number  to  the  Integer  mode.  The  CALPLT  routine  Is  the  actual 
plotting  routine.  In  Its  call,  the  third  parameter  specifies  pen  up  (a3)  or  pen 
down  (*2). 
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FORTRAN  implementation  of  point-thinning  logic 


TYPICAL  POINT-THINNING  PROGRAM  OUTPUT 


The  point-thinning  logic  has  been  implemented  in  an  independent  plot  proqram 
used  to  plot  from  real-time  analysis  runs  in  the  batch  mode.  In  this  example, 
there  are  twenty-five  data  points  per  second.  Thus  somewhat  over  1000  points  were 
examined  for  each  curve.  Obviously  some  curves  are  more  active  than  others  and 
required  more  data  points  to  be  plotted.  Visual  comparison  with  methods  plotting 
every  point  showed  no  difference  in  curve  representation. 
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RESULTS  AND  CONCLUSIONS 


Typical  results  from  the  real-time  data  plotting  program  showed  that  about 
90  percent  of  the  points  could  be  skipped  In  the  actual  plotting  process.  This 
resulted  in  significantly  reduced  CPU  time  and  reduced  I/O  activity.  Similar 
benefits  are  anticipated  for  interactive  graphics  applications. 
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RESULTS  AND  CONCLUSIONS 


NORMALLY  MORE  THAN  90  PERCENT  OF  THE  DATA  POINTS  COULD 
BE  IGNORED. 

CPU  TIMES  WERE  REDUCED  TO  ABOUT  30  PERCENT  OF  THEIR 
FORMER  REQUIREMENTS  IN  THE  PLOTTING  PHASE. 

RESULTING  PLOT  VECTOR  FILES  (INCLUDING  TITLES,  GRIDS  ETC  ) 
WERE  REDUCED  IN  SIZE  BY  OVER  50  PERCENT. 

THE  METHOD  NEEDS  TO  BE  APPLIED  TO  THE  INTERACTIVE  GRAPHICS 
DISPLAYS.  COMPUTER  RESPONSE  SHOULD  IMPROVE  DUE  TO  BETTER 
CPU  UTILIZATION  AND  REDUCED  PLOT  VECTOR  FILE  SIZE. 
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PARAMETRIC  CUBIC  SURFACE 
REPRESENTATION 


BY 
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NASA  SUBFACE  REPRESENT AT  I ON  UORKSHO  1-2 
AliES  RESEARCH  CENTER  MOFFETT  FIELD, 


t 


HARCH  1378 
CALIFORNIA 


Foil 

Ames  Research  Center  has  developed  a seomet 
oriented  towards  interactive  computer  sra^h 
interface  to  linearized  panel  aerodynamics 
This  system  provides  a aeometm  rer  res-'v 
realistic  aircraft  configurations  from  wh: 
tical  mathematical  models  can  be  created, 
configurations  can  he  assembled  Quick jy  f 
pendent  geometric  components. 


B.  P.  Roland 
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3* 
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r«--  system 
lts»  to 
*-  * os  rams . 
*■  avion  of 
ch  analy- 
Various 
tom  inde- 


GEOMETRY  SYSTEM  GOALS 


* GEOMETRY  DEFINES  PANELLING 

* ASSEMBLE  CONFIGURATIONS  QUICKLY 

* DATA  FORMAT  USEFUL  FOR  INTERACTIVE 
GRAPH  I CS 


Foil 

Parametric  cubic  eauatioris  were  selected  s<-  rie  data 
format  because  they  had  several  important  tea-ures* 

* a single  mathematical  format  can  be  used 
to  represent  curves  and  surfaces  of  el' 
kinds . 

% parametric  eauatioris  do  not  experience 
numerical  difficulties  with  infinite 
slopes . 

* it  is  a mature  technology  with  a larsp 
base  of  software  available  both  in  in- 
dustry and  within  NASA. 


D»  P . Roland 


FEATURES  OF  PARAMETR I C CUB  I CS 


SINGLE  FORMAT  FOR  ALL  CURVES/SURFACES 
HANDLES  INFINITE  SLOPES 

»■ 

MATURE  TECHNOLOGY 

SOFTWARE  AVAILABLE  IN  INDUSTRY 


Foil  ^4 

It  has  been  recognized  that  parametric  curies  have 
some  limitations*  The  ability  to  specify  r Lope  con- 
tinuity in  the  general  case  end  higher  or:.--  deriva- 
tives continuity  was  lacking*  The  maximum  number  of 
sides  in  a patch  is  four  and  some  waviness  occurs  with- 
in a surface.  However?  it  has  been  determined  that 
these  restrictions  will  not  impact  the  reeui r aments  of 
the  aerodynamics  programs  employed. 


D.  F'.  Roland 


LIMITATIONS  OF  PARAMETRIC  CUB  I CS 


NO  GENERAL  SLOPE  CONTINUITY 
NO  HIGHER  ORDER  CONTINUITY 
MAXIMUM  OF  4 SIDES  PER  PATCH 
WAV  I NESS  POSSIBLE  WITHIN  SURFACES 


Foil 


The  geometric  entities  utilized  in  this  system  ere 
the  point*  the  parametric  cubic  curve  and  the  bicubic 
surface  patch.  A point  in  three-space  is  a triple 
of  the  component  coordinates*  A curve  in  three-space 
has  each  coordinate  defined  as  a cubic  polynomial 
of  the  parameter  u.  Similarily  for  a surface*  which 
is  a bicubic  in  u and  w.  Since  the  parameters  are 
limited  to  the  range  0-1*  each  entity  is  bounded 
and  has  a sense  of  direction. 
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PARAMETRIC  CUBIC  GEOMETRY  FNT |. T | ES 
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The  format  of  parametric  cubic  eouations  is  well-suited 


to  matrix  notation* 


In  the  algebraic  form) 


is  com- 
putationally efficient r although  the  coefficients  are 
difficult  to  interpret*  In  he  geometric  form?  the 
coefficients  become  the  end  or  corner  prints-  arid  the 
parametric  tangent  vectors » from  which  the  slope  at  the 


ends  can  be  determined  by  division, 
useful  to  the  engineer  as  a simple 
tion  and  slope  continuity. 


This  form  is  auite 
reference  for  Fosi- 
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CONVENIENT  FORMATS 


* 


P<  u, 
* 

PC  u ) 


ALGEBR I AC  MATRIX  FORMAT  * 


PCU  ) = [u3  U2  U l" 


“>=MO]Mt 


A 

B 

C 

D 


GEOMETRIC  MATRIX  FORMAT  * 
=[f1<  u)  F2<  u)  F3(  u)  F4C  u)] 


F 03  1 #7 

Parametric  cubic  approximations  of  other  commonly 
used  curves  can  be  crested  via  st raishtf orward 
techniRues . A conic  approximation  is  developed  by 
the  beginning  arid  end  points*  slopes  and  heel  point. 
Rho  determines  the  type  of  conic  created!  0 i.  .5  is 
an  ellipse*  .5  is  a parabola  (which  is  exact*  not 
approximate)?  .5  i .75  is  a hyberbola.  The  circle 
approximation  is  a special  case  of  the  conic  having 
the  correct  end  points  and  slopes*  and  the  midpoint 
on  the  circle.  The  graphical  layout  describing  the 
point  of  u = .5  can  be  used  to  Quickly  define  the 
curve  from  its  coefficients  or  to  approximate  even 
hand-drawn  curves.  This  technique  is  alsc  applic- 
able to  space  curves. 


Ti.  P.  Roland 


CREATE  CURVE'S 


CONIC  APPROXIMATION 

P'-  4s*3 

r~4/°0>-n) 

■g 

CIRCULAR  APPROXIMATION 

J°r  0=  4so-te&) 

graphical  layout 

Pf.5)  = + 


£-/? ' 


Genersl  space  curves  can  be  defined  analytically  via 
f^r-^int'  transformation . Co»Plex  shaaas  .r.  <£ 

fined  as  Fiec^ise  continuous  se nen  ♦ segmented 
tinuity  can  be  speciTiea  ir  ieuu'LI': 

curves  with  curvature  continuity  car.  be  created  with 
spline  fittind  techriicaues* 


CREATE  CURVES 


* THROUGH  4 POINTS 

+ 

R, 

* SEGMENTED  CURVES 

POSITIONAL  CONTINUITY  ALLAYS 
SLOPE  CONTINUITY  IF  REQUIRED 


* SPLINE  FIT 

CURVATURE  CONTINUITY 


P(U3) 


A 


Foil  ^ ^ 

Many  aircraft  surfaces  are  specified  analyically. 
This  wind-body  conf iduration  (with/without  canard) 
has  been  used  at  NASA  Landley  for  a series  of  wind 
tunnel  tests i and  it  is  currently  beind  used  to  ver- 
ify the  advanced  panel  aerodynamic  codes  at  NASA 
Ames.  It  has  been  modelled  as  a set  of  components 
utilising  parametric  bicubic  surface  patehs.  The 
surfaces  were  created  from  the  analytical  descrip- 
tion of  the  deometry  in  the  Landley  report. 


D*  F‘.  Roland 


CREATING  SURFACES 

CUfTHOUT  SPLINES} 


< COMPOUND  CURVATURE.) 


Foil  t/0 

In  other  areas*  complex  compound  surfaces  mas  only 
be  defined  as  cross-sections  on  a drawing.  A ten- 
sion-spline  fit  prodram  is  utilized  to  create  patch 
networks  from  digitized  cross  sections  that  have 
positional  continuity  and  regions  of  curvature  con- 
tinuity while  simultaneously  creating  apparent 
slope  discontinuities ♦ 


H.  P.  Roland 
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The  parametric  cubic  form  is  useful  in  mie -active 
graphics  because  a limited  set  of  software  s offices 
for  Geometric  shapes.  The  patch  form  a]]r.~s  each 
entity  to  be  stored  on  disc  as  a separate  record. 
Lerse  numbers  of  patches  need  not  be  kept  ir  core. 
Trar-sf  ormations  on  the  patches  sip  simple  matrix 
multiplications.  The  surfaces  can  have  their  eddes 
connected  and  new  patches  created  between  patches. 


Ii.  F'.  Flo  land 
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USEFUL  IN  INTERACTIVE  GRAPHICS 


* SINGLE  SET  OF  SOFTWARE 

* EACH  ENTITY  HANDLED  SEPARATELY 

« • 

* ABILITY  TO  MODIFY  COMPONENTS 

TRANSLATE,  ROTATE,  SCALE,  SUBDIVIDE 

* ABILITY  TO  CONNECT  COMPONENTS 

CONNECT  BOUNDARIES 
CREATE  RULED  SUFACES 
FILLET  BETWEEN  COMPONENTS 


Foil  #/?• 


The  final  product  of  the  system  is  the  panel  defini- 
tions* Panel  edges  are  created  by  specifying  the 
number  and  distribution  of  parameter  values  at  which 
points  are  to  be  extracted  from  the  parametric  equ- 
ations* Eeual*  sine*  cosine*  half  cosine*  or  user 
defined  spacing  of  the  parameters  may  be  specified* 


Roland 


Foil  J # 

Another  use  of  the  geometric  definition  is  to  create 
the  inputs  to  Level  I aerodynamic  prodrams.  These 
include  surface  areas  and  body  volumes.  Tie  ability 
to  intersect  the  surfaces  with  plane  eouat ions  is 
useful  for  cross  sectional  distributions.  If  the 
pressure  data  is  made  into  a surface  analog*  plane 
intersections  provide  isobars  for  display.  The  inter- 
section between  surfaces  can  be  computed  to  locate 
the  interface  of  components. 


EVALUAT  I NG  SURFflPF<? 


SURFACE  ( WETTED ) AREA. 

VOLUME 

4f  • 

PLANE  INTERSECTIONS 
CROSS  SECTIONAL  AREA.  D I STR I BUT  I ON 
PRESSURE  SURFACE  ISOBARS 

SURFACE-SURFACE  INTERSECTIONS 
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In  conclusion*  parametric  cubic  geometry  is  a useful 
technique  for  defining  realistic  aircraft  confidura- 
tions  for  linear  panel  aerodynamic  methods*  The 
limitations  do  not  currently  impact  the  reeui cements 
and  improved  methods  will  he  developed* 
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CONCLUS  I ONS  - PARAMETRIC  CUBIC  GEOMETRY  i £ ; 


* USEFUL  NOW 

* capable  of  defining  realistic  conf i gurat i ons 

* MEETS  THE  REQUIREMENTS  OF  CURRENT  ACTIVITIES 

* IMPROVED  METHODS  ARE  COMING 


fioCOAJO 


SURFACE  1 


INTERSECTION  CURVE 


G(u,v,s,t)  * rr(u,v)  - r2(s,t)  * 0 
THREE  EQUATIONS  IN  FOUR  UNKNOWNS 

SOLUTION  HAY  BE  EXPRESSED  AS  r(o),  OR  EQUIVALENTLY  {u(a) ,v(o) ,s(a) WHERE 
a IS  SOME  ARC  LENGTH  RELATED  PARAMETER 


y 'mm&Z  2j 


HUNTING  PHASE 


• THE  HUNTING  GRID  IS  ESTABLISHED  ON  THE  PRIMARY  SURFACE  BY  Ng. 

• EACH  LINE  IN  THE  GRID  IS  A CURVE  DEFINED  BY  EITHER  u OR  v BEING  CONSTANT. 

• THE  INTERSECTION  OF  A CURVE  WITH  THE  SECONDARY  SURFACE  IS  FOUND  BY  SOLVING  THE 

CONSISTENT  SET  OF  EQUATIONS  £(v,s,t)  = 0 OR  H(u,s,t)  = 0. 

• THE  SOLUTIONS  OF  THESE  SETS  OF  EQUATIONS  PROVIDE  THE  INITIAL  CONDITIONS  FOR  THE 
TRACING  PHASE. 


CURVE-SURFACE  INTERSECTION  CALCULATION 


CURVE  CRAWLING  PROCEDURE  TESTS  . Vi+1  = * 

THE  POINT  ON  THE  SURFACE  (X)  APPROXIMATES  THE  CLOSEST  POINT  ON  THE  SURFACE 
GIVEN  POINT  ON  THE  CURVE  (•). 


NOTE:  MINUS  SIGNS  INDICATE  A REVERSED  DIRECTION  FOR  THE  SEGMENT 


1 . Generation  of  Input  Data  for  a Three-Dimensional  Potential -Flow  Proora^ 

This  talk  describes  work  performed  at  Douglas  Aircraft  Company  order 
contract  to  NASA  Lar.gley  to  develop  a geometry  package  to  s''~ol  ify  the 
task  of  preparing  the  input  data  for  a potential -f lev;  program.  Tr.e  par- 
ticular program  for  which  the  geometry  package  was  developed  is  the  three- 


dimensional  lifting  version  of  the 


-ann  program,  which  was 


developed  by  John  Hess  under  Navy  sponsorship  and  whico  Is  ir.  use  at 
numerous  companies,  universities,  and  government  agero'es.  T'"e  geometry 
package  is  sufficiently  general,  however,  to  be  of  use  ir.  ct^er  applications. 
A major  portion  of  the  expense  involved  *r‘  using  the  notent'ai -flow  program, 
is  consumed  by  the  man-hours  required  to  prepare  tne  input  data  (virtually 

all  rx-f  u/hir-h  ir  nr-nmni-vi  r-  rl  -» 4-  ->  > 


Hal sey/Hess 
Slide  1 


<c 

QC 

o 

UL_ 

<c 

I — 

<c 

a 


LU 

o 


o 

> — c 

h- 

<c 

DC 

UJ 

sc 

LU 

CD 


h- 

o 

CL. 


I 

LU 

LU 

dc 

zrz 


s: 

<c 

DC 

CD 

O 

DC 

Q_ 


O 

! 

LL. 


i 


Typical  Element  Distribution  for  a Wing-Fuselage  Case 

Preparing  the  geometric  input  is  so  time-consuming  because  the  coter, 
program  requires  a large  number  of  points  for  good  accuracy.  A t.-pi 
fuselage  case,  shown  here,  required  approximately  800  points. 


n 


WING-FUSELAGE  CASE 
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elements  used  by  bell  aerospace  for  a 

FLUSH  HATER  INLET  ON  A SURFACE 
EFFECT  SHIP  (-3500  ELEMENTS) 


DUCT 

CROSS  SECTION 


mm 

fciifcB 

wmit 
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ACCOUNTED  FOR 
BY  SYMMETRY  — 


rr";/  -A4 
7 ''///, 
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4.  Factors  Which  Influence  Point  Spacing  Requirements 

In  addition  to  shear  number  of  points,  the  user  must  take  care  in  the  way 
he  distributes  the  points.  On  a wing , for  example,  the  proximity  of  leading 
or  trailing  edges,  the  tip  or  root,  breaks  in  the  planform  or  any  other 
corners,  all  influence  the  point  spacing  which  should  be  used.  Cn  more 
general  bodies,  regions  of  high  curvature  or  any  factor  which  causes  the 
solution  (either  the  source  density  or  the  velocity)  to  vary  rapidly,  have 
to  be  considered. 
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factors  which  influence  point 


SPACING  REQUIREMENTS 


0 

0 

o 

o 

o 

o 


LEADING  EDGES 
TRAILING  EDGES 


CORNERS 

PLA.NFORM  BREAKS 


AREAS  OF  HIGH  CURVATURE 

STRF>jfTU°noRAPID  variation  01 
RENGTH  OR  SURFACE  VELOCITY 


SINGULARITY 


5.  Definition  of  Frequently  Used  Terms 

In  addition,  the  user  must  organize  the  points  in  a man 


■er  consistent  with 


the  logic  of  the  potential -flow  program.  .his  os  ids  illustrates  sorr.^  c-f  ~ne 
loaical  considerations  and  also  defines  some  terms  to  be  used  i«ter  th«s 
talk. 


cener; 


; or : 

con~i 


.ration,  seen  as  ar 


rcraft. 


is  divided  ir. 


c one 


noer 

^ r 


of  simpler  components,  such  as  a winq  or  a fuselage,  or  a porti 
these.  Each  component  is  defined  by  specifying  points  on  a section  curve, 
then  on  an  adjacent  section  and  so  on,  until  all  sections  have  been 
defined.  Points  on  a section  curve  are  ordered  as  if  one  were  traversing 
the  perimeter  of  the  section  (not  generally  fore-and-aft).  Such  section 
curves  are  designated  N-lines.  The  curves  connecting  correspcnding  points- 
on  the  N-lines  are  designated  M-lines.  The  area  bounded  by  two  adjacent 
N-lines  is  termed  a strip  and  the  area  bounded  by  two  N-lines  and  two  M-lines 
is  termed  an  element. 
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STRIP  ON 
INTERSECTED 
COMPONENT 


STRIP  ON 
ITERSECTING 
COMPONENT 


.EMENT  ON 
ERSECTING 
COMPONENT 


N-LINES  ON 

INTERSECTED 

COMPONENT 


M-LINES  ON 

INTERSECTED 

COMPONENT 


ELEMENT  ON 
INTERSECTED 
COMPONENT 
N-LINES  ON 
INTERSECTING 
COMPONENT 


M-LINES  ON 

INTERSECTING 

COMPONENT 
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6.  General  Features  of  the  Geometry  Package 

This  geometry  package  simplifies  the  user's  job  in  several  ways.  First, 
it  greatly  reduces  the  number  of  points  for  which  the  user  must  specify  f~e 
coordinates.  It  is  often  possible  to  reduce  this  number  by  an  order  of 
magnitude,  or  more,  from  the  number  required  for  an  accurate  potential -flow 
solution.  Second,  it  reduces  the  care  that  he  must  devote  to  spacing  the 
points,  since  it  provides  a number  of  automatic  spacing  options.  Th'rd,  the 
very  difficult  task  of  calculating  intersection  curves  (if  any)  between  the 
components  is  done  automatically.  Finally,  if  there  are  intersections,  sr~e 
repaneling  is  done  to  obtain  compatible  paneling  distributions  between 
- adjacent  components. 


i 
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GENERAL  FEATURES  OF  THE  GEOMETRY  PACKAGE 


o INITIAL  GEOMETRY  DATA  CAN  BE  VERY  SPARSELY  DEFINED 

o GEOMETRY  DATA  AUGMENTED  AND  REDISTRIBUTED  (SEVERAL 
OPTIONS) 

o INTERSECTION  CURVES  CALCULATED  (USING  P,  C,  PATCHES) 

o FINAL  REPANELING  MAKES  ADJACENT  ELEMENTS.  LINE  UP 
ACROSS  INTERSECTION  CURVES 
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7.  Use  0*  the  'Geometry  Package  for  a Wing-Fuselage  Configuration 

This  is  a typical  example  of  v/hat  the  geometry  package  can  do.  Because  of 
the  simplicity  cf  the  configuration,  a very  small  number  of  points  is  suf- 
ficient to  define  the  geometry.  The  geometry  package  has  added  ar.d  redis- 
tributed points  in  both  chordwise  and  spanwise  directions  on  the  wing  and 
in  both  axial  and  circumferential  directions  on  the  fuselage.  It  has  also 
. found  the  wing-fuselage  intersection  curve  and  matched  the  spacing  cf  the 
points  across  the  intersection  curve. 


ENRICHED  OUTPUT 
DISTRIBUTION 


8.  Paneling  of  Isolated  Components 

The  first  operation  performed  by  the  geometry  package  is  the  paneling  of 
isolated  components.  In  this  operation,  the  program  ta?.es  the  sparse  input 
data  and  augments  it  to  a specified  number  of  points  distributed  according 
to  one  of  several  possible  spacing  algorithms.  At  this  stage  intersections 
are  completely  ignored;  each  component  is  paneled  as  if  it  were  completely 
independent  of  the  others  (except  that  proximity  of  other  components  should 
be  considered  before  deciding  on  the  desired  spacing  algorithm  and  a desire 
for  continuity  between  adjoining  components  may  dictate  that  these  adjoin- 
ing components  be  paneled  using  the  same  algorithm) 


This  slide  shews  a wing 


various  stages  of  the  panel -'ng  operation. 


Initially,  it  is  defined  tv  a small  rum 
(tip  and  root).  The  program  first  dist 
N- lines  and  then  adds  more  N-l ir.es.  Tr. 
by  redistributing  points  cm  the  "-lines 
cedure  is  accomplished  by  two  sets  of  i 


her  cf  oc’nts  sn  ; u 
.ributes  points  or.  t 
.e  adding  of  N-l  inc-s 
. Tr.us,  the  entire 
nterpolaticns  on  ou 


st  two  .1-1  .res 
r.e  initial 
is  acoc"ol is^ed 
paneling  pre- 
rva: ; no  general 


surface-f i tting  is 


required. 


INITIAL  ELEMENT  DISTRIBUTE 


9.  Interpolation  Procedure  for  General  Curves  in  Space 

The  interpolation  procedure  is  slightly  unusual.  Cubic  splines  were 
initially  used,  but  with  the  very  sparse  defining  data  (which  generally 
may  not  have  been  smoothed)  the  waviness  could  get  out  of  control.  The 
method  finally  chosen  can  be  very  accurate,  but  it  is  very  simple  to 
implement  and  less  likely  than  a true  spline  to  cause  problems.  This  is 
because  a slight  waviness  introduced  in  one  portion  of  a curve  cannot 
propagate  to  other  portions  of  the  curve. 

The  independent  variable  of  the  interpolation  is  the  straight-line  arc  length 
between  the  defining  points  (normalized  by  the  total  length  of  the  curve). 
When  normalized  in  this  manner,  the  straight-line  arc  length  is  a very  good 
approximation  to  the  true  arc  length,  which  is  perhaps  the  most  general 
quantity  tnat  can  be  chosen  as  the  independent  variable.  Separate  calcula- 
tions are  performed  to  interpolate  x,  y,  and  z coordinates  versus  arc 
length.  In  each  interpolation,  the  first  derivatives  of  the  dependent  var- 
iable (x,  y,  or  z)  with  respect  to  straight-line  arc  length  are  first  found 
by  taking  a weighted  average  of  the  angles  of  the  straight-line  segments. 
These  are  then  used,  together  with  the  coordinate  values,  to  determine  the' 

. coefficients  of  cubic  interpolating  polynomials. 

10.  Comparison  of  Curve-Fit  Methods 


INTERPOLATION  PROCEDURE  FOR  GENERAL  CURVES  IN  SPACE 


o THE  STRAIGHT-LINE  ARC  LENGTH  BETWEEN  INPUT  POINTS 

»»  >»E  TOTAL  LENGTH)  IS  T H E , » D E P E « D E ■ T 

0 SEPARATE  INTERPOLATIONS  ARE  USED  FOR  (X  VS  S) 

(Y  VS,  S)  AND  (Z  VS.  S). 

o PROCEDURE  FOR  EACH  INTERPOLATION; 


0 


FIND  THE  1ST  DERIVATIVES  OF  THE  FUNCTION  AT  THF 
INPUT  POINTS  BY  TAKING  A WEIGHTED  AVERAGE  OF  THE" 
ANGLES  OF  THE  STRAIGHT-LINE  SEGMENTS 

0 FIT  A CUBIC  CURVE  OVER  EACH  SEGMENT  (GIVEN  2 
POINTS  AND  2 FIRST  DERIVATIVES  ON  EACH  SEGMENT) 

THIS  IS  NOT  A TRUE  SPLINE  METHOD  SINCE  THE  2ND 
DERIVATIVES  ARE  DISCONTINUOUS, 

RESULTS  ARE  CONSISTENTLY  LESS  WAVY  THAN  RESULTS  OF 
A TRUE  SPLINE  METHOD  (PERHAPS  B E C A U S E T H E I HD E P *HD 
cNT  VARIABLE  IS  HOT  A CONTINUOUSLY  VARYING  QUANTITY) 
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1 1 . Options  for  Spacing  of  Points  on  N-Llnes 

This  slide  shows  the  available  options  for  the  spacing  of  points  on  tf-lines 
( i . e . , in  the  churdwise  direction  on  wings  and  usually  the  circumferential 
direction  on  fuselages).  The  first  option  is  to  make  no  change.  The  second 
and  fifth  options,  input  distribution,  augmented  in  number  and  curvature- 
dependent  distribution,  are  described  further  in  later  slides.  The  third 
and  fourth  options,  constant  increments  in  arc  length  and  the  so-called 
* cosine  spacing  (or  constant  increments  on  a superscribed  circle),  are  very 
common  and  need  no  further  explanation.  The  sixth  option  allows  the  user 
to  specify  any  arbitrary  arc-length  distribution  he  desires. 
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OPTIONS  FOR  S P A C I N fi  OF  POINTS  ON  N - L I N E S 
(CHORD  WISE  PANELING  ON  A 1/ 1 N fi  ) 


INPUT  DISTRIBUTION,  UNALTERED 

INPUT  DISTRIBUTION,  AUGMENTED  IN  NUMBER 


CONSTANT  INCREMENTS  IN  ARC  LENGTH 
COSINE  SPACING 


CURVATURE-DEPENDENT  DISTRIBUTION 
USER-SPECIFIED  DISTRIBUTION 


12.  Curvature-Dependent  Soacng  **ethgd 

A method  in  which  the  point  spacing  is  a function  of  the  local  surface 
curvst'ire  should  be  useful,  since  the  results  of  the  potential -flow  method 
shou  d vary  rapidly  in  regions  of  high  curvature.  In  this  method,  the 
spacing  is  linearly  dependent  cn  the  absolute  value  of  the  curvature  and 
some  control  over  the  variation  of  the  spacing  ever  fe  curve  is  provided. 
Artificial  curvatjre  is  added  in  the  generally  flat  trail ing-edge  regions  on 
lifting  sections  in  order  to  bunch  points  there,  as  recuired  Dy  the  potential 
flow  method.  ~he  implementation  of  tne  method  requires  an  iterative  procedure 
in  which  the  arc  length  /alues  are  initially  estimated,  the  curvature  at  the 
. center  of  each  segrent  of  the  curve  is  calculated  and  tnen  used  (with  the 
specified  equation;  to  update  tne  estimated  arc-iengtn  aistrioution. 
Convergence  generally  reouires  only  three  or  four  iterations.  Examples  of 
the  method  for  a lifting  and  a nonlifting  section  are  shown  in  the  slide. 


* * 
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CURVATURE-DEPENDENT  SPACING  METHOD 

Asi  * ^0  ” Asmir/Asmax^  “ k/kmax^  + Asmir/Asmax^  Asmax 
where 

as  « ARC  LENGTH  BETWEEN  ADJACENT  POINTS 
k = CURVATURE  (AT  ELEMENT  MIDPOINT) 


WHEN 

ki  " kmax  ’ 

AS . = AS. 
i min 

WHEN 

k,  = 0 

AS.  = As 
• 1 max 

i * 

SPECIFY 

4smin/is.,v  = 0.25 
min  max 

as.  - [0.75(1  -ki/kmax)  + 0,25]  Asmax 
0 SOLVE  FOR  AS,  BY  AN  ITERATIVE  PROCEDURE 


13.  Method  of  Enriching  Number  of  Points  While  Maintaining  a Similar  Distribution 

It  is  impossible  to  provide  an  automatic  point  spacing  method  that  will  be 
appropriate  for  the  large  variety  of  cases  which  could  be  considered.  This 
method  provides  a great  deal  of  flexibility,  however.  The  user  must  input 
properly  distributed  points,  but  he  need  not  input  an  extremely  large  nurr.cer 
of  points. 

To  enrich  the  number  of  points  while  maintaining  a similar  distribution  to 
the  input  distribution,  it  is  necessary  to  define  a normalized  point  number 
(as  shown)  and  construct  a curve  of  arc  length  (at  the  input  points)  as  a 
function  of  this  normalized  point  number.  Arc  lengths  at  the  output  points 
are  determined  simply  by  interpolating  on  this  curve  to  find  the  values 
* corresponding  to  the  output  values  of  the  normalized  point  number.  An 
example  of  a section  of  a supercritical  wing  enriched  from  fifteen  to 
twenty- five  points  is  shown. 


NORMALIZED  ARC-LENGTH  - (S/P) 


Spacing  nf  n.i 

^uslfa P°-‘nfS  °n  N‘lines 

more  explanation,  however/  dr"sectlon  "ode  and  arc-length 


disfibutio"  Of  P0i..„  „„ 

d re,Wre  more  explanation 


* 


OPTIONS  FOR  SPACING  OF  N-LINES 
(SPAN  WISE  PANELING  ON  A WING) 

1.  INPUT  DISTRIBUTION,  UNALTERED 

2..  INPUT  DISTRIBUTION,  AUGMENTED  IN  NUMBER 

3.  CONSTANT  INCREMENTS 

'I.  USER-SPECIFIED  DISTRIBUTION 

A . PLANAR-SECT  ION  MODE 


B.  ARC-LENGTH  MODE 


Halsey/Hess 
Slide  14 


15.  Spanwise  Redistribution  of  Elements  on  a Supercri tical  Wing 

N-lines  on  this  wing  have  been  distributed  using  the  planar-section  mode  of 
operation.  In  this  mode,  points  on  all  N-lines,  except  possibly  the  first 
and  last  N-lines,  lie  in  parallel  planes.  The  specified  distribution  param- 
eters refer  to  distances  between  the  planes,  rathe'*  than  to  arc  lengths. 

This  option  is  important  for  lifting  components,  such  as  wings,  since  the 
potential -flow  program  requires  elements  on  lifting  components  to  be 
trapezoidal.  (Elements  on  ncnlifting  components,  such  as  fuselages,  may  be 
more  general  quadrilaterals.) 

16.  Comparison  of  Planar-Sertinn  xnH  irr.ion"^  *■ — ' — -- 


SPA.'IWISE  redistribution  of  elements  on 

a s U P E R C R I T I C A L WING 
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16.  Comparison  of  Planar-Section  and  Arc-Length  Modes  of  Distribution  of  N-Lines 

In  some  cases,  such  as  this  strut  on  a thick  wing,  planar  distributions  of 
N-lines  leave  undesirable  sparse  areas.  This  is  prevented  by  use  of  the  afc- 
. length  mode  of  distributing  N-lines.  In  this  mode,  the  specified  distribution 
parameters  are  the  arc  lengths  along  the  M-lines  and  the  distribution  procedure 
is  identical  to  the  procedure  for  distributing  points  on  N-lines. 


COMPARISON 


OF  PLANAR-SECTION  AND  ARC-LENGTH 

OF  distribution  OF  N — lines  — 

STRUT  OH  a thick  wing 


MODES 


(a)  Planar-Section  Mode. 


(b)  Arc-Length  Mode. 


Assumptions  of  r.he  Intersection_^thol 


the  next  step  performed  by 


Having  paneiea  me  , . ..  n rho  rurvos  of  inwr«i.u«ii  \ • 


Having  paneled  the  cOTponent^as  Holated^o^es,^^  intersection  (if 
jmetry  pac<a' 
any) . This  is  not  a 
are  nut  too  restrictive  for  the  typical  cases  encountered  in  aerodynamic 
applications.  First,  it  is  assumed  that  a distinction  can  be  made  between 
intersecting  and  intersected  components.  The  M-lines  of  the  intersecting 
component  pierce  the  elements  of  the  intersected  component.  These  v-lines 
must  extend  completely  through  the  intersected  component  and  must  not  blend 
in  tangent  to  the  surface  of  the  intersected  component.  A component  may 
intersect  only  one  other  component  and  may  be  intersected  by  only  one  other 
component.  Intersecting  and  intersected  components  are  identified  in  the 
input  data. 


Steps  in  Calculatinq  Intersection  Curve 


An  intersection  curve  is  defined  as  a series  of  intersection  points  between 
the  M-lines  on  the  intersecting  component  ^nd  the  surface  of  the  intersected 
component.  For  each  M-line  a search  is  conducted  to  determine  wnich  segment 
of  the  M-line  intersects  which  element  of  tre  intersected  component.  For 
the  purposes  of  this  search,  each  element  is  represented  by  two  triangular 
planes,  but  each  M-line  segment  is  assumed  to  be  curved.  When  the  intersected 
element  is  found,  its  geometric  data  is  used  to  determine  its  parametric 
cubic  patch  coefficients.  These  are  not  stored  permanently  because  of  the 
large  amount  of  redundancy  in  storing  P.C.  coefficients  of  patches  with 
cormon  boundaries.  The  final  calculation  of  the  intersection  point  reouires 
the  solution  of  three  simultaneous  nonlinear  equations  and  is  described  in 
the  next  slide. 


STEPS  IN  CALCULATING  INTERSECTION  CURVE 


FOR  EACH  INTERSECTING  M - L I N E : 

1.  SEARCH  FOR  THE  INTERSECTED  ELEMENT 

2.  FIND  APPROXIMATE  INTERSECTION  POINT 
(ASSUMING  PLANAR  ELEMENTS) 

5 ' F 1 M D P * C • PATCH  COEFFICIENTS  OF  ELEMENT 

''  • pointLATE  f1  0 R E P R E C 1 S E INTERSECTION 
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1 9 . Intersection  of  a Cubic  Curve  and  a Parametric  Cubic  Surface 

Knowledge  of  the  P.C.  patch  coefficients  of  an  element  allows  the  equation 
of  the  surface  to  be  written  in  terms  of  the  parametric  variables,  u and 
w,  as  shown  in  this  slide  (equations  for  y and  z coordinates  are  of 
exactly  the  same  form  as  the  equation  hown  for  the  x coordinates).  Equa- 
tions for  the  x,y,  and  z coordinates  of  points  on  the  cubic  curve  are 
. also  known.  Equating  coordinate  values  on  the  curve  and  on  the  surface 
gives  three  nonlinear  equations  for  the  three  unknowns  (u,  w,  and  s). 
bolution  by  Newton's  method,  starting  from  an  approximate  solution  obtained 
in  the  searching  operation,  generally  converges  in  less  than  about  six 
iterations. 


INTERSECTION  OF  A CUBIC  CURVE  AND  A 

CUBIC  SURFACE 


PARAMETRIC 


o REPRESENTATION  OF  SURFACE: 

xs  = w3(Axu3  + Bxu2  + Cxu  + Dx) 

* "2(Ex“3  + Fxu2  + Gx“  + V 

+ w (Ixu'1  * v2  * V * s,) 

♦ (M  u3  + H u2  + 0 u + P ) 

x x x x' 

o REPRESENTATION  OF  CURVE: 

*c  = Qxs3+  Rx$2  + Sys  + T 
^ x x x X 

0 3 EQUATIONS:  x - v = n 

S xc  0 

ys-yc  - 0 
*S-2C  = 0 

o 3 UNKNOWNS:  / . 


(AND  SIMILAR  EXPRESSIONS  FOR  Y AND  Z 
COORDINATES) 


o 


SOLVE  BY  NEWTON'S  METHOD 
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20.  Illustration  of  Intersection  Method 

Various  analytic  cases  (involving  spheres,  cones,  cylinders,  and  ellipsoids) 
have  been  used  to  verify  the  accuracy  of  the  intersection  method . Various 
cases  involving  realistic  aircraft  components  have  also  been  calculated. 

Few  of  these  cases  have  sufficient  character  to  be  interesting  to  view 
graphically.  A less  realistic  case  that  does  have  more  character  is, 
therefore,  shown  in  this  slide. 


* 


. use  or  tne  ueometry  package  tor  a wmg-mse laqe  configuration 

After  calculating  the  intersection  curves,  sore  repaneling  takes  place. 
Various  options  for  different  sorts  of  cases  are  provided.  In  the  interest 
of  brevity,  only  the  wing-fuselage  case,  oerhacs  the  most  common  apolication 
of  the  geometry  package,  will  be  described.  Most  of  the  calculations  in  this 
portion  of  the  method  involve  procedures  very  similar  to  those  already 
described  (for  example,  interpolating  along  curves  or  the  intersection  of 
curves  and  planes). 

First,  the  intersection  curve  between  the  wing  and  the  fuselage  is  made 
an  N-line  on  the  wing  and  all  N-lines  outboard  of  this  are  shifted  to 
maintain  a smooth  distribution  of  N-lines.  The  area  of  the  wing  inside 
the  fuselage  becomes  a single  strip  which  functions  as  a means  of  making 
the  vorticity  continuous  between  left  and  right  wings.  Planar  N-lines  on 
the  fuselage  are  passed  through  the  leading  and  trailing-edge  points  on  the 
intersection  curve  and  both  N-  and  M-lines  are  shifted  cn  the  forward  and 
the  after  portions  of  the  fuselage  in  order  to  maintain  smooth  point  distri- 
butions. N-lines  on  the  fuselage  are  passed  through  points  on  the  inter- 
section curve  ( :*her  through  every  point  or  every  other  point)  and  then 
points  on  these  N-lines  are  redistributed  to  provide  a smooth  distribution. 
The  area  covered  by  the  wing  is  not  paneled.  This  final  repaneling  breaks 
up  a fuselage  (assumed  here  to  initially  consist  of  only  one  component) 
into  four  components  — one  forward,  one  aft,  one  above,  and  one  below  the 

wino 


* * 
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PANEL  2:  METHODS 

1.  What  are  the  unsolved  (or  not  reported  in  the  literature)  motl  -matical 
p rob  lens: 

(a)  Intersections  and  blending  of  arbitrary  splined  surfaces? 

(b)  Filleting  with  arbitrary  radius? 

(c)  Design/specification  of  surface  front  known  aerodynamic 
solutions  or  pressure  distributions? 

2.  Which  method  of  surface  representation  are  most  conservative  of  computer 
storage  and  processing  power? 

3.  Can  we  list  the  advantages  and  di sadvantages  of  each  method?  There  are 
several  bi products  inherent  in  each  approach.  Can  we  summarize  these7 

*».  Will  "connectivity  tables"  that  represent  the  relative  location  of  various 
aircraft  components  become  cumbersome? 


\ 


•/ 
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DISCUSSION  VUGRAPH  1. 


This  paper  discusses  an  analytical  method  for  computing 
a body  fitted  coordinate  system  for  an  arbitrary  three-dimensional 
flow  field.  This  research  has  been  carried  out  at  General  Dynamics, 
Fort  Worth  Division  under  a continuing  NASA/Ames  Contract  to  extend 
General  Dynamics  2-D/axisymmetric  finite-difference  flow  field 
computation  procedure  to  three  dimensions. 


A BODY  ORIENTED  MESH-GENERATION 


TECHNIQUE  FOR  3-D  FLOW  FIELDS 


BY 


Ishwar  C.  Bhateley 
& 

Leroy  L.  Presley 


PRESENTED  AT 

1)  AIAA  6TH  Mmisimposium.  Arlington, Texas  25  February  1978 

2)  NASA  Workshop  on  Aircraft  Surface  Representation  for  Aerodynamic 
Computation,  March  1 2, 1978 


B1 t 124 

~SJ-/Ar€csy  / 


DISCUSSION  VUCRAPH  2. 


Extensive  research  has  been  done  in  developing  two- 
dimensional  body-fitted  coordinate  systems  at  Mississippi  State 
University  and  the  University  of  Cincinnati.  Two  transformation 
techniques  for  mapping  an  arbitrary  doubly  connected  r gion  into 
a rectangle  are  shown  in  this  slide.  In  the  first  technique  the 
inner  and  outer  boundaries  of  the  doubly  connected  region  map  into 
two  opposite  boundaries  of  the  rectangular  region.  An  arbitrary 
curve  connecting  the  two  boundaries  in  the  physical  plane  maps 
into  the  other  two  sides  of  the  rectangle  in  the  transformed 
plane.  In  the  second  technique  the  inner  boundary  is  mapped 
into  a horizontal  or  vertical  slit  in  the  transformed  plane, 
while  the  exterior  boundary  is  broken  into  four  arcs  each  of 
which  maps  into  a side  of  the  rectangle  in  the  physical  plane. 

The  Laplace  equations  are  used  as  the  transformation  functions. 

The  functions  on  the  right-hand  side  of  the  equations  can  be 
chosen  such  as  to  provide  desired  coordinate  system  control. 


TWO  DIMENSIONAL  METHODS  FOR  GENERATION 
OF  BODY  FITTED  COORDINATES 


PHYSICAL  PLANE 


TRANSORMED  PLANE 


( 


B 1 1 osc 
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DISCUSSION  VUGRAPH  3. 


The  first  transformation  technique  described  for  two- 
dimensional  doubly  connected  regions  could  not  be  readily 
extended  to  arbitrary  three-dimensional  doubly  connected 
volumes.  However,  the  slit  transformation  technique  can  be  extended 
to  three  dimensions.  In  this  approach,  the  inner  surface 
transforms  to  a planar  slit  parallel  to  one  of  the  coordinate 
planes, while  the  exterior  surface  is  divided  into  six  parts, 
each  of  which  maps  into  a side  of  a rectangular  solid  in  the 
transformed  field.  It  is  important  that  each  part  of  the 
exterior  surface  and  interior  surface  be  approximately  in 
the  same  relative  position  in  the  transformed  plane  as  in  the 
physical  plane  to  obtain  reasonable  cell  distribution.  Again, 
the  Laplace  equations  are  used  as  transformation  functions. 

The  functions  on  the  right-hand  side  can  be  used  to  provide 
coordinate  system  control.  The  inverse  transformation 
equations  for  three  3-D  flows  have  been  derived  and  will  be 
published  soon  as  the  contractors  report.  Numerical  solu- 
tions are  obtained  using  finite-difference  approximation  to 
the  various  partial  derivatives  and  successive  over  relaxation 
iteration. 


THREE  DIMENSIONAL  SLIT  TECHNIQUE  FOR  GENERATION 

OF  BODY  FITTED  COORDINATES 


PHYSICAL  PLANE  TRANSFORMED  PLANE 


DISCUSSION  VUGRAPH  4. 


The  application  of  the  3-D  slit  transformation  technique 
to  generate  a body  fitted  cell  arrangement  for  an  inlet  flow 
field  is  shown  in  this  slide.  The  curved  surface  of  the 
cylindrical  exterior  boundary  is  divided  into  four  parts  each 
of  which  map  into  four  sides  of  a rectangular  volume  in  the 
transformed  field.  The  two  circular  faces  of  the  cylindrical 
exterior  surface  map  into  the  other  two  sides  of  the  rectangular 
volume.  The  inlet  is  also  divided  into  four  segments  each  of 
which  map  into  a rectangular  planar  slit  parallel  to  the  coordinate 
planes.  These  rectangular  areas  are  connected  and  themselves 
describe  a rectangular  box.  Care  is  taken  to  partition  the 
external  boundary  and  the  inlet  in  such  a manner  that  the  same 
relative  position  of  the  various  surface  components  is  maintained 
in  the  physical  and  transformed  fields.  The  volume  enclosed 
by  the  inlet  in  the  physical  field  maps  to  a volume  enclosed 
by  the  planar  slits  in  the  transformed  field. 


DISCUSSION  VUGRAPH  5. 


Since  the  cell  definition  on  both  the  exterior  and  interior 
boundaries  serves  as  a boundary  condition  for  the  generation  of 
3-D  mesh  arrangement,  the  two-dimensional  slit  method  was  used  to 
obtain  a cell  definition  for  the  upstream  circular  face  of  the 
exterior  cylindrical  boundary.  The  circular  face  boundary  was 
divided  into  four  equal  a^eas,  each  of  which  were  mapped  into 
a side  of  a square  in  the  transformed  plane.  The  coordinate 
system  obtained  by  this  process  is  shown  in  this  slide. 
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DISCUSSION  VUGRAPH  6. 


The  starting  solution  on  the  downstream  face  of  the 
exterior  volume  is  also  generated  using  the  two-dimensional 
slit  transformation  technique.  The  circular  exterior  boundary 
is  mapped  into  four  sides  of  a square  in  the  transformed  plane. 
The  annular  inlet  cross-sectional  area  is  divided  into  four 
parts  each  of  which  map  into  four  connected  slits  parallel  to  the 
coordinate  axis  forming  an  embedded  square  as  shown  in  the  slide. 
The  resulting  solution  shown  in  this  slide  was  used  as  a starting 
solution  for  3-D  cell  generation  program 


STARTING  SOLUTION  FOR  INLET  FLOW  FIELD 

DOWNSTREAM  FACE 

2 D SLIT  TRANSFORMATION  TECHNIQUE 
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DISCUSSION  VUGRAPHS  7 THROUGH  13. 


Typical  starting  intermediate  and  converged  solutions  for 
the  3-D  circular-symmetric  inlet  flow  field  using  the  tr.thod 
discussed  previously  are  shown  in  the  following  seven  vugraphs. 
Cell  arrangement  for  four  transverse  and  three  lateral  cuts 
are  shown.  No  coordinate  system  control  was  used.  These  solu- 
tions were  generated  using  the  GD  interactive  graphics  facility. 

A 47  x 31  x 31  mesh  arrangement  was  generated.  As  can  be  seen, 
a satisfactory  cell  structure  was  generated.  The  cells  from 
one  cut  to  the  next  were  blended  with  no  sharp  discontinuities. 
Methods  to  improve  the  cell  arrangement  by  incorporating 
coordinate  system  control  and  redistribution  of  points  on  the 
boundaries  is  being  investigated.  The  extention  of  this 
technique  to  generalized  3-D  flow  fields  is  also  being  undertaken. 


COMPUTED  MESH  GEOMETRY  FOR  SYMMETRIC  CIRCULAR  INLET 

WITH  ELLIPTIC  LIP  SHAPE 


COMPUTED  MESH  GEOMETRY  FOR  SYMMETRIC  CIRCULAR  INLET 

WITH  ELLIPTIC  LIP  SHAPE 
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Intermediate  Solution  C15  iterations) 


Convereed  Solution  OO  iterations) 
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Title  Page: 

A primary  protlem  in  Computational  Aerodynamics  is  obtaining  flew  field  solutions 
about  irregular  geometries.  An  effective  approach  to  this  protlem  is  to  transform 
the  governing  equations  and  boundary  conditions  zo  a coordinate  system  where  -he 
protlem  is  most  easily  attacked.  For  regular  geometries  there  are  analytical 
transformations  which  are  most  appropriately  used,  however,  for  the  general  case 


numerical  approaches  car.  be  employed. 


BOUNDARY  FITTED  CURVILINEAR  COORDINATE 


SYSTEMS  USING  TENSION  SPLINE 


FUNCTIONS 


Figure  1: 


Consider  the  physical  domain  for  a simple  two-dimensiooal  aerodynamics  Frotlan  s.Jch 
as  the  flow  about  an  airfoil  as  shown  at  the  top  of  this  figure.  Assume  that  it  is 
desirable  to  numerically  obtain  the  solution  of  the  fluid  flow  about  the  airfoil  using 
a rectangular  unifonn  grid  where  the  boundary  of  the  airfoil  r,  transfonns  to  the  unit 
inter s-0  on  the  computational  domain  and  the  free  stream  boundary  Fg  transforms  to 

the  Unit  ir‘terVal  3=1  °n  the  COmpUtational  domain.  The  cut  in  the  physical  domain  r 
and  r4  transforms  to  the  unit  intervals  at  t=0  and  w on  the  left  ar.d  right  sides  of 


the  computational  domain. 


In  order  to  solve  a physical  problem  using  the 


computational 


domain  there  must  be  a relationship  between  each 
ar.u  each  point  (t^  in  the  computational  domain 
between  corresponding  points  must  be  known. 


point  (x^y^  in  the  physical  domain 
. AlsOj the  differential  relation 


Figure  2: 

There  is  a set  of  transformation  equations  relating  the  physical  domain  to  the 
computational  domain.  For  a doubly  connected  region  in  two-dimensions  the  trans- 
formation equations  ar.d  boundary  conditions  are  shown  in  this  figure.  The  equations 
are  two  nonlinear  coupled  elliptic  partial  differential  equations  relating  the 
computational  domain  to  the  physical  domain.  A technique  for  numerically  solving 
this  system  of  equations  is  found  in  reference  1.  The  approach  is  to  select  a set 
of  points  on  each  boundary  which  are  to  be  connected  by  grid  lines,  choose  initial 
guesses  of  all  the  grid  points  in  the  physical  domain  corresponding  to  grid  points 
in  the  computational  domain,  and  select  forcing  functions  F and  G that  will  yield 
the  desired  concentration  of  grid  points  . This  is  followed  by  the  application  of 
Successive  Over— Relaxation  of  the  discretized  partial  differential  equations  'until 
convergence . Initial  guesses  to  the  nodes  in  the  computational  domain  are  necessary 
because  of  the  nonlinearity,  and  these  guesses  should  be  relatively  close  to  the 
desired  converged  va'  *s . 
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Figure  2 TRANSFORMATION  EQUATIONS 
2-0  D°ubly  Connected  Region) 


Figure  3: 


The  approach  presented  herein  is  algebraic  and  is  based  on  parametric  cubic  poly- 
nomial equations  ar.d  tension  spline  approximation  to  the  boundary  data.  The 
technique  works  for  singly  and  doubly  connected  regions  in  two-dimer. si on s and, 
it  can  be  < tended  to  three-dimensions.  Since  the  process  is  algebraic,  it  is 
computationally  fast  and  amenable  to  interactive  computer  graphics.  This  tech- 
nique is  a special  case  solution  to  the  general  transformation  equations . The 
primary  functions  are  shown  in  this  figure.  The  algebraic  transformation  equations 
which  relate  t and  s to  x and  y are  parametric  cubic  polynomials  in  the  variable  s 
and  depend  on  eight  parameters  which  are  functions  of  the  variable  t.  (see 
reference  2 . . The  parameters  are  position  and  derivatives  with  respect  to  s on 
the  boundaries  r±  ar.d  rg.  These  p-.ameters  are  obtained  by  tension  spline  fitting 
to  sets  of  data  defining  ^ • Jd  rg.  The  relationship  between  derivatives  with 
respec.  to  s ar.d  t is  shown  at  the  bottom  of  the  figure.  The  variable  t is 
effectively  the  percentage  of  accumulated  cord  length  on  the  two  boundaries. 


Figure  3 PARAMETRIC  CUBIC  TRANSFORMATION  2-D 
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Figure  4: 

The  tension  spline  functions  are  represented  symbolically  in  this  figure.  The 
independent  variable  (t  or  t)  of  the  spline  fits  to  the  boundaries  is  adulated 
cord  length  along  the  boundaries.  The  variable  t «Xt<L)  is  the  percent  of  cord 
length.  The  spline  functions  consist  of  sets  of  piecewise  continuous  hyperbolic 
functions  and  the  parameters  ^ and  ^ are  tension  parameters  which  govern  the 
degree  of  damping  of  the  fit.  Reference  3 describes  the  tension  spline  functions. 


Figure  4 r'ARAMETRIC  SPLINE  APPROXIMATION 
TO  BOUNDARY  DATA  WITH  TENSION 


Figure  5: 

This  figure  illustrates  the  reasons  behind  using  tension  spline  functions.  Initially 
cubic  splines  were  used,  but  wiggles  occurred  in  some  of  the  fits  such  as  near  the 
leading  edge  of  the  airfoil  shown  in  the  figure.  The  tension  spline  routine  that  has 
been  used  degenerates  to  a cubic  spline  when  the  tension  factor  is  near  sero.  Increasi 
the  tension  factor  to  a large  number  in  effect  increases  the  damping  which  forces  the 


ng 


fit  between  data  points  to  be  nearly  straight  lines  as  illustrated  at  the  bottom  of 
the  ligure.  Note  that  the  tension  factors  are  user  chosen. 


Figure  6: 

*“■  flgUre  m“£trateS  the  fits  to  the  boundary  data  start. 

Points  on  the  two  born,**  - * ±ng  with  th-  data 

•e  two  boundarres  and  iodicattng  the  initial  points  cord  1.  tw 

-d  tension  spline  tits  to  the  two  sets  of  data  are  perked 

a-  to  he  closed  nor  of  i*e  see  arc  length.  On  T * “* 

anywhere  on  the  fits  can  he  interpolated,  ^ ^ Perf°™d'  ^ 

°°rd’  P°SitiODS  “ - and  derivatives  with  respect  ^ 

Equation  (2)  is  Bm,v  . P to  t can  be  evaluated. 

(2.  is  applied  to  compute  any  point  between  ttecw.. 
value  of  t.  aen  the  two  boundaries  for  the  given 
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Figure  J: 

This  figure  illustrates  the  transformation  from  a uniform  rectangular  coordinate  system 
to  a coordinate  system  about  a Karman  Trefftz  airfoil.  The  technique  has  worked  as  well 
for  other  airfoils  including  supercritical  airfoils.  Also,  it  has  worked  for  singly 
connected  domains  in  two  dimensions.  It  is  noted  again  that  the  technique  is  a special 
case  solution  of  the  general  transformation  equations,  however,  it  is  algebraic  and 
computationally  fast  yielding  dense  data.  Also^the  use  of  the  technique  is  simple  as 
compared  to  numerically  solving  the  differential  equations. 
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Figure  8: 

The  density  of  the  grid  lines  near  a boundary  can  be  controlled  by  the  transformation 
shown  in  this  figure.  Large  negative  values  of  k concentrate  grid  lines  near  Pg. 
Large  positive  values  of  k concentrate  grid  lines  near  r1-  When  solving  the  differ- 
ential system  (equation  l)  this  is  controlled  by  the  forcing  functions. 
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Figure  8 CONTRACTION  OF  THE  GRID  NEAR  A BOUNDARY 


Figure  9: 


The  effects  of  the  parameter  k are  shown  in  this  figure.  A value  of  k = 5 is  used  with 
the  Karxnan  Trefftz  airfoil.  Four  normal  grid  lines  are  shown  corresponding  to  F = .05, 
.1,  .15,  and  .2.  At  the  bottom  of  the  figure  one  grid  line  is  shown  for  s = .05  and 
k = 2.  There  are  other  ways  of  controlling  the  grid  line  in  both  the  t and  s directions 
while  At  and  As  remain  constant,  however,  this  illustrates  the  ease  of  adding  such 


transformations . 


BLOWUP  OF  LEADING  EDGE 


K-5 


Figure  10: 

Thir  figure  illustrates  the  proposed  extension  to  three-dimensions . The  parametric 
cubic  equations  are  exactly  like  those  on  figure  3 except  for  the  addition  of  one 
equation  for  the  third  coordinate  and  the  replacement  of  the  variable  S by  W.  The 
position  and  derivative  parameters  are  now  obtained  from  surface  definition  of  the 
:,iimer"  and  "outer"  boundaries . The  derivatives  with  respect  to  W are  obtained 
from  the  cross  product  relation  and  U and  V correspond  to  percent  of  cord  in  the 
longitude  and  latitude  directions. 
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FIGURE  10  EXPANSION  TO  THREE  DIMENSIONS 


Figure  11: 


This  figure  illustrates  how  the  inner  and  outer  boundaries  aught  be  conceived . 
Although  the  complexity  of  the  inner  boundary  may  be  over  optimistic,  each 
component  part  such  as  wing,  fuselage,  or  nacelle  could  be  used  to  generate  a 
three-dimensional  grid.  The  boundaries  presented  here  were  generated  with  the 
computer  program  "A  Computer  Program  for  Fitting  Smooth  Surfaces  to  an  Aircraft 
Configuration  and  Other  Three-Dimensional  Geometries"  described  in  reference 
Positions  and  derivatives  are  available  on  both  boundaries . 


Conclusions : 

Although  it  is  not  anticipated  that  this  technique  will  supplant  the  solving  of  the 
differential  transformation  equations,  it  offers  a simple  and  rapid  solution  for  a 
transfoimation  in  two-dimensions  for  a large  nuriber  of  cases  and  is  extendable  to 
three-dimens ion s . 


CONCLUSIONS 


1.  Curvilinear  coordinate  systems  can  be  generated 
algebraically  for  singly  and  doubly  connected  regions 
in  2-  D using  parametric  cublcs  and  tension  spline 
functions. 

2.  Density  of  grid  lines  near  a boundary  is  easily 
controlled. 


3.  The  technique  should  be  extendable  to  three-dimensions. 
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SLIDE  3 

The  intent  of  the  approach  is  to  not  only  provide  an  algorithm 
which  is  intuitively  correct  but  is  also  systematically  correct 
Hence,  here  we  attempt  to  obtain  a mathematical  representation 
not  only  accurate  relative  to  the  data,  but  also  correct  in  its 
mathematical  form. 
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SLIDE  4 


This  "goodness  of  fit"  criteria  (the  Tchebycheff  Criteria)  allows 
us  to  approximate  with  this  error  in  mind.  We  attempt  to  deduce 
a mathematical  representation  constrained  by  minimizing  this 
maximum  relative  error. 


SLIDE  5 1 
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Shown  here  is  a typical  pressure  profile  for  a NACA  wing  section. 

Typically,  the  point  to  be  noted  is  that  the  pressure  data  is  | 

recorded  to  four  digits.  Hence,  we  intuitively  deduce  that  this  | 

data  is  accurate  to  four  significant  digits.  Secondly,  we  note  ] 

the  rapidly  changing  pressure  gradient  from  the  stagnation  point  ! 

at  the  leading  edge  and  likewise  on  the  trailing  edge.  Now,  we 

want  to  deduce  a mathematical  representation  which  is  only  | 

constrained  by  such  a relative  error,  i.e.,  accurate  approximately  ] 

to  four  significant  digits. 


SLIDE  6 


Four  basic  mathematical  forms  are  considered:  the  polynomial 
fonn.  the  rational  form  (i.e.  ratio  of  two  polynomials),  and 
therr  piecewise  counterparts:  piece., ise  polynomials,  piecewl: 
rational  and  with  or  without  smoothness  constraints  imposed. 
Likewise,  the  four  basic  mathematical  forms  in  two-dimensional 
space  are  the  double  polynomials/rationals/wlth  or  without 
smoothness  constraints,  i.e. 

f<x,y)  - t T «tc, 

j«i  i-i 
for  example. 


SLIDE  7 

Analgously,  these  basic  forms  can  be  represented  in  terms  of 
Tchebycheff  Polynomials.  Vie  perform  our  computation  in 
terms  of  such  expansions  for  the  purpose  of  maintaining 
computational  accuracy.  This  advantage  is  due  to  their 
orthogonality  property. 
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ITVr  U Sh0Kn  th'  TC',e,,VChef,  Po'^-'s  expanded  in 
terns  of  the  ordinary  polynomials  and  vice-»ers..  The  detail 
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SLIDE  9 


Shown  on  this  slide  is  the  pievewise  construction  of  smoothness  at 
assumed  knots  x . If  we  have  a mathematical  representation  y^(x) 
for  a set  of  prescribed  data  over  the  first  subinterval  expressed 
as  shown,  we  can  impose  the  value  y(l)(X  ) for  the  second  approxi- 
mation at  its  left-hand  end  point.  Likewise,  the  value  y^^(X  ) is 
imposed.  Inductively,  we  can  do  so  for  any  such  subinterval.  What 
we  get  is  a over-determined  coupled  system  of  linear  equations  which 
we  eliminated  for  a j simultaneously.  Likewise,  the  two-dimensional 
counterparts  follow  analogously. 
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The  Tchebycheff  Problem  restated  i„ 
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Now  since  in  general,  we  are  interested  in  not  only  determining 
one  dimensional  curves  f(x)  but  also  two  dimensional  surfaces 
f(x,y),  we  have  arrived  at  an  apparent  impasse  because  the 
classical  Tchebycheff  theory  d 3s  not  apply.  However,  if  we 
reformulate  our  Tchebycheff  problem  in  terms  of  the  Linear 
Programming  context,  we  are  able  to  solve  our  problem.  The 
opposite  slide  shows  this  reformulation. 


SLIDE  12 


Now  computationally,  we  do  not  wUh  , , 

,,  not  wish  to  solve  the  primal  linear 

rcrr-  bai rather  ,ts  *"•  «. 

des.reabl.  because  th.  number  of  constraints  is  typically  much 

less  in  the  dual  formulation.  Also  we  seek  to  H r 

order  of  crh  seek  t0  detennine  the  least 

order  " apprOX,matio"-  ™*  Possible  since  the  least 

~ ‘S  »*—  « - kth  approximation  by  the  «t 
t at  on.  additional  constrain  (k-1  th)  can  be  added  to  the  d„, 

simplex  tableau  and  we  can  resume  the  oral-  . 

an  resume  the  ordinary  simplex  operations. 
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Some  numerical  results  are  shown  on  the  opposite  slide.  The  test  case  is 
y * 1 - X2.  Seven  data  points  are  inputted  with  a 10%  error.  These  data 
points  could  reflect  measured  data  with  some  error  and  are  marked  on  the 
chart  with  X's.  The  exact  solution  is  the  solid  curve.  The  dashed 
curve  is  a cubic  spline  fit.  The  absolute  Tchebycheff  solution  is 
depicted  by0.  The  relative  Tchebycheff  solution  is  approximately  on 
top  of  the  exact  solution.  The  results  show  waviness  of  the  cubic  spline 
solution  due  in  part  to  interpolation  thru  the  data  points.  The  absolute 
Techebycheff  approximation  failed  near  the  roots  of  y - 1 - X2  and  was 
affected  by  the  10%  error  resulting  in  a .12  maximum  residue.  The 
relative  Tchebycheff  approximation  not  only  showed  no  waviness,  but 
was  the  most  accurate  (1.2%  maximum  relative  error).  It  also  was  able 
to  identify  the  curve  y « 1 - X*.  That  is,  it  characterized  it  by 
determining  it  was  in  fact  c*  pirabola  by  determining  its  coefficients 

y * .987623  - .987623X2.  Finally,  the  approximation  was  able  to  reduce 
the  effect  of  the  error. 
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A major  objective  at  NASA-Ames  Research  Center  is  to  develop  the 
ability  to  integrate  aerodynamic  theory  with  experiment.  Wind 
tunnel  test  results  and  theoretical  predictions  of  the  aerodynamic 
configurations  will  be  compared  on  the  local  minicomputer  system. 

In  order  to  accomplish  this  objective,  it  is  necessary  to  develop 
geometry  models  that  are  as  detailed  and  as  accurate  as  the 
physical  wind  tunnel  model,  in  order  to  develop  such  a mathematical 
model,  the  Interactive  Parametric  Equation  Geometry  System  (IPEGS) 
was  developed. 
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The  IPEGS  System  can  be  broken  down  Into  six  major  steps: 

STEP  1 : Digitize  Drawings 

The  cross  sections  of  a particular  component  of  the  wind  tunnel  model,  e.g.,  the  tail, 
nose,  or  upper  fuselage,  are  digitized  using  an  optical  digitizer.  The  digitized  points 
of  the  cross  section  are  displayed  as  they  are  being  digitized  on  an  IMLAC  CRT.  This 
display  ensures  that  the  operator  doesn't  digitize  a bad  point  and  also  that  he  has 
sampled  the  cross  section  sufficiently  to  get  its  representative  shape.  Some  of  the 
components  of  the  model,  e.g.,  the  wing,  canard,  or  vertical  tail,  are  not  digitized 
but  are  input  analytically  into  the  IPEGS  System. 

STEP  2:  Create  Surfaces 

The  digitized  points  are  transmitted  to  the  local  PDP-11  minicomputer  where  they  are 
transformed  into  parametric  bicubic  surfaces  through  the  use  of  tension  splines. 

STEP  3:  Review  Surfaces 

The  parametric  bicubic  surfaces  are  then  examined  on  the  Evens  & Sutherland  Picture 
System.  The  picture  system  allows  the  operator  to  rotate,  translate,  or  scale  the  object 
in  all  3 dimensions.  The  operator  can  also  display  the  object  in  four  views  simultaneously, 
reflect  the  object  about  any  axis  or  display  cross-sectional  views  of  the  object. 

STEP  4:  Modify  Surfaces 

The  operator  can  interactively  translate  or  scale  the  entire  object  or  any  component  of 
the  object.  In  this  way  he  can  easily  exchange  components  of  the  model  or  modify  any 
component.  He  can  also  "pick"  any  bicubic  surface  or  patch  and  then  operate  on  that 
surface.  He  has  the  ability  to: 

(a)  Split  a patch  into  two  patches. 

(b)  Delete  a patch. 

(c)  Force  two  patches  to  connect  together  with  or  without  slope  continuity. 

(d)  Create  a fillet  patch  between  two  existing  patches. 

(e)  Create  a ruled  surface  patch  between  two  existing  patches. 

( Co*V ) 

J.  Ashbaugh 
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STEP  5-  Create  Networks 


After  the  model  has  been  reviewed  and  modified  to  the  specifications  of  the  aerodynamist, 
the  paneling  information  required  by  a particular  aerodynamic  program  can  be  extracted. 

The  distribution  (sine,  cosine,  half  cosine,  even  spacing)  and  the  density  of  the  paneling 
can  be  changed  interactively  to  emphasize  the  critical  areas  of  the  model. 

STEP  6:  . Review  Output 

The  paneling  information  is  sent  to  the  CDC  7600,  operated  on  by  the  aerodynamics  codes 
and  the  output  plot  information  is  sent  back  to  the  PDP-11  and  the  EfiS  Picture  System. 

The  output  plots  can  be  quickly  scanned  on  the  Picture  System  and  a hardcopy  made  of  any 
plot. 
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OBJECTIVES 

o ACCURATE  MODEL  - GEOMETRY  DEFINED  BY  MATHEMATICAL  EQUATIONS. 
SURFACE  COORDINATES,  SLOPES  AND  CURVATURES  ARE  CALCULATED. 

o GENERATE  MODELS  AT  ALL  LEVELS  OF  CONFIGURATION  DEFINITION  STARTING 
WITH  THE  INITIAL  “BACK  OF  THE  ENVELOPE'  SKETCH  THROUGH  THE  FINAL 
THREE-VIEW  DRAWING. 

o CAPABILITY  OF  ADDING  DETAIL  TO  THE  MODEL  AS  THE  DESIGN  IS 
DEVELOPED. 

o COMMON  GEOMETRY  - SINGLE  SOURCE  OF  GEOMETRY  INPUTS  FOR  A 
VARIETY  OF  PROGRAMS. 
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Interactive  Input  For  The  QUICK  Geometry  System 


In  order  to  compute  the  flow  around  any  body  in  detail,  the  body  surface 
description  must  be  sufficiently  smooth  to  avoid  generating  disturbances  that 
would  not  occur  on  the  prototype.  Also,  many  of  the  methods  for  flow  comput- 
ation require  points  on  the  surface  to  be  defined  without  restrictions  imposed 
by  the  geometry  method.  These  requirements  can  be  met  by  defining  the  surface 
analytically.  An  added  benefit  of  analytic  geometry  definition  is  that  it 
allows  derivatives  of  the  surface  contour  to  be  determined  analytically  and 
therefore  exactly. 

The  QUICK  geometry  system5  fills  that  need  for  an  analytic  surface 
definition  method  for  a wide  range  of  moderately  complex  aircraft  geometries. 

It  has  been  applied  to  such  codes  as  a supersonic  shock-fitting  finite  difference 
method2  and  a transonic  wing-body  flow  code15.  A system  for  generating  the 
inputs  to  QUICK  interactively,  using  a graphics  terminal  connected  to  a time- 
sharing computer  system,  will  be  described.  When  fully  developed,  this  system 
will  make  QUICK  much  easier  to  use  and  therefore  more  readily  accessible  to 
anyone  requiring  its  capabilities. 


1.  Vachris,  A.  F.,  Jr.;  and  Yaeger,  L.  S.:  QUICK-Geometry  - A Rapid  Response 

Method  for  Mathematically  Modelling  Configuration  Geometry.  NASA  SP-390, 
1975,  pp.  49-61. 

2.  Marconi,  F. ; and  Yaeger,  L.  S.:  Development  of  a Computer  Code  for  Calculating 

the  Steady  Super/Hypersonic  Inviscid  Flow  Around  Real  Configurations. 

NASA  CR-2676,  1976. 

3.  Boppe,  C.  W.:  A Computational  Method  for  Transonic  Wing-Body  Flows.  AIAA 

Paper  78-104,  January  1978. 
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In  the  QUICK  geometry  system  concept  the  aircraft  surface  is  enveloped  by  a series 
of  body  lines.  Each  of  these  lines  is  a mathematically  defined  curve  in  space,  con- 
sisting of  a sequence  of  linked  curve  segments  (generally  conics).  The  intersection  of 
these  body  lines  with  any  desired  cross-section  plane  defines  a set  of  control  points 
in  that  plane.  Elliptical  arcs  fitted  to  these  control  points  according  to  a logically 
defined  cross-section  model  determine  the  surface  shape  at  that  axial  location. 


*4f 


CONCEPT  OF 


CONFIGURATION 


BODY- LINES 

MATHEMATICAL  DEFINITION 


LOGICAL  DEFINITION 


/OUJ/US£K)b 


z 


A principle  barrier  to  the  use  of  QUICK  has  been  the  difficulty  in 
understanding  the  concept  and  relating  it  to  the  required  program  inputs. 
Especially,  having  to  begin  the  inputs  with  logical  definitions  of  cross 
section  models  has  turned  some  people  away  without  giving  the  program  a 
fair  trial.  Once  the  concept  has  been  mastered,  further  difficulties  arise 
in  trying  to  accurately  match  any  even  moderately  complex  configuration. 

Some  of  these  difficulties  are  related  to  the  geometric  limitations  of  QUICK 
itself,  particularly  the  requirement  that  the  surface  be  single-valued  in 
polar  coordinates.  But,  aside  from  these,  there  is  often  difficulty  in 
finding  the  appropriate  locations  of  body  lines  controlling  surface  slopes 
or  in  choosing  which  of  the  many  possible  shapes  for  a body  line  segment 
gives  the  best  surface  fit  with  the  desired  configuration.  Making  these 
choices  often  comes  down  to  an  iterative  situation  for  which  batch  mode 
operation  is  too  cumbersome  to  allow  a sufficient  number  of  trials  to 
completely  determine  the  optimum  model. 

What  is  needed  is  a new  mode  of  operation  which  will  immediately  display 
the  results  of  a choice  graphically  to  allow  its  evaluation,  which  will  allow 
new  choices  to  be  made  interactively  with  the  computer  as  required,  and  which 
will  lead  the  user  through  the  process  of  making  choices  until  the  whole 
configuration  has  been  designed.  This  mode  of  operation  would  also  do  much 
to  avoid  the  barrier  to  conceptual  understanding  of  QUICK. 
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The  concept  being  developed  steins  from  the  need  to  work  interactively  with  the 
computer  using  a graphics  terminal  for  display  and  user  input.  It  will  generate  an 
input  file  (or  card  deck)  in  the  formats  required  for  the  inputs  to  QUICK.  It  will 
operate  within  the  geometric  limitations  of  QUICK  (using  the  same  equations  and  sub- 
routines where  possible)  so  that  the  resulting  configuration  should  be  acceptable  by 
the  QUICK  system. 

In  this  concept  the  control  point  locations  in  the  cross  sections  are  defined 
(numerically)  as  the  cross  section  models  are  being  defined  (logically).  These  con- 
trol point  locations  from  the  cross  sections  are  then  "strung  together, " using  the 
QUICK  curve  segments  interactively,  to  define  the  body  lines.  A data  base  system 
(SPAR4)  is  used  for  mass  storage. 


4.  Giles,  Gary;  and  Haftka,  Raphael:  SPAR  Data  Handling  Utilities.  Proposed  NASA  TP 
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These  next  few  slides  were  made  directly  from  the  screen  of  the  graphics  terminal. 
They  show  some  features  of  the  interactive  input  concept  as  it  has  been  implemented  so 
far. 


Slide  5 shows  the  first  displays  on  the  screen.  The  circles  indicate  user  responses. 
"Restart"  provides  for  the  option  to  continue  working  on  a configuration  previously 
started  but  not  completed.  The  list  of  cross  sections  below  the  response  "3“  refers  to 
input  cross  sections  to  be  matched.  These  are  sets  of  (y,z)  surface  points  obtained  by 
digitizing  from  drawings  or  by  taking  cross  section  cuts  through  some  other  surface 
description  (e.g.  Harris  inputs).  They  are  previously  stored  on  the  data  base.  There 
might  be  a hundred  of  these,  but  only  the  four  used  for  the  check  case  are  shown. 

QUICK  requires  as  input  for  the  cross-section  logical  definitions  Hollerith  control 
point  names,  which  later  become  Hollerith  body-line  names  in  the  body  line  definition 
phase  of  input.  The  24  control  point  names  shown  were  pre-selected  so  that  the  user  can 
refer  to  them  by  number,  rather  than  by  typing  names  into  the  key  board.  The  model  names 
help  the  user  keep  track  of  which  cross-section  logical  definitions  have  been  made.  The 
arc  shapes  and  types  are  also  Hollerith  inputs  required  by  QUICK  and  referred  to  by 
number. 
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Slides  6A,  6B,  and  6C  show  three  stages  in  the  logical  definition  of  a cross  section 
model.  The  cross  section  to  be  matched  was  obtained  from  Harris  inputs  for  a supersonic 
fighter  proposal.  The  cursor  is  used  to  locate  the  map  axis  (center  for  polar  represen- 
tation of  cross  section  by  QUICK)  and  each  control  point  as  it  is  called  for.  A list  of 
control  points  used  is  displayed  at  the  upper  right.  After  the  user  has  located  the  two 
or  three  control  points  needed  to  define  an  arc,  he  enters  three  numbers  corresponding  to 
the  arc  type,  the  arc  before  and  the  arc  after  (0  if  none).  The  arc  so  defined  is  immedi- 
ately displayed  on  the  screen  to  be  accepted  or  rejected.  Since  the  arcs  are  drawn  using 
the  same  equations  as  QUICK  uses,  any  satisfactory  arc  will  be  satisfactory  when  done  by 
QUICK;  conversely,  if  an  arc  is  unsatisfactory  it  would  not  be  done  correctly  by  QUICK  and 
must  be  done  again  with  changes  in  control  point  location  or  type.  When  all  of  the  arcs 
have  been  defined,  the  arc  numbers  are  input  in  order,  thus  completing  the  logical  discre- 
tion of  the  cross  section  model.  The  completed  cross  section  is  shown  as  slide  6C.  Note 
that,  several  arc  types  were  used  in  its  definition. 

The  above  process  has  not  only  logically  defined  the  model  but  has  also  located  the 
control  points  at  one  x-station.  The.*;  locations  will  later  be  used  to  define  the  body 
lines. 
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Slide  8 shows  the  same  cross  section  as  on  the  previous  slide.  Because 
of  the  reverse  curvature  it  was  impossible  to  match  the  top  of  this  cross 
section  using  the  same  cross  section  model  logical  definition.  Therefore 
a new  logical  definition  was  made  with  two  additional  arcs.  Note  that 
arc  7 was  not  done  correctly  the  first  time  and  had  to  be  done  over.  On 
a graphics  terminal  having  a refresh  capability  the  erroneous  arc  would  be 
deleted  so  as  not  to  clutter  the  screen.  Note  also  that  the  last  two  arcs 
could  have  been  done  over  by  moving  the  control  points  slightly  so  as  to 
represent  the  desired  surface  more  closely. 
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Slide  9 shows  the  proposed  method  of  defining  body  lines  to  complete  the 
geometric  definition  of  the  confiquration.  As  was  noted  previously  the 
control  point  locations  defined  in  doing  the  cross  sections  are  "strung 
together"  using  the  QUICK  curve  segments.  The  side  view  (xz  plane)  and 
top  view  (xy  plane)  of  each  body  line  is  defined  separately.  "Aliasing" 
refers  to  the  QUICK  provision  for  defining  a body  line  as  exactly  matching 
a previously-defined  body  line  when  this  occurs  rather  than  re-doing  it. 
"Scaling"  is  a provision  for  changing  the  vertical  scale  of  the  plot  on  the 
screen  to  accommodate  the  variety  of  body  lines  which  may  occur  for  some 
configurations. 
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After  all  body  lines  have  been  defined,  all  the  information  collected 
will  be  output  in  the  formats  required  for  input  to  QUICK. 

Experience  so  far  has  indicated  that  the  system  being  developed  will 
be  easy  to  learn  to  use,  even  by  those  who  have  never  used  QUICK.  The 
results  shown  indicate  that  with  care  good  matches  can  be  made  with 
moderately  complicated  cross  sections. 

The  program  is  being  written  using  ANSI  standard  FORTRAN,  and  is  being 
made  machine-independent  as  much  as  possible  in  order  to  enhance  its 
portability.  (The  hardware  being  used  are  a PRIME  400  computer  and  a 
TEKTRONIX  4014  graphics  terminal  with  interactive  buffer.)  The  bodyline 
part  of  the  program  is  still  being  written,  and  the  whole  program  will 
continue  to  be  developed  as  experience  is  gained  with  it. 
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FUTURE  SPACE  TRANSPORTATION  STUDIES 


Since  the  early  seventies,  the  Space  Systems  Division  has  been  studying  advanced  space  trails* 
portation  systems  to  define  critical  technology  areas  which  need  to  be  developed  in  order  to 
direct  present  funding  plans.  In  order  to  evaluate  the  impact  of  a technology  advancement  on 
the  total  vehicle  system,  all  of  the  technical  disci  pi  ines— structures,  propulsion,  subsystems, 
aero  thermodynamics,  and  cost— must  be  integrated  into  a complete  design  synthesis.  Due  to  the  many 
concepts  being  studied,  a general  computer-aided  design  system  was  developed  to  handle  the  analysis. 
Aerodynamics  are  an  integral  part  of  the  design  process  since  aerodynamic  surface  mass  comprises 
approximately  25  percent  of  the  total  vehicle  dry  mass.  This  percentage  can  vary  greatly  depending 
upon  center-of-gravity  position  and  operational  mode  of  the  vehicle. 
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future  space  transportation  concepts 

space  transportation  excepts  have  ranged  fn»  single-stage  _ _ , 
concept  with  twin  turbojet  haters  and  , rocket  seconds  *°  * 
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speed  repine  basically  for  hypersonic  trlw  and  L/Oband  s ha  a in9  "alwtl  throuph  the 
landing  speed.  A main  problaa  area  Is  predict!  " SUb,mf‘  *rtm,  t/D,  and  design 

«rtng  conbinatlons.  '*  Subs,w1c  sUfc,1,tJ'  for  the  large  bluff  body 


3 


W.  WILHITE 


FUTURE  SPACE  TRANSPORTATION  CONCEPTS 
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HISTORY  OF  COMPUTER-AIDED  DESIGN 


The  development  of  , new  colter-elded  design  (CAD)  system  et  Langley  was  initiated  due 
to  , acic  of  capacity  of  any  one  present  method.  0,  to  the  variation  of  concepts-, aonch 
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HISTORY  OF  COHPUTER  AIDED  DESIGN 
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avid  system 
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AEROSPACE  VEHICLE  INTERACTIVE  DESIGN  (AVID)  SYSTEM. 
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MENU  OF  AVID  PROGRAMS 


TO  execute  any  program,  the  user  only  has  to  h,t  the  appropriate  hey  at  the  terminal 
At  the  end  of  a design  session,  the  data  base  and  gentry  data  can  he  saved.  At  a later 
time,  the  design  can  be  retrieved  for  future  analysis  In  any  technical  area. 
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ALA*  W.  WILHITE 


avid  programs 

INPUT  LETTER  OF  PROGRAM  TO  EXECUTE 

A - RESTORE  PREVIOUSLY  SAVED  SIMULATION 
B - DIGITIZE  BODY  SHAPE 


D - 

E - 
F - 

G - 
H - 

I - 

J - 

K - 
L - 

M - 

N - 

0 - 
P 


rLUi  BODY  SHAPE  (AWW  IMAGE) 

VOLUMES  AND  AREAS  (WAB) 

HYPERSONIC  AERODYNAMICS  (AWW  NEWTONIAN! 
HYPERSONIC  AERODYNAMICS  (HA^CR 

SUBSONIC  AERODYNAMICS  (DATCOM) 

SUBSONIC  AERODYNAMICS  (VORTEX  LATTICE) 
PROPULSION  SYSTEM  CHARACTERISTICS 
TRAJECTORY  CALCULATION  (REHDER  MINI-TRAJ) 

MASS  BREAKDOWN  (MARTIN 'TASK  II  BASELINE) 
MASS  BREAKDOWN  (MARTIN  TASK  II  W/CG) 

INTERACTIVE  DATA  BASE 
COST  (JAM  VERSION  OF  WILCOX) 

LIFE  CYCLE  COST  (JAM) 

SCREEN  REPORT  (GEOMETRY  SIZED) 


Q - SAVE  THIS  SIMULATION 
R ’ SEQUENCE  MODE 
S - ENDS  EXECUTION  OF  AVID 


I 


* i..  .1 


i • j *#t 


ALAN  W.  WILHITE 


i. 
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BODY  DIGITIZING 


*">  - «.  1’wTw”  Z'nZnT  “ “f ■"“* 

:r::r„rrr:rLz 

to  give  the  desired  cross-sectional  creature  To  inte  te"Si0"  *"  t,9hte,,ed  or  ,00s«*< 

case  with  ce  cross  section  both  the  how  n’°  between  cross  sections  or  for  the 

section,  both  the  body  planter,  and  side  view  are  digitized. 


ALA*  W.  WILHITE  , 


MINS  DIGITIZING 

For  the  wfnq  the  planfom  Is  digitized  , 

standard  NACA  or  arhltrary,  are  Input  hy  TTY;  the  J^U  T ^ 

mot  airfoil  thickness  In  order  to  oia«  po1nt  of  tb*  ■*rtn« 

Horizontal  and  vertical  surfaces  are  handled'  snarly!*  " dt"eetl0B  ^ “*  ten,,”,  “rS#r- 


i V.  ... 

1 -.1.  , /*>V 


' 13  *'  ( ( V'- 


M.  WILHITE' 


WING  DIGITIZING 


* 

INPUT  CR  - INPUT  IS  CORRECT 
RD  - RE-DIGITIZE 
GM-  GLOBAL  MOVE 
PM-  POINT  MOVE 
GS- GLOBAL  SCAIE 
fX*X  SCAIE 
PY:Y_SCAI£ 


I 62.53799 

TTY  INPUT 

i- 

• INCIDENCE 

• AIRFOIL 

' • STANDARD  NACA  ' 

• SPECIFY  ABSCISSA  AND  ORDINATES 

* 


* 


I 

1 
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4LMI  M.  WILHITE 


VEHICLE  MODELING 


Both  the  externa,  and  interna,  gentry  can  he  modeled  with  the  1nterlctt 

system.  The  externa,  geometry  was  modeled  with  S ' <"fr,cti,e  peo^try 

j ras  moaeiea  with  3 cross  sections  I nrat«H 

of  the  payload  bay,  and  wing  j„„ct, on.  This  externa,  conf,  „ the  nose,  start 

1"  .ess  than  ,5  minutes.  rf  every  point  on  th T “*  - «*•«- 

be  generated  in  , to  2 hours  ,„Z  ""  dl‘rft,"d-  ' “"«**"«•"  ««,d 

nours.  Internal  geometry  is  diqitized  in 

external  geometry.  The  internal  geometry  is  used  for  t * *""* 

and  rocket  engine  placement.  r " ^ arran9e"ent»  volu"e  allocations. 
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AW*  W.  WILHITE 


AERODYNAMIC  PROGRAMS 


AVID  is  mostly  used  for  conceptual  and  preliminary  design.  Very  fast  programs  (>5  sec) 
which  use  only  a gross  definition  of  the  vehicle  are  used  for  conceptual  systems.  For  prelimi- 
nary design,  programs  which  use  X,  Y,  Z coordinates  are  used  (>5  min).  The  more  detailed 
programs  have  not  yet  been  integrated  into  the  AVID  system. 
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ALAN  W.  WILHITE 


AERODYNAMIC  PROGRAMS 


% 


DESIGN  LEVEL 

LEVEL  1 
(CONCEPTUAL) 

LEVEL  II 
(PRELIMINARY) 

LEVEL  III* 
(DETAILED ) 

PROGRAMS 

• DATCOM  1 

• DATCOM  II 

• HYPER 

• SKIN  FRICTION 

• HABACP 

• VORTEX  LATTICE 
t WAVE  DRAG 

• STEIN 

• HESS 

• WOODWARD 

- ‘ 

INPUT 

• GROSS 
• sref.  _A _b. 
-A-te’ 

• NUMERICAL 

• HARRIS 

• X,  Y,  Z 

• NUMERICAL 
AND  LOFTING 

* NOT  AVAILABLE  WITH  AVID 
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ALAN  V.  WILHITE 


PROGRAM  INTERFACING  TECHNIQUES 


transfect  ’nterfaCi"9  Pr0Sra,S  W°  "»*  ^ ^ ««r«t  data 

transfer  to  the  prog™.  This  method  Is  the  «,st  efficient  but  Is  not  used  for  ae™,^e 

17^"!  ” t0  ^ C,”P,eXity  °f  ^ inPUt-  ^ ^ "eth°d  — • «.  base  processor  ' ’ 

baseV  r inPUt  ^ * P‘rt1CU,ar  Pr°9ra"-  The  *-  the  data 

• replace  data  base  co-ands  In  the  skeleton  Input  to  create  an  Input  file  that  the 

107  7 The  USt  15  ' ~ — — operates  with  the  data  base  to 

ZZ  X "IT  Ut‘er  ~ - "-^r  to  glue  greater  „ ^ 

flexibility  in  the  type  of  configurations  that  can  be  analyzed 

• s'.  * ..  . 


I 1 


• ! 


• M 
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ALAN  H.  HILHITE 


AVID  PROGRAM  INTERFACING  TECHNIQUES 


AERO 


PROGRAM 


t DIRECT  DATA 
TRANSFER 


SKELETON 


INPUT 


DATA  BASE 

S'  prograSK 

PROCESSOR 

— f INPUT  ) 

AERO 


PROGRAM 


• DATA  MANIPULATION 
WITH  DATA  BASE 
PROCESSOR 


CONVERSION 

A progra«K 

INTERFACE 

— ( INPUT 

- 

\J_J 

AERO 

PROGRAM 


• SEPARATE  CONVERSION 
PROGRAM 


AIM  W.  WILHITE 


* 


SUMMARY 


AVID  (S  USED  FOR  INTERACTIVE  VEHICLE  SYNTHESIS 


SEVERAL  AERODYNAMIC 
LEVELS  OF  DESIGN 


PROGRAM  OPTIONS  ARE  AVAILABLE  FOR  VARIOUS 


THE  AVID  SYSTEM,  WITH  ITS  INTERACTIVE 
DRASTICALLY  REDUCED  DESIGN  CYCLE  TIME 


GRAPHICS  SYSTEM,  HAS 

* i * . • , 


(=30  MINUTES  TO  DIGITIZE  VEHICLE) 

10-30  SEC  FOR  LEVEL  1 AERODYNAMICS) 
(~  * MINUTES  FOR  LEVEL  2 AERODYNAMICS) 


PROBLEM  AREAS 


• DIFFERENT  GEOMETRY  DEFINITION  FOR  EVERY  AERODYNAMIC  PROGRAM 

• SENSITIVITY  OF  RESULTS  TO  MODELING 

' SrTI0N  0F  SUBS0NIC  STABILITY  OF  VEHICLES  WITH  NONS LENDER 
BOD  I ES 

• GEOMETRIC  DEFINITION,  AERODYNAMIC  COMPUTATION  AND  AERODYNAMIC 

GRAPHICAL  OUTPUT  INTEGRATED  INTO  ONE  program  c 

• QUICK  ARBITRARY  HYPERSONIC  ANALYSIS  PROGRAM 
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RECOMMENDATIONS 

• STANDARDIZED  GEOMETRY  FOR  AERODYNAMIC  PROGRAMS 

• DEVELOP  INDEPENDENT  PROGRAMS  FOR  GEOMETRY  INPUT, 
AERODYNAMIC  COMPUTATION,  AND  GRAPHICAL  OUTPUT 

• DEVELOPMENT  SIMPLIFIED 

• MODIFICATIONS  ARE  EASIER 

• ALLOWS  MINI-COMPUTER  GRAPHICS  AND  HOST 
COMPUTER  COMPUTATIONS 
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ALAN  M.  WILHITE  - 


uMQINAL.  PAGE  IS 
OF  POOR  QUALITY 


PANEL  3; 
1. 


SYSTEMS  in  use 


2. 


3. 


i*. 

5. 


>S  there  3 need  for  be5te«’  display  of  surfaces? 

(« suxj'-  '— >» 

(c)  Color? 

•s  it  time  to  standardize? 

overlooked?  'tems  ,n  surface  representation  that  the  Workshop 


WORKSHOP  COMMENTS  & OBSERVATIONS  DURING  THE  SUMMARY  SESSION 


1.  STANDARDIZATION 
"Not  needed  (yet)." 

"Needed  when  there  Is  significant  information  exchange  (later?)". 

"List  needed  functions  instead". 

"Observe  CAM I and  ANSI  Standards  Committee  progress". 

2.  INFORMATION  EXCHANGE 

"Several  programs/systems  are  available,  e.g.,  QUICK,  ICAD,  iPEGS,  etc.  but  require 
dedication  to  adopt/install  elsewhere". 

"Parts  of  systems  are  available,  e.g.,  point  thinning,  patch  intersection,  etc.". 
"Minicomputers  and  graphics  terminals  are  an  .important  part  of  a system?'. 

"Advantages  and  limitations  of  various  methods  were  surveyed,  e.g.  Staley's  paper". 

"Items  missed  in  the  Workshop... "Fairing?*  methods". 

3.  UNDERSTANDING  THE  ISSUES 

"Conics,  cubics,  quintics  all  have  a place  (a  defender)". 

"Interpolation  with  polynomials  causes  spurious  waviness". 

"Non- Interpol at ivc  methods  are  an  alternative". 

"Wind  tunnel  model  builders  want  equation  representation  of  surfaces  not  coordinate  points". 
"Interpolation  with  polynomials  is  the  basis  for  most  current  systems". 

"Separate  geometry,  aerodynamics,  graphics". 

"Aerodynamics  can  be  sensitive  to  small  changes  in  the  surface  in  special  regions/flow 
conditions". 


-2- 


NASA'S  B0U  UL  SURFACE  REPRESFHTATinu 


t.;r'de  a S”a"  paclta9e  for  those  who  don't  have  their  own". 

AM  major  colonies  wll,  have  their  own  system  no  matter  what  NASA  does" 

:cZ7  t COmPaniCS  ^ 3 MSt  °f  *«»«■  (utilities,"  ‘ 

Create  aircraft  Reentry  hy  Poi„ts...the  «».  denominator". 

Need  to  communicate  slopes  and  curvature  too" 

"Need  to  know  the  difference  (tolerance)  between  one  form  „f 

De tween  one  form  of  representation  and  another' 


National  Aeronautics  and 
Space  Administration 

Aims  Research  Center 

Moffett  Field,  California 

94035 


Btply  to  Attn  ol:  YAE 1 22  7-2 


NASA 


December  29,  1977 


On  March  1-2,  1978,  NAS A- Ames  Research  Canter  will  host  a workshop  on 
Aircraft  Surface  Representation  for  Aerodynamic  Computation*  The  pur- 
pose of  the  workshop  is  to  exchange  information  on  converting  aircraft 
geometry  into  the  form  needed  by  aerodynamic  computation  programs. 
Primary  emphasis  will  be  on  aircraft  surface  specification  for  linear 
aerodynamics  paneling  programs,  but  sense  time  will  be  allotted  to  dis- 
cuss areas  of  commonality  with  aerodynamic  flow-field  mesh  generation 
and  possibly  with  computerized  lofting  systems. 

We  propose  to  include  presentations  in  three  topic  areas:  (1)  Geometry 
requirements  in  aerodynamic  computation,  (2)  Current  or  proposed  geom- 
etry methods,  and  (3)  Use  of  interactive  graphics.  The  presentations 
will  be  followed  by  panel  discussions  designed  to  explore  the  user's 
common  desires  and  concerns.  A summary  of  each  discussion  and  the 
visual  material  used  will  be  published  by  NASA. 


You  and/or  members  of  your  staff  are  invited  to  participate  in  the  work- 
shop by  presenting  material  and/or  entering  the  discussions.  If  you 
wish  to  attend,  please  address  all  proposed  presentation  material, 
comments  or  questions  by  January  23,  1978  to: 


Thomas  J.  Gregory  or 
Captain  John  Ashbaugh,  USAF 
NASA-Ames  Research  Center,  M.S.  227-2 
Moffett  Field,  CA  94035 
Telephone:  (415)  965-5881 


We  are  looking  forward  to  an  open  communication  of  ideas  at  the  workshop 
and  will  be  pleased  if  you  can  attend. 

rely  yours 


/UUi Bm 


C.  A.  Syvertson 
Acting  Director 


' * '&m.-  ' 

* ' National  Aeronautics  and 

{ Space  Adminiat  ration 

i nMtl  n^UBTOn 

' Moffett  FiaW,  California 

; > 94035 


NASA 


• ' { XwMyloAnnol  FAE:  227*2 

i 
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,! 
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I 
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I Otar  Attandae: 

j*“ 

Thank  you  for  your  interest  In  the  Ames  Workshop  on  “Aircraft  Surface 
Representation  for  Aerodynamic  Computation."  On  the  basis  of 
conversations  and  draft  material  from  each  presenter  we  have 
prepared  a preliminary  program  (Enclosure  I)  that  includes  both 
presentations  and  panel  sessions.  Possible  questions  or  topics  for 
the  panel  sessions  are  listed  In  Enclosure  2. 

• The  Workshop  will  be  held  In  the  Ames  Research  Center  Auditorium 

\J  (Building  N-201 ) (Enclosures  3 and  4)  and  start  promptly  at  8:30  a.m. 

Please  register  early  at  the  Visitor  Reception  Bui Iding  (N-253)  and 
, proceed  to  the  Workshop  parking  lot  as  shown  on  the  Enclosure  4.  The 

Auditorium  Is  a short  walk  (approximately  100  yds.).  Lunch  is  planned  ?" 

at  the  Ames  Cafeteria,  another  short  walk  of  approximately  200  yds.  . 

, , A list  of  area  motels  and  restaurants  Is  on  the  back  of  the  Enclosure  3. 


There  will  be  an  incoming  message  board  (telephone  (415)  965-5256) 
and  pay  telephone  available  In  the  Auditorium  Lobby. 


Thank  you  agatn  for  your  Interest  and  please  call  (415)  965-5681  if 
you  have  questions.- 


Thomas  /.  Gregory  /' 
Chief,  Aircraft  A&rod 
Branch 


Aerodynamics 


j ,t  ALi 

Capt.  John  Ashbauijh,  USAF 


t 


Enclosures: 

1.  Prel iminary  Program 

2.  Topics  for  Panels 

3.  Hap  of  Sunnyvale  & Ht.  View  Area 

4.  Ames  Research  Center  Map 


% 
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National  Aeronautics  and 
Spaca  Administration 

Ameeftaeearoh  Center 

Moffatt  FteW,  California 
94035 


NASA 


itapty  to  Mm  oi:  FAE; 227*2 


Osar  Prasantar: 

Thank  you  for  agreeing  to  praaant  matartal  at  tha  Amas  Workshop  on 
"Aircraft  Surfaca  Representation  for  Aarodynamlc  Computation."  On 
tha  basis  of  conversations  and  draft  material  suppllad  by  aach 
prasantar,  we  have  prepared  a preliminary  program  (Enclosure  I)  that 
identifies  the  prasantar,  his  organization  and  a title  or  topic  to 
be  emphasized.  White  much  of  the  Initial  matartal  supplied  was  vary 
broad  in  scope,  subsequent  discussion  suggested  that  tha  attendees  at 
this  Workshop  will  benefit  most  from  presentations  that  are  focused 
on  key  topics.  Therefore,  we  are  requesting  that  you  concentrate 
your  presentation  on  the  titles  listed  in  the  preliminary  agenda. 

Panel  sessions  are  planned  after  each  session  and  will  address  issues 
or  questions  that  appear  to  be  of  interest.  Enclosure  2 is  a list  of 
questions  which  you  may  want  to  consider  prior  to  participating  In  a 
panel  or  commenting  from  the  audience. 

in  preparing  your  presentation  please  keep  in  mind  that  a workshop 
environment  encourages  dialog  and  Interaction  between  the  audience,  the 
presenters  and  the  panels.  To  provide  adequate  time  for  the  panel 
sessions,  It  is  very  important  to  stay  within  presentation  time  limits 
so  that  we  derive  this  important  benefit  from  the  Workshop. 

Please  bring  one  xerox  copy  of  your  presentation  visual  material  and 
a companion  paragraph  for  each  vugraph  or  slide.  Please  insure  that 
the  last  name  of  the  presenter  and  a page  number  appears  on  each  page. 
NASA  will  provide  a copy  of  this  material  to  each  attendee  near  the 
end  of  the  Workshop.  These  may  be  mailed  to  attendees  after  the 
Workshop  If  printing  is  delayed. 

The  Workshop  will  be  held  in  the  Ames  Research  Center  Auditorium 
(Building  N-201)  (Enclosures  3 and  1») » and  start  promptly  at  8:30  a.m. 
Please  register  early  at  the  Visitor  Reception  Building  (N-253)  and 
proceed  to  the  Workshop  parking  lot  as  shown  on  the  Enclosure  k.  The 
Auditorium  is  a short  walk  (approximately  100  yds.).  Lunch  Is  planned 
at  the  Ames  Cafeteria,  another  short  walk  of  approximately  200  yds.  A 
list  of  area  motels  and  restaurants  Is  on  the  back  of  the  Enclosure  3. 


FAE:227-2 


2. 


There  will  be  an  Incoming  message  board  (telephone  (415)  965*5256) 
and  pay  telephone  available  In  the  Auditorium  Lobby. 


Thank  you  again  for  your  response  and  please  call  (415)  965*5681  if 
you  have  questions. 

Thomas  Gregory  J 
Chief,  Aircraft  Aerodynamics 
•ranch 


Enclosures: 

1.  Preliminary  Program 

2.  Topics  for  Panels 

3.  Hap  of  Sunnyvale  and  Ht.  View  Area 

4.  Ames  Research  Center  Map 
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NASA-Ames  Research  Center 


W. 

J. 

L. 

J. 

J. 

R. 


Arild  Bertel rud,  M.  S.  229-1,  Ph.  (415)  965-6156 
Reese  Sorenson,  M.  S.  202-1,  Ph.  (415)  965-6416 
Herbert  Hoy,  M.  S.  210-5,  Ph.  (415)  965-5320 
Joe  Steger,  M.  S.  202-1,  Ph.  (415)  965-6417 
Harve  Lomax,  M.  S.  202-1,  Ph.  (415)  965-5124 
F.  R.  Bailey,  M.  S.  202-1,  Ph.  (415)  965-6419 
P.  Jones,  M.  S.  233-10,  Ph.  (415)  965-6482 
P.  Mendoza,  M.  S.  227-2,  Ph.  (415)  965-6133 
King,  M.  S.  227-8,  Ph.  (415)  965-5856 
Rakich,  M.  S.  229-1,  Ph.  (415)  965-6192 
Murphy,  M.  S.  227-8,  Ph.  (415)  965-5655 
Carmichael,  M.  S.  227-2,  Ph.  (415)  965-6216 
Terry  Holst,  M.  S.  202-1,  Ph.  (415)  965-6032 

Tom  Pul  1 ian,  M.  S.  202-1,  Ph.  (415)  965-6417 

John  Viegas,  M.  S.  229-1,  Ph.  (415)  965-6200 

Steven  Mathews,  M.  S.  210-5,  Ph.  (415)  965-5943 
Vernon  J.  Rossow,  M.  S.  247-1,  Ph.  (415)  965-6681 
Chon-Yin  Tsai,  M.  S.  247-i,  Ph.  (415)  965-6681 
R.  Furey,  M.  S.  237-11,  Ph.  (415)  965-5888 
R.  Medan,  M.  S.  227-2,  Ph.  (415)  965*5880 
D.  Firth,  M.  S.  227-2,  Ph.  (415)  965-5855 
A.  Barlow,  M.  S.  227-2,  Ph.  (415)  965-5880 
Paul  Kutler,  M.  S.  202-1,  Ph.  (415)  965-5194 
S.  227-8,  Ph.  (415)  965-6396 
227-8,  Ph.  (415)  965-6116 
S.  237-9,  Ph.  (415)  965-5673 

S.  247-1,  Ph.  (415)  965-668! 

M.  S.  227-2,  Ph.  (415)  965*6216 
M.  S.  202-1,  Ph.  (415)  965-6032 


Roy  Presley, 

Gil  Chyu,  M. 
Jeff  Bowles, 
Larry  Olsen, 
Larry  Erickson, 
W.  F.  Ballhaus, 


M. 

S. 

M. 

M. 


C.  K.  Lombard,  M.  S.  202-1,  Ph.  (415)  965-6607 
S.  Ghose.  M.  S.  227-8,  Ph.  (415)  965-6396  „ 
Marcel  Vinokur,  M.  S.  202-1,  Ph.  (415)  965-6607 
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NASA-Lanqley  Research  Center.  Hampton,  VA.  23669 

James  C.  Townsend,  M.  S.  462,  Ph.  (804)  827-3181 
Carson  Yates,  M.  S.  340,  Ph.  (804)  827-2611 
Jim  Thomas,  M.  S.  280,  Ph.  (804)  827-3611 
Charles  Swanson,  M.  S.  280,  Ph.  (804)  327-2673 
Mary  Adams,  M.  S.  402,  Ph.  (804)  827-3181 
David  Miller,  M.  S.  402,  Ph.  (804)  827-3181 
C.  L.  W.  Edwards,  M.  S.  402,  Ph.  (804)  827-3181 
T.  R.  Rau,  M.  S.  36$,  Ph.  (804)  827-3911 
Carl  E.  Freeman,  M.  S.  286,  Ph.  (804)  827-3611 
Robert  E.  Smith,  M.  S.  125,  Ph.  (804)  827-2981 
A.  W.  Wilhite,  M.  S.  365,  Ph.  (804)  827-3911 


NASA-Lewis  Research  Center,  21000  Brookpark  Road,  Cleveland,  OH  44136 
Hsiao  Kao,  M.  S.  86-1,  Ph.  (216)  433-5520 


David  Taylor  Naval  Ship  and  Research  Development  Center.  Bethesda,  MD  20084 
Tsze  C.  Tai,  Code  1606,  Ph.  (202)  22J-1462 


Naval  Coastal  Systems  Lab..  Panama  City,  FL  32407 
Neill  S.  Smith,  Mail  Code  794,  Ph.  (904)  234-4348 


Airforce  Systems  Command  Lla.snn  Off.,.  M.,  R«earch  Fie,d>  „ 


Carl  W.  Tusch,  M.  s.  206-3,  (415)  965-5833 
D.  S.  McRae,  M.  S.  206-3,  (415)  965-5870 
John  Ashbaugh,  M.  S.  227-2,  (4 15)  965-5990 
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ATTENDEES 


~ S-  Afr  F°r?-e  Academy , Colorado  Springs,  CO. 
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